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AB55TRACT 


The purpose of thin renearch wan to inventiqate information extraction 
and tranamiasion techniouea for aynthetic aperture radar (SAP) imaqery. Four 
interrelated problems were addressed. An optimal tonal SAP imaqe classifica- 
tion algorithm was developed and evaluated. A data compression technique was 
developed for SAP imagery which is simple and provides a 5:1 compression with 
acceptable imaqe quality. An optimal textural edge detector was developed. 
Several SAP image enhancement algorithms have been proposed. A study was 
underta)cen to quantatively compare the effectiveness of each algorithm. 
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1 .0 INTHOIMK’TION 

Tli*« valijf* of rowot.** m*nninq havn U>fln cl*»arlv d»*inon- 

Htrat«>) h*/ th« nucc«SH of th« LANiiSAT of Hia tachnoloqy 

for axtractinq UHoful Information from tha data raturnad from thaaa satallltaa 
ia wall davalopad. Pecantly, a naw claaa of apacaborna imaqing ayatama hava 
haan taatad, Synthatic Aparturea Radar (SAR). Procaaainq the data from apaca- 
borna SAR ayatawa to form imaqea la a non-trlvial *aak and a qraat daal of 
raaaarch ia currently underway to develop high apead SAR proceaaora. Unfor- 
tunately, a aimilar effort ia not bainq devoted to the development of informa- 
tion extraction techniquea for SAR. The reaearch reported here developed 
information extraction and tranamiaaion techniquea for SAR. 

Four interrelated problems have been addresaed: 1) classification; 2) 

data compression; 3) textural modeling and edge detection; and 4) image en- 
hancement. An appendix is attached for each of these topics which describes 
the results of our investigations. 
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2.n H\rKnPMllNI) AND OVKWVIV'W 

ni'ilt'.il |>roc«>HHin<i t*»chni<iut*t h.tv»« ht»nn miroMSHf ul ly 

wid« v/iri»*ty of prohl*!iii ranginq from th»» «n.»ly«i»i ot X-r/iv ingiiHH to 
ficvition of handwriting to th** «^xtractior. of information from Bat*»llit»* imag»*« 
(1, 2, 3, 4, S). Kach nf^ problem in image procetjning present* new technical 
challenges which must he addressed. This research addressed im|X)rtant prob- 
lems involved in digital process!. tg of RAR images. 

Spaceborne imaging radars are important because they provide a unique 
view of the earth's surface. Their imaging geometry, spectral characteris- 
tics, and all-weather capability give active microwave systems ‘^n advantage 
over conventional imaging techniques, e.q., photography. For many years, 
imaging radars have been successfully used for military reconnaissan :e and 
geologic mapping. Operational systems arqcurrently in use for ice surveys and 
the detection of oil spills on the world's oceans. 

Until very recently, all informatin was extracted from radar images by 
human interpreters. Other techniques were nat considered because radar data 
was disseminated as photographic products, e.g., paper or film positives. Not 
only was it awkward for an interpreter to convert these data into a digital 
format for automatic or machine-aided information extraction, but this conver- 
sion process degraded the quality of the data. The advent of the digital 
signal correlator (e.g., 6) for synthetic aperture radars (SAP) has changed 
this, and now radar images are commonly distributed in a digital format. In 
addition, quantitative information is now desired from radar images, e.g., 
soil moisture and crop type estimates are being sought. The volume of data 
collected from proposed spaceborne imaging radar missions can be enormous, as 
shown by the Seasat-A SAP mission. If imaging radars are to meet their full 
potential, then automatic information extraction is needed, or at least ma- 
chine-aided analysis to ease the work load on the interpreters is required. 

An initial approach for automatic information extraction from radar 
imagery would he to employ the well developed technology associated with other 
sensors (1-51 (e.g., LANDSAT). This method was unsuccessful. The failure of 

We are concerned here with extracting information from RAP images, not 
.jrocessinq the received inphase and quadrature voltages to form the RAP image. 
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thin af>proflch occurred cxiiitinc] prons«Rin<i fllnorithnn w**r<* df*«iun**d 

Asnuminq ip«cific nnd HtAtiHticnl mod«»l. Th** moxt common model 

aRRumed waB additive, white riauRRian noiee (AWGN). Thin model doee not apply 
to radar imaqea. Therefore, procesRim alnorithma haaed on an AWCN did not 
f>erform aatiaf actor i ly when applied to ladar imaqeR. 

It has been shown that [7, R) radar imaqee can be modeled b>' a multipli- 
cative noise process which is non-Gaussian and has a poor sic,nal-to-noiae 
ratio. The standard technique for improvinq the interpretahi lity of SAR 
imaqes has been to use noncoherent inteqration 19-11). Noncoherent averaqinq 
is used extensively for noise reduction -n coherent systems [12-14] and can be 
implemented in irwiny ways. The basic idea behind this method is that indepen- 
dent samples (or looks) of the terrain are gathered and averaged after detec- 
tion. Independent samples can beobtained by polarization or frequency diver- 
sity. Most SAR systems use frequency diversity, i.e., nonover lappino sub- 
areas of the spectrum of the complex received signals are used to form the 
independent images. 

Continuous scanning of the spectrum is also used [11, 13, 14). The net 
effect is to reduce the bandwidth (degrade the spatial resolution), while 
improvinq the siqnal-to-noise ratio, S/N. Several studies [10, 15, 16) have 
shown the advantages of noncoherent processing for the interpretation of SAR 
imagery. However, there are several disadvantages to this approach: 1) the 

technique is spatially invariant and thus does not account for the multiplica- 
tive nature of the noise and the nonstationarity of the signal; 2) the tech- 
nique was developed for coherent optical processors and thus it is easily 
implemented with such a processor but is not necessarily optimum for digital 
processors; and 3) the technique is only aimed at improvinq the S/N and not 
for direct extraction of information. 

Several digital image enhancement algorithms have recently been developed 
[17, 18, 19) to treat the multiplicative nature of SAR images. Also, several 
hueristic approaches have been applied [20) . The technique we developed in 
[17) is an adaptive minimum mean square error spatial filter which preserves 
edges while smoothing homogeneous areas, e.g., agricultural fields. Homomor- 
phic filtering has also been used [18) for enhancement of speckled imaqes. A 
linear approximation to the multiplicative noise process is used in [19) to 
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df»v«*lon rtn .iil.iptiv*' 1 on fhi» lr>c<il n^an and varlanc«, Wh*T«»a*< (17) 

and M*>l dr*nonatrat»« t!i«* »>**rrornanr^ of th»»ir al<ior I fh»nfi on aotiinl RAR ima- 
»iarv, only a.mulatlona ar** us«*i1 in MR). All of alqorithwB on artual 

SAP iiBaq«rv, only Riimilationa art* ua«d in I IS). All of thana aloorithma hav«» 
only been diractod toward imatit* enhancamant aaauminq that the Information 
extraction procesa would he {>erformed manually. While thaae techniquen are of 
qreat value for manual interpretation of SAP imaqery, alqorithma for automat- 
inq the information extraction tanka are needad. 

Appendix A deacrihea a quantative evaluation of aeveral of theae enhance- 
ment algorithma. The comparison was based on both an edqe quality measure and 
computer execution time. 

What is the information to be extracted from spaceborne SAR imagery? 
Active microwave remote senainq is beinq used for several applications (e.q., 
qeoluqy, ice mapping, and moni torinq wave conditions in the ocean); several 
more applications have been proposed, e.q., soil moisture meastirement, crop 
classification, monitoring crop growth and harvest progress, and snow mapping 
( 211 . 


In each application, the SAR is used to measure different properties of 
the earth's surface. The problem of information extraction can be viewed as a 
signal analysis problem where different properties of the signal (here the SAR 
image) are used to imply certain properties of one target. For example, in 
geology SAR images are analyzed for their large scale spatial structure. 
Different image structures imply different geologic structures (221. The 
backscattered power (which is mapped into the SAR image intensity) is used to 
estimate the soil moisture of crop type (23). For ocean applications, it is 
the Fourier spectrum of the SAR image which is used to estimate ocean condi- 
tions. 

Thus, the successful extraction and handling of information for SAR 
images requires: 1) identifying unique image properties; 2) relating those 

image properties to the ratget characteristics of interest; 3) finding ways of 
efficiently measuring those image properties; and 4) finding ways of effi- 
ciently representing (for storage or transmission) these image properties 


App#'niiix r« roMtrtinx a <;ar f i on nlnorithm, Tn thi« ^Imo- 

rithm ir\rt«i»», ron»> w-tn fh*« proiwrtv of th** nlMnnl (RAW Imn-io) unert to oh^Ain 
information about tho ••arth, Thr mul t i pi irati v» for BAR imauan wax unwri 

to davalop a maximum liHalihojd clacaif ication algorithm. It wa« firat aa- 
Bumad that tha target feature information wa« known a-priori. A probability 
of incorrect claasif ication waa then determined for thia algorithm aa a func- 
tion of tl.e BAR paramatera (e.g., the degree of noncoherent averaging). A 
technique waa also developed to extract the target feature information from 
the supplied image. Thia classification algorithm was tested on RKASAT-A SAR 
imagery. 

Appendix C contains a description of a data compression developed for SAR 
images. This technique was based on the multiplicative noise SAR image model 
and is designed to preserve the local atatiatica of the image. The technique 
is an adaptive variable rate modification of the block truncation coding 
technique develo;>ed in (24). A data rate of approximately 1.6 bits/pixel is 
achieved with the technique while maintaining the image quality and cultural 
(point like) targets. The algorithm requires no large data storage and is 
computational ly simple. 

In some applicationa , the image attribute, desired cignal characteristic, 
which is needed is clearly defined aa in the case of image intensity. And tne 
classification algorithm described in Appendix C provides the maximum likeli- 
hood technique to separating image features baaed on intensity (or tone), 
however, for others, a distinct signal property has yet to be identified. One 
promising attribute is texture. Texture is known to be important for the 
manual interpretation of SAR images for geology and these large scale texture 
differences have been shown to be separable using digital techniques [221. 
Appendix D contains a description of an optimum textural edge detection fil- 
ter.* Thia filter will he of value in separating regions of different tex- 
tures in SAR images. The filter described in Appendix D is optimal in the 
sense that for the given texture model, a maximum amount of output signal 
energy is concentrated within a given resolution interval about the textural 
transition for a set filter bandwidth. The filter developed here is also a 

* 

The development of the textural edge detection filter was also partially 
supported by NASA Contract No. NASA 9-16664. 


qlohAl ofMirAtor, Otrtn a IoraI operator. In othnr wxjr'i*, th»* nntirR 

imAq* 1 a trA»4f onv»<l An<1 wortlfind hy th« filter, InAt^Art of hrAAklnq th«* 
procAHAinq <1<jwn into many IoraI operAtionn, Thw optimum tmxtiirAl ajIur dmtmc- 
tor ia an axtonaion of th« optimum tonAl wdqa d«t«ctor (2^1 • 

With tha incraasad availnhility of diqital RAR imaqaa, information ax- 
traction and tranamia.* Ion alqorithma will incraaaa in importanca. Thia ra- 
aaarch affort addraaaad aavaral critical problama in RAR imaqa procaaainq. It 
ia hopad that thaaa alqorithma will ha both uaaful and provida a baaia for 
furthar davalopmanat of radar imaqa procaaainq tachniquaa. 
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ABSTRACT 


Many algorithms exist for the enhancement of synthetic 
aperture radar (SAR) images. These algorithms improve the 
signsl-to>noise ratio of an image in order to make the 
image more useful for an observer. The trade-off for this 
improvement is a decrease in the image resolution and is 
seen as a blurring of fine detail. For many applications, 
most notably for the analysis of agricultural scenes, this 
blurring effect is only critical in edge areas where a 
boundary exists. 

In this study several different enhancement algorithms 
were implemented. These algorithms were then compared based 
on an edge quality testing procedure. This procedure 
established a method with which the edge-preserving 
abilities of the various algorithms could be compared. The 
amount of computer processing time required by each 
algorithm was also recorded. Using the results of these 
comparisons, recommendations are made as to which 
algorithms are best suited for various applications. 
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1.0 INTRODUCTION 


The usefulness of synthetic aperture radar (SAR) imagery 
is dependent on the ability of an observer to recognize 
detailed features in an image. This ability is often 
greatly decreased by the presence rf noise. Various 
algorithms have been developed to suppress noise in SAR 
imagery. (1,2, 3,4) The problem with noise reduction 
algorithms is that they tend to suppress the desired signal 
as well as the noise. This is particularly a problem for 
edge areas. The purpose of this paper is to evaluate 
several of the algorithms in order to determine how well 
each algorithm preserves edge information while suppressing 
noise. 

The principle type of noise in SAR imagery is speclcle 
noise, which is multiplicative in nature. Spec)cling effects 
are due to the fact that SAR generates images by the 
coherent processing of the reflected signals, resulting in 
more noise in those areas where the signal is greater. This 
can be modeled mathematically by the equation: 

Zi.j'Xij Vi,j (1.1, 


where is the observed power at a particular range and 
azimuth, Xi j is the signal that would ideally be observed 
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and Vi.,) is the noise. The subscripts are present in order 
to emphasize that the image is to be processed digitally 
wit'.) i and j corresponding to a row and column in the SAR 
image to be analyzed. Hereafter, for simplicity, the 
subscripts will be omitted. 

In order to improve the interpretability of an image, it 
is necessary to suppress speckle noise while enhancing the 
desired signal. Many different enhancement algorithms have 
been developed for this purpose. (1,2, 3, 4) Most algorithms 
suppress noise by averaging the surrounding points. That 
is, a pixel is replaced by an average of its neighbors, 
producing a smoother, less noisy image. However, since 
neighborhood averaging is often applied to all points 
within an image, the desired signal points are also 
averaged, leading to a degradation of the image resolution. 

In an edge area (defined as a boundary between two areas 
of differing average power return), the edge will also be 
smoothed, causing the boundary to. appear "blurry". In many 
cases tills retards interpretability as effectively as 
speckle noise. Obviously, it is desireable to suppress 
speckle while also preserving the edge information. 

A useful application of edge preservation is analysis 
of agricultural areas in which relatively large homogeneous 
areas of differing radar reflectivity have distinct 
boundaries. While smoothing in the uniform regions 
suppresses speckle, a loss of resolution in these areas is 
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rot <18 critical, since there should be little variation in 
signal strength for a level field containing only one type 
of crop. (Variables such as soil mosture or plant disease 
would cause the field to appear as more than one region of 
uniformity. ) However, no smoothing is desired in the edge 
area in order to preserve the distinct boundaries. The 
type of filter chosen for processing of noisy images 
significantly affects the amount of edge degradation in the 
processed image. 

The extent to which an enhancement algorithm preserves 
edges can be evaluated by an "edge figure of merit" (EFM) 
algorithm. (5) The EFM establishes a means of image 
comparison by detecting the amount by which an edge was 
smoothed. In addition to detecting smoothed edges, the EFM 
may detect "false" edges due to noise. These combined 
effects contribute to the calculation of a relative EFM 
value which is used in drawing conclusions about the 
effectiveness of each algorithm. 

, This report provides a quantitative comparison of 
several different noise suppression algorithms. Sectio i 2 
provides a brief description of the development of a 
digital noise filter. Section 3 explains the manner in 
which the filter performances were quantitatively compared, 
and section 4 outlines the development of each filter 
compared in this study. A discussion of the results is 
given in section 5. The appendices contain the listings 
for all computer programs used in this study. 
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2.0 DIGITAL SPATIAL FILTERING TECHNIQUES 

The purpose of this section is to briefly outline the 
considerations involved in the development of a spatial 
digital filter. Although many different types of image 
processing algorithms exist (1), this study focuses only on 
spatial filtering techniques. A large number and variety 
of algorithms have been developed to perform this type of 
filtering. The SAR image to be digitally enhanced is 
contained in a two-dimensional array of values representing 
the reflected power at each discrete area of the target. 
Each element of the array is a pixel in the image. Usually 
the image is been scaled to accomodate the display system. 

The term "image enhancement" is the process of 
suppressing the image noise while retaining the signal. To 
best illustrate the procedure used in developing a specific 
filter for a particular application, it is helpful to 
examine the development and implementation of a simple, yet 
effective, filter — the "equal-weighted" or "box 
filter. " 

Consider an NxM image f(i,j) containing both the 
signal and multiplicative noise. The enhancement procedure 
is to generate a smoothed image g(i,j) in which the gray 
level at every point (i,j) is the average of the gray 
levels of f contained in a predefined neighborhood of 
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(i,j). (For thif> study, the neighborhood of (i,j) includes 
the point (i,j).) That is, 


(n,m}£S 


f (n/fn) 


(?.i) 


where S is the set of coordinates of points in the 
neighborhood of the point (i,j), and T is the total number 
of points defined by the coordinates In S (7). 
Computationally this process involves calculating an 
average for each nxm region of the image. The nxm region 
is referred to as a "window" because of the way in which 
the entire image is viewed nxm pixels at a time. 

To generalize the filter, the filter window is defined 
in terms of an nxm array of weighted values (fig. 2.1). 

Each cell in the window contains a value which determines 
the degree to which the image gray level at that coordinate 
influences the average. By changing the window weight, it 
is possible to change the filter characteristics. 

Subsequent sections of this report examine the 
determination of window weightings for several different 
filters . 

As an example, consider a 5x5 window with a gaussian 
weighting to be used in filtering a 100x100 pixel image. 

By dividing each element in the window by the sum of the 
window elements, the window is normalized so that it sums 
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to one. The operation to be porformed is a "moving 
average.” Visualize placing the window on top of the 
image, starting with the upper left corner and moving 
across the image. At each position of the window, the sum 
of the products for the overlapping cells can be 
calculated. Since the window has been normalized, this sum 
will equal the weighted average of the values for the 
pixels "covered" by the window. The values in the window 
determine how much each covered image pixel is weighted in 
the average. For the normal weighting used in this 
example, the points farther from the center of the window 
have less influence on the average. Each sum is a point in 
the output image. Moving the window across generates a 
line in the output image. Once a line is output, the 
window is moved back to the left side cf the image, down a 
row, and then across to generate the next line of output 
(fig. 2.2 illustrates the operation). At this point it is 
clear that the window will not be able to cover enough 
points at the end of a column or row for determining an 
average (fig. 2.3). These end points can either be copied 
directly into the filtered image or discarded, making the 
filtered image smaller by an ^UT)ount equal to one less than 
the window size. (To prevent false edges from being 
created, the filters used in this study discard the end 
points ) . 

With a basic idea of why windows are used in image 
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processing, the computational requiremi*nts (or the 
implementation of a window operation can be examined. One 
of the most fundamental considerations for image processing 
is memory management. Typically a SAR image is quite 
large; usually there are over 500x500 pixels, while an 
image containing 2000x2000 pixels is not uncommon. To 
store an entire image in memory while it is being processed 
can take up many megabytes of storage space. Obviously 
some thought must be given to how much of the image is 
needed at any one time for processing. Most images are 
stored sequentially. That is, the top record or line of 
data is read first while following lines are read in the 
order that they appear in the image (fig. Z.4). Therefore 
the algorithm must also take into account the order in 
which the data is read. 

One of the most efficient algorithms to deal with 
window operations is a two-dimensional circular queue. 

This data structure is a " f irst-in-f irst-out ** construct. 
Since only the rows covered by the window are being 
processed at any one time, this is the only data which 
needs to be in memory. 

All of the programs presented in this paper utilize a 
circular queue to implement various types of window 
operations. A straightforward illustration of the queueing 
technique as applied to a spatially invariant 
two-dimensional convolution is presented in Appendix A. 
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This filter is better shown in the section on the 
equal-weighted filter, but this simple example program 
makes the queuing operation more obvious. In the next 
section of this report, the performance criteria for the 
digital filters are presented, and the development of the 
edge figure of merit algorithm is outlined. 
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3.0 Performance Criteria 


All of the filters coMpared in this study improve the 
signal-to-noise ratio of SAR imagery. However, the degree 
by which edges are degraded differs from one filter to the 
next. This section outlines the procedure used in 
developing a "fair" test of the edge preserving qualities 
of each filter. 

The most difficult question asked when developing a 
"fair" test deals with the amount of filtering done in the 
homogeneous areas of the image. This is an important 
consideration, since a filter might appear to perform much 
better based solely on the results of the edge figure of 
merit (EFM) test. The EFM algorithm used in this study is 
based on tne mean-square distance EFM developed by Pratt 
(S) defined as 


F = 


1 


h 

E 
1 • 1 


1 

l + ad2(l) 


(3.1) 


where Ij and represent the number of ideal and actual 
edge map points, a is a scaling constant, and d is the 
separation distance of an actual edge point normal to a 
line of ideal edge points. The scaling constant was chosen 
to be a *1/9 to provide a relative penalty between smeared 
edges and isolated, but offset, edges. 
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Since filter parameters differ widely from one 
algorithm to the next, a filter may not do as much 
smoothing in the homogeneous areas as another filter. If 
too little smoothing is performed, then the presence of 
noise will cause the detection of false edges, leading to a 
lower EFM. To validate the results of this study, each 
filter was applied so that the amount of filtering 
performed in homogeneous regions was approximately equal. 
Equal amounts of filtering were most important in the 
adaptive filters where a great deal of flexibility exists 
in the filter characteristics. For more rigidly defined 
filters, the window size is used as a common factor, while 
it still might be true that a filter operating on a 5 x 5 
local area of one type does as much filtering as a 7 x 7 
filter of another type. 

The input images used in this study were 144 rows by 
144 columns, containing only one vertical edge near the 
center of the image. The edge ratios used in this study 
were 3, 6, and 9 dB, with the number of independent looks 
equal to 1 and 4. A total of 6 images were processed by 
each filter for various window sizes. 

Once the filtered images were produced, they were 
processed using a differential operator (see App. A). For 
this study the 2x2 Robert's gradient was utilized. The 
gradient images became the input images for the edge figure 
of merit program which generates a binary image by 
thresholding the gradients. That is, if a gradient point 
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was found to be greater than or equal to an edge threshold, 
that point was set equal to a vaiue representing the 
existence of an edge. Otherwise the point was set equal to 
zero. This operation resulted in the generation of an 
"edge map" which, when displayed as an image, shows edges 
as bright areas and homogenecis areas as dark. Since the 
original images contained only one edge, the optimum result 
would be a single vertical line at the edge's location in 
the image. In practical cases, however, the actual edge is 
not completely identified. 

The threshold chosen in the generation of the edge 
map radically influences where edge points are identified, 
since too low a threshold highlights smaller gradients 
which might be false edges, while too high a threshold may 
miss some edges which actually exist. Therefore the edge 
figure of merit routine has two distinct, yet related 
tasks. The best threshold is sought in order to obtain the 
best edge map and the highest EFM, while the EFM is used as 
a way of determining the best threshold. 

The EFM routine implemented for this study uses a quadratic 
search for finding the best threshold in order to minimize 
the number of passes through the images, since each EFM 
calculation for a particular threshold requires processing 
the entire image. Once the best threshold was calculated 
for the 3 dB image, the same threshold was applied to the 6 
and 9 dB images. These calculations allowed comparison of 
the filter's ability to preserve low contrast edges. The 
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optimum threshold was also found for the 6 and 9 dB images 
so that further comparisons could be made. As an 
additional factor of comparison, the amount of computer 
processing time was recorded, yielding a relative measure 
of cost. The results are presented and discussed in 
Section 5. In the next section, the development of each 
filter is examined separately. 


4.0 FILTER IMPLEMENTATION 


In this section, each filter compared in this study is 
discussed separately. Complete program listings are given 
at the end of the report in Appendix A. The filters 
presented in this section were written, tested and 
implemented using a HARRIS 230 minicomputer and FORTRAN 77. 

Some routines were converted from existing programs, while 
others were designed from information supplied in reports. 
(Specific source references are made where applicable. ) In 
each case, an attempt was made to make the algorithms as 
general as possible with emphasis on readability. To 
achieve both generalizability and readability, the programs 
should be viewed as a package, since program layout and 
structure are common for most of the filters. 

Computational efficiency was an important consideration in 
the original design, and while specific algorithms could be 
optimized slightly, the common design approach does yield a 
usable product. 

4.1 The Equal-Weighted Filter 

The equal-weighted filter, often referred to as a box 
filter, is probably the most widely used spatial filter. 

The term "equal-weighted" describes the filter weighting 
function applied over each local area. In this filter, all 
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points in the area of the window are each weighted equally, 
resulting in the average for the area being the output 
point. The second term, "box filter," views the filter as 
a discrete two-dimensional convolution with a box function. 
This is given by the relation 

M- 1 N- 1 

f(U)*9(ij)»E E f (m^n) 9(x - niy y - n) (4.i) 

m ■ 0 n ■ 0 

for i«0 , 1 , 2 , . . . ,M-1 and j-0 , 1 , 2 , . . . ,N-1 where g(i,j) 
represents the SAR image as a two-di.mensional function and 
f(i,j) is the box function. The MxN array given by this 
equation is one period of the discrete, two-dimensional 
convolution (7). The equal-weighted filter performs a 
moving average over the entire image, smoothing edge points 
as well as homogeneous areas. This procedure makes the 
implementation of the filter straightforward, since only 
one pass through the filter window is required for each 
local neighborhood. It should be noted that this filter 
program is a straightforward application. The code can be 
optimized by moving some operations out of the innermost 
loop, resulting in a slight increase in efficiency. This 
optimization was not done here in order to maintain program 
readability. 

4.2 The Median Filter 
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The median filter is another commonly used filter, but 
in many respects it is quite different from other filters. 
The operation performed by the median filter is to replace 
an image point by the median value of the points in the 
local area. That is, for the points contained within the 
filter window, the point is found for which there are an 
equal number of pixels with lower intensities as there are 
pixels with higher intensities. 

In developing this filter, the question of window size 
must be addressed. Completely different programming 
approaches are required in order to optimize the filter for 
either small or large window sizes. This consideration is 
important, since the median filter, though simple in 
concept, requires a significant amount of computational 
overhead. For a small window size, it is easiest to sort 
the points and retreive the median. For larger window 
sizes, the increasing number of points makes this 
multi-pass approach very inefficient, and it becomes more 
desirable to generate a histogram of the points and sum 
counts until a median is found. The trade-off in central 
processing unit (CPU) time is summarized below for the 
relatively small 144 x 144 pixel images used in this study: 


WINDOW SIZE 

SORT METHOD 

HISTOGRAM METHOD 

3x3 

23.16 sec 

67.08 sec 

5x5 

90.64 

79. 14 

7x7 

271.37 

96.55 
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Larger images increase the need for computational 
efficiency. Since a more general approach was desired for 
this study, the histogram method was implemented. 

4.3 Lee's Edge Filter 

Most filters developed for the enhancement of SAR 
images take a general approach to the suppression of 
speckle; no special consideration is given to the filtering 
of edge areas. The local statistics algorithm developed by 
Lee (2) attempts to identify edges so that less smoothing 
can be done in these areas. 

To determine if an edge is present for a local area of 
the image defined by the filter window, the local 
statistics are first calculated. An edge is defined as a 
point of transition between two areas of differing 
properties. The property examined by the local statistics 
algorithm is the pixel intensity. If an edge exists within 
the local area, then a transition is present between 
relatively high and low pixel intensities. This results in 
the local area having a higher pixel variance. In this 
manner edges may be identified as being present within the 
local area by establishing a threshold for the variance. 

Lee also presented a statistical model based on the 
signal dependency of speckles (3). Using the model given 
by (Eq. 1.1), an extension may be made to find the a priori 
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mean and variance 


X ■ 



( 4 . 2 ) 


and 


Var(z) + 

Var(x)- vor(V)^v2 


where z and Var(z) are approxir.ated by the local mean and 
variance of the speclcle corrupted image. The parameters ){ 
and V are the means of the desired signal and the noise, 
while var(v) is the variance of the speckle noise. By 
linearizing the observed pixel z using the first-order 
Taylor series expansion about (X,v): 

Z»VX + X(V-V) (4.4) 


and 


Var(z) = E[ (XV - xv)^] 

= E(x^l EIv^] -x^v^. 


( 4 . 5 ) 


which can be further simplified if the window is assumed to 
cover an area of constant average intensity. Therefore, 
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EIx^l -x^ 


(4.6) 


and 


Var(z) - x^(El v^l - v^) 
-x^ Var(v) 

Var(z) 
Var(v) - — — 


Equations (4.5) through (4.8) are used to justify the 
multiplicative noise model, while also providing a 
simplified expression for Var(v) which is needed to 
complete Eq. (4.4). From Eq. (4.7) it can be further 
determined that 


Var(v) - Y 


(4.9) 


where N is the number of loo)cs for the SAR image. This 
information, along with the assumption that v>l and Eq. 
(4.2), leads to the final equation for the estimation of x, 

X»X + k(Z-X) (4.10) 


with 
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Vor(x) 

“ x2 Vor(v) > Var(x) 


This allows the actual signal to be estimated for a pixel 
by knowing the local mean and local variance for the filter 
window and also knowing che average number of looks for the 
SAR image. The best results were achieved by setting the 
number of looks parameter equal to an exact value 
calculated from the statistics for the SAR image being 
processed. The edge preserving quality of Lee's filter is 
achieved through his use of the local statistics of the 
homnaeneous areas of the image. If the local area is 
determined to be homogeneous based on the variance teut, 
then the statistics for the entire window area are applied 
to the model in obtaining the estimate for the signal. If 
an edge is present, then its orientation is found by 
applying a 3x3 gradient mask to a 3x3 array of subareas 
calculated from the window. Fig. 4.1 shows the steps taken 
in obtaining the edge orientation. Once the edge position 
and orientation are known, statistics are calculated for 
the homogeneous portion of the local area. The new 
statistics are calculated by masking the edge points in the 
window when calculating the new mean. These masks are 
given in fig. 4.2. These statistics are applied to the 
noise model to obtain a new estimate for the signal. Since 
this new estimate is only dependent on the homogeneous 
area, there is less degradation of the edge information. 
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Formation of 3x3 subarea meant from a 7x7 window. 



After a 3x3 gradient mask has been applied to determine the 
edge orientation, the pixels perpendicular to the edge are 
compared to one another to determine on which side of the edge 
the center pixel lies. 


Fig. 4.1 


Steps used in determining edge orientation. (2) 
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Fig. 4.2 Edge masks used in calculating new window 
statistics. (2) 








The entire procedure is flowcharted in fig. 4.3. 

For the purposes of this study, the implementation of 
Lee's filter was designed with an emphasis on flexibility. 
The window size was allowed to vary so that comparisons 
could be made based on this parameter. The edge threshold 
value, used as an input parameter for determining the 
presence of edges, is modified based on the results of 
previous executions of the program which report the 
percentage of the image assumed homogeneous. Since the 
images used in this study contain an exactly defined number 
of edge points, this allowed the results to be kept 
consistent for all applications of the filter. 
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4..4 The Adaptive Filter 


The operation of the adaptive filter is as the name 
implies; the filter is adapted to the local area based on 
some criteria. The filter developed in this study uses the 
local statistics of the area defined by the filter window 
in order to determine the filter window weightings. (1) 

For each local area of an image, the local number of looks 
is calculated using the relationship 


N 


Z 


Vor(z) 


(4.12) 


where Nj is the local number of looks, f is the mean of the 
area defined by the window, and var(z) is the variance of 
the local area. Using the value , an index into an array 
of filters is chosen. The filter weightings range from an 
equal-weighted filter to a filter which does no averaging. 
The filters defined within this range are weighted using 

F(x)=e’®'*' ( 4 . 13 ) 
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The rate of decay of the exponential determines how heavily 
surrounding pixels are weighted in the local average and, 
hence, how much smoothing is done for the pixel being 
processd. Figure 4.4 illustrates this for the 
one-dimensional discrete case, but an extension can easily 
be made to two dimensions. 

li. adaptive filter preserves edge information by 
applying a filter window which affects less averaging for 
the local areas with a lower Nj (higher Var(7,) ), while 
homogeneous areas receive more averaging, since a more 
uniformly weighted window is chosen for those areas. 

Several different parameters are involved in the 
definition of the adaptive filter, allowing a great deal of 
flexibility in the characteristics of the filter. The 
window size is variable, while the number of filters to be 
used is also variable. A larger number of filters allows a 
more continuous smoothing effect. To generate the filters, 
the first and last filters are generated, with the first 
filter equal-weighted, and the last filter weighted so that 
no averaging is done. The weightings of the filter windows 
within this range of filters were specified by the value, 
alpha, which relates to the rate of decay of the 
exponential. Alpha is first calculated, based on the 
relationship , 
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(4.M) 


a.l 

w 


where w is the equivalent resolution for a box filter. 
Alpha is then found by evaluating the integral 


W 


2 / 

0 


,-ax 


dx 


(4.15) 


Quantizing alpha gives 

a » K q Index (4.i6) 


where Kg is a constant evaluated fcr the case where a«.5 
when Nj-N, the number of looks for the image. The 
motivation behind these choices of values is based on the 
constraint that when the local variance (described by N^) 
is equal to the average of all the local variances 
(described by N), then the filter applied to this area 
should be the filter in the middle of the range of filters. 
This yields 

1/ _ ^ /A i-jy 

Aq ^ (4.17) 


with 


A 


N 


Max(Nt) -Mln(Ni) 
ff Filters 


(4.18) 


23 


Collecting terms produces 


2 Aw 

a - — Index 


( 4 . 19 ) 


which describes the filter shape for each filter. In order 
to generate a complete range of filter shapes, the program 
uses the user-supplied parameters, w and . The other 

filter characteristic which may be modified is the rate of 
filter usage. By calculating the local statistics, the 
local number of looks may be found from 

”«’v5F(T) 


which is used to select the filter to be applied to the 
area , 


Filter = ft Filters - 




( 4 . 21 ) 


Here the only parameter to be varied is , which also was 
used in the filter generation. The parameters, w and , 
are both factors used in the generation of the filter 
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is needed for the determination 


weightings, while only a 
of the filter usage. To accommodate this parameter 
dependency, the filter usage should first be determined. 
Once this characteristic is resolved, then the amount of 
averaging desired for the image may be regulated by varying 
the parameter w. 

4.5 The Edge Adaptive Filter 

This filter combines the edge-locating attributes of 
Lee's edge filter with the filter flexibility of the 
adaptive filter. Once Lee's edge filter determines the 
presence of an edge within the local area, the edge 
orientation can be found. With this information, it is 
possible to develop an adaptive filter algorithm which uses 
non-isotropic filter windows so that the edge pixels within 
the local area receive less averaging than those of the 
homogeneous portion within the window. 

Assume that for a particular local area of an image, 
a vertical edge has been found to exist on the left side of 
the window area (Fig. 4.5). Depending on the local 
statistics for the window region, the adaptive filter of 
Section 4.4 would apply an isotropic filter window to the 
area. However, since it is Icnown which pixels are part of 
the edge, a non-isotropic filter window may be defined 
which does less averaging for the edge points. An 
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Figure 4.5 A vertical edge is present on the left side 
of the local area. 
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exponential weighting function similar to the one used for 
the adaptive filter is utilized. However, the edge 
adaptive filter applies an exponential weighting function 
with a steeper decay on the pixels containing the edge 
(Pig. 4.6). This weighting is accomplished by first 
calculating a new set of statistics for the homogeneous 
porti ")f the window, establishing the weighting factor 
for the right-hand side of the window. For the remaining 
pixels within the window, filter weightings are applied 
using weightings for the filter defined by 

Fg ■ ^ Filters - F^^ + 1 (4.2?) 

where Fj| is the filter index chosen for the homogeneous 
portion of the local area, and F^ is the filter to be 
applied to the edge area. 

The parameters needed to define the edge adaptive 
filter are very similar to those used for Lee's edge filter 
and for the adaptive filter. 

4.6 The Sigma Filter 

The chief attraction of the sigma filter developed by 
Lee (6) is its simplicity and speed. The pixel to be 
processed is replaced by the average of those neighboring 
pixels having their gray level within two standard 
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Weight 



Figure 4.6 Non-Isotropic window weighting 


deviations from that of the concerned pixel. 

The sigma filter is baaed on tho multiplicative noise 
model for speckle in SAR images as presented in Eq. (1.1) 
with the assumption that the multiplicative noise, v<j , 
has a mean of 1 and a variance Var(v). From this it follows 
that 


2 ■ X V ■ X 


(4.23) 


and 

Var(z) » E( (XV - xv)^l 

- EIx^j EIv^I - x^v^ 


(4.24) 


For a small local area, the signal may be assumed nearly 
constant, allowing E(x2)»x^ ', which reduces Var(z) to 


Var(z) = x^ Var(v) 


(4.25) 


or 



(4.26) 


Equation 4.26 describes the standard deviation of 


27 


nmltiplicAtive noise in SAR imagery as the ratio of the 
standard deviation of z and the mean of z. 

It is assumed that z^ j is the a priori mean of 
and also that pixels in the window with gray level within 
two standard deviations from Z(j are from the same 
distribution. Since the speckle noise is multiplicative in 
nature, the two-sigma intensity range from Eq. (4.26) is 
(zi,j-20v z<j, z<j + 2oy The average of pixels in 

this intensity range replaces the center pixel as the 
smoothed value of zij (6). 

The implementation of the sigma filter requires that 
the relative limits for the distribution of Oy be 
established at the onset of program execution. In the 
filter implemented by Lee, a normal distribution was 
assumed. For comparison, a chi-square distribution was 
also modeled in this study. Additionally, a threshold K is 
established to deal with the presence of spot noise in the 
filtered image. The retention of spot noise is due to the 
fact that the gray level of the spot noise is significantly 
different from its neighborhood pixels. If the total 
number of pixels within the two-sigma range is less than or 
equal to K, the center pixel is replaced by the average of 
its four neighbors. 

The sigma filter was applied to the test image usina 
a wide range of configurations. Both chi-square and normal 
distributions were modeled, while the threshold K was set 


28 


to 3 for the 7x7 window and 2 for the 5x5. Reaulta are 
preaented for a aingle paaa and for a aecond paaa of the 


5.0 RESULTS, CONCLUSIONS AND RECOMMENDATIONS 


This final section presents the results of the filter 
performance comparisons. The results are shown in plots so 
that comparisons can be more easily made. The edge figure 
of merit value is expressed as a percentage shown on the 
vertical axis. The EFM is plotted versus the edge contrast 
ratio for varying filter window sizes and number of looks. 
The EFM is also shown with respect to the window size for 
varying edge ratios and number of looks. These results are 
summarized in fig. 5.25 which ranks the six filters using 
several different criteria. A value of one indicates the 
best performance by a filter. 

There were two methods used in obtaining the EFM 
values. The first method calculated the best edge 
threshold for the 3 dB image and applied the same threshold 
to the 6 dB and 9 dB images to determine their EFM. The 
second method found the optimum threshold at each edge 
ratio and used this threshold in calculating the EFM. 
Results are plotted for both methods. 

Figures 5.23 and 5.24 show the central processing 
unit (CPU) time required of each algorithm. This is first 
shown in terms of the actual time (fig. 5.23) reported for 
each algorithm using various window sizes. The values are 
then scaled and shown in fig. 5.24 as a relative time 
versus window size. The relative time factor is determined 
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by dividing the time required for each algorithm by the 
time required for the equal-weighted filter of the same 
size. The scaling factor varied with window size. 

The results should be analyzed carefully, since the 
wide diversity of filtering methods leads to a large number 
of factors which should be considered in the evaluation of 
each filter. Several general observations can be made 
about the results. As expected, the higher contrast edges 
y< u a higher edge figure of merit. Also, the four-look 
images, since they are not as noisy, show a higher EFM. 
These variations are independent of the filter type. 
However, allowing for this inherent EFM improvement, it is 
still possible to draw some corclusions about how well a 
particular filter type performs relative to edge quality. 
Another general observation is that all filters yielded 
results which were, to varying degrees, better than those 
obtained from *-he unfiltered images. As a final overall 
observation, it should be noted that for most filters an 
increase in the filter window size led to an increase in 
the edge figure of merit. This improvement can be 
explained by the fact that, since a larger window usually 
allows more averaging, fewer false edges are detected. For 
reasonably sized windows, the suppression of false edges is 
usually enough to outweigh the contrast loss caused by a 
larger window performing more smoothing. 
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5.1 Evaluation of Equal-Weighted Filter Performance 

The results for the equal-weighted filter are 
important, since it is the most widely used spatial filter 
and may be used as a standard for comparison. The 
popularity of the equal-weighted filter is justified by the 
results. The filter, even though it makes no assumptions 
concerning edges, averaging all pixels equally, still 
yields a fairly good EFM when compared to the other 
algorithms. This performance is clearly evident from the 
results. For the low contrast edges, the equal-weighted 
filter provides results which are only a few percent higher 
than those of the other filters. (A percentage value given 
for comparison is a value difference read from the 
graph — not a ratio of the EFM for each filter. ) For 6 dB 
edges, however, the equal-weighted filter is superior by as 
much as 20 percent. As the edge sharpness is increased to 9 
dB, the results of the other algorithms improve 
dramatically, though the equal-weighted filter still earns 
a comparable EFM. 

A comparison based on filter window size also shows 
the effectiveness of the equal-weighted filter. A 5x5 
pixel window gives an EFM very nearly equal to that for the 
more complex filters, while the 7x7 window gives results 
which are as much as 20% higher than those of the other 
filters (fig. 5.2). A 9x9 equal-weighted window was not 
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applied in this study, so a definite comparison cannot be 
made for large window sizes. 

Fig. 5.25 shows that the equal-weighted filter 
performed best for the N>1 images. This conclusion is not 
as evident from the plots, since fig. 5.25 was based on an 
average for all the N»1 data. The fact that the 
equal-weighted filter obtained top ranking in this catagory 
is not surprising when it is considered that in order to 
prevent false edged from being detected, more averaging is 
required. The equal-weighted filter performs more 
averaging than any other type of filter using the same 
window size. 

As a further benefit, when processing cost is 
considered a factor, the equal-weighted filter algorithm 
turns in the best times of all the filters considered. In 
some cases, the algorithm is over three times faster. This 
filter has another advantage in that it is very easy to 
apply, since the only parameter is the size of the window. 
This can sometimes be an important consideration when 
compared to the more complex edge filters utilizing several 
different parameters. 

5.0 Evaluation of Median Filter Performance 

The median filter is also quite common and simple to 
use. However, the results of the EFM were nearly the 
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lowest of those from all the filters. This is evident from 
the plots which show that the median filter gives better 
results for larger window sizes and greater edge ratios. 
This improvement with window size and edge ratio is 
probably related to the nature of the EFM as discussed 
earlier in this section, since the improvement is not 
nearly so dramatic as for the other filters. Fig. 5.25 
shows that in every category compared, the median filter 
e,:;rn8 a low ranking. In addition, the large overhead in 
processor time tends to discount this filter as a practical 
alternci..lve to the other filters. 

5.3 Evaluation of Lee's Filter Performance 

Lee's edge filter compares well to the other filters. 
This was expected, since special processing is performed 
for edge areas. However, for' smaller windows the algorithm 
does not do as well as the box filter. In fact the lowest 
overall EFM values are for Lee's filter using a 5 x 5 
window (fig. 5.1 and 5.25). This poor performance is 
probably due to the fact that with a smaller window, there 
is a greater chance that a false edge will be detected by 
Lee's filter. This speculation receives some justification 
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when it is noted that for the less noisy N«4 images and for 
more distinct edges, the equal-weighted filter is not as 
significantly superior to Lee's filter. 

As the window size, edge ratio or number of looks was 
increase', the performance of Lee's algorithm improves 
dramatically as shown in the plots. In some cases, Lee's 
filter had an EFM as much as 40 percent greater than that 
of the other filters (fig. 5.5). A conclusion that may be 
made about this filter is that it is prone to the detection 
of false edges, but by using less noisy images (larger N) 
with higher contrast edges and a larger window size, quite 
good results can be attained. 

The time for processing Lee's edge filter might bp 
considered too long — especially since a larger window size 
and more CPU time is required in order to achieve the best 
results of the filter. It should also be noted that the 
CPU time shown for Lee's filter is somewhat misleading, 
since the parameters for Lee's filter were adjusted for the 
mostly homogeneous test images. These images are not very 
realistic, since most SAR imagery contains more edge 
information. Between 3% and 7% of each test image was 
known to contain edges while roughly 20% to 60% of an 
actual image may have edge information. For actual 
imagery, the required CPU time is much closer to that for 
the more time consuming filters compared in this study. 
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5.4 Evaluation of Adaptive Filter Performance 

The adaptive filter also yields a high EFM, but unlike 
Lee's filter, achieves good results for the smaller window 
sizes. The adaptive filter earned the best ranking for low 
contrast edges (fig. 5.25), and when the rankings for all 
the comparison criteria are averaged, the result is 2.75 
(The result is 2.0 if the computer time is not included in 
t.he av *rage. ) . 

Again, noisy images and low contrast edges tend to 
lead to inferior performance, since the adaptive filter 
also uses the local statistics method in order to determine 
the amount of smoothing to apply to an area. Larger window 
sizes improve the performance of the adaptive filter for 
these images. 

The amount of CPU time 'required for the adaptive 
filter is quite large. This is especially true when 
compared to that of the equal-weighted filter, but also 
true when compared to the time required by Lee's filter. 
However, since Lee's algorithm requires a larger window 
size to get similar results, this disadvantage is slightly 
offset. For the adaptive filter, the processing time 
required for larger windows becomes even harder to justify, 
but some of the best results were achieved using this 
filter configuration (fig. 5.9). 
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A disadvantaqe of the adaptive filter is that it is 
relatively difficult to establish the filter 
characteristics. This disadvantage can be minimized by 
gaining experience with the filter, but often it can still 
take at least one initial pass with the filter before 
optimum filter characteristics may be determined. 

5.5 Evaluation of Edge Adaptive Filter Performance 

The edge adaptive filter gave some of the best results 
of all filters compared. For the N«4 images with a 5x5 
window, this filter earns the best overall EFM (fig. 

5.10). The improvement is slight, and in general the 
unmodified adaptive filter produces better results. As is 
shown in figure 5.25, the edge adaptive filter earned the 
best ranking for small window sizes, but the filter also 
obtained the worst ranking for noisy images. CPU time is 
not significantly higher than that for the adaptive filter, 
while the smaller window size allows good results without 
using a great deal of processor time. 

5.6 Evaluation of Sigma Filter Performance 

The results for the sigma filter are somewhat 
discouraging. This is shown in figure 5.25. Even though 
the filter is quite fast, two passes are required in order 
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to receive results comparable to those of the other 
filters. Two passes more than doubles the amount of 
processing required, yet the results are still inferior. 

In fairness, it should be noted that the poor results are 
probably due to the presence of spot noise which the EFM 
routine detects as false edges. Lee gives several 
techniques for reducing spot noise, which if applied, might 
lead to better results for the EFM test. 

5.7 Recommendations 

Having reviewed each filter separately, several 
conclusions and recommendations can be made. Figure 5.25 
cannot by completely relied on, since it was compiled from 
averages for all values given for each particular 
parameter. Optimum combinations of parameters are not 
shown by these rankings. In addition, rankings were not 
calculated for large window sizes. 

If processing time is not a factor, then it can be 
assumed that a large filter window should be used, since a 
higher EFM results by using the larger windows. Lee's 
algorithm gives the best results for the noisier images, 
while the adaptive or the edge adaptive would be better for 
less noisy images. If a smaller window size is preferred, 
then the equal-weighted filter should be used for noisy 
images. If processing time is an important consideration. 
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then the equal-weighted filter is clearly the best choice. 

The main basis of comparison in this study was the 
edge figure of merit. Though the desirability of an 
algorithm should not be determined solely from this 
criterion, some indication is given as to how well each 
algorithm will filter edge ar:;as. 
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FILTER PERFORMANCE RANKINGS 


Technique 

Low 

Contrast 

High 

Noise 

Small 

Window 

Computer 

Time 

1. Equal Weighted 

3 

1 

3 

1 

2. Median 

6 

5 

4 

4 

3 . Lee ' s 

4 

2 

6 

3 

4. Adaptive 

1 

3 

2 

5 

5. Edge Adaptive 

2 

6 

1 

6 

6. Sigma 

5 

4 

5 

2 


NOTE; The rankings for each colunm were determined by 
averaging the results for that parameter. The results for 
optimized thresholds were not included in the averages. 


Fig. 5.25 
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APPENDIX A 


Listings for computer programs 

Program Description 

Name 

CONVO Example of an equal-weighted filter 

DIFOP Discrete differentiation of an image 

EFM Calculate the edge figure of merit 

EQUFLT Equal-weighted filter 

MEDFLT Median filter 

LEEFLT Lee's edge filter 

ADPFLT Adaptive filter 

EADFLT Adaptive filter with a non-isotropic filter window 

SIGFLT Sigma filter 


All required subroutines are listed after each mainline. 
However, modules that are common to several routines are given 
at the end of the appendix. 


C 'PHIS IS A PROGKAM TO PKRKOHM A TWO-DI MKNS lONAL SPATIAL 

C CONVOLUTION ON A REAL ARRAY. 

C 

C 

SUBROUTINE CONVO (OUEUK, WINDOW, OUT, WNDSIZ, SIZE, 
fc DUMMY 1, DUMMY 2, OUTSIZ) 

C 

C THE DUMMY ARGUMENTS ARE USED IN THIS IMPLEMENTATION IN 

C ORDER TO ALLOW VARIABLE ARRAY DIMENSIONS. 

C 

C 

IMPLICIT INTEGER (A-2) 

BEAL TOTAL, QUEUE (WNDSIZ, SIZE), WINDOW (WNDSIZ , WNDSIZ ) 
KEAL OUT (OUTSIZ) 

C 

C • 

C INITIALIZE CIRCULAR QUEUE WHICH WILL STORE A HORIZONTAL 

C STRIP OF THE INPUT IMAGE. 

C 

DO 10 X - 1, WNDSIZ 

READ (1) (QUEUE (X, WRD) , WRD - 1, SiZE) 

10 CONTINUE 

C 

C 

C INITIALIZE RECORD COUNT AND QUEUE POINTER. 

C 

REC « 0 
QREC - 1 
C 

C BEGINNING OF OUTERMOST LOOP. SET THE TEMPORARY QUEUE 

C POINTER EQUAL TO THE FRONT OF THE QUEUE. 

C 

20 TMPQRC » QREC 

C 

C PRCXESS A ROW OF THE IMAGE 

C 

DO 60 START - 1, OUTSIZ 
C 

C PROCESS THE CONTENTS OF THE WINDOW. INITIALIZE THE SUM 

C OF THE WINDOW PRODUCTS TO ZERO. 

C 

TOTAL - 0.0 

DO 80 WREC » 1, WNDSIZ 
WWORD - 1 

DO 50 QWORD - START, START + WNDSIZ - 1 
TOTAL - TOTAL + QUEUE (TMPQRC, QWORD) 

& * WINDOW (WREC, WWORD) 

WWORD ■ WWORD + 1 
50 CONTINUE 

C 

C UPDATE THE QUEUE POINTER 

C 

TMPQRC « MOD ( TMPQRC, WNDSIZ ) + 1 
C 

80 CONTINUE 


J) 


ooo o o o o n non 


C 

C 

60 


C 

200 

C 




OUT (STArtT) - TOTAL 
CONTINUE 

WRITE THE OUTPUT RECORD AND READ IN A NEW ONE 
WRITE (2) (OUT (WRD), WRD - 1, OUTSIZ) 

REC ■ REC + 1 

READ (1, END-200) (QUEUE (QRLC, WRD), WRD - 1, SIZE) 
UPDATE THE POINTER TO THE FRONT OF THE QUEUE 
QREC - MOD ( QREC, WNDSIZ) + 1 

LOOP BACK FOR ANOTHER RECORD UNTIL AN EOF IS REACHED 
GOTO 20 

PRINT,'* ** ALL DONE***' 

PRINT, 'THERE WERE', REC, 'RECORDS WRITTEN' 

STOP 

END 
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CONTINUE 
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REC « REC + 1 
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1 F ( ( :l T T Yf»r; .Nt;. 1) .A7.C. (CLTTYPF .K". ^ ) ) THf. 

WR I T ? ( T r YCl.T , ) 

62C FC?KA r ( IX , • EP'xCK - - - CUTRLT IKACF TYi»t A?c 1 ANC ) 

G C T C I C 1 C 
£^G IF 

I? C.-LTSIZ » SIZE - i.NDSl? ♦ 1 
ARRSIZ = /..NCSL’ UNOGIZ 
REGS » n 
C 
C 

r. CAU SfHSCLTINE TC DC Thj] -..CRK 

f. 

call CIFSUei INFrH, aUTFf.y, ARRSIZ. bltn, RFCD, 

^ £ hNCSIZ. CUTSIZt TIYCUT. ILF'JF, Ur,DCN, T)uF, CP ►CM . C L T T Y F p ) 

C. 

C 

VsPITE {TTYCL'T, G31) 

631 FCP. FAT (IV, '<•>>> ALL CCNt »«><>•) 

► RITE ( TTYLLT .t.AO) flUTSlZ, REGS 

6tC FCKf-AT (lX.'TFt CUTPUT II«AGE IS 'tl',' ‘/C’CS 'Y ' , 1 *; , • RLCf.PCS’) 

C 

ICIC STLP 
ENC 

ECF . . 
o 


ORIGINAL PAGE IS 
DE POOR QUALITY 


o o n o 


C lMVLHSITY CF K/ikSA*j TFLECCPf'LM F at ICNS ANH I N F C? ^ A T I C N ' CFS LAB 


FPCCRAK suit : NflSE FILTLPS REF. « : 


poCGRAH KAF- C : C I FSL> AU7MC:^:JtFF V.A1'.C\ LATF;/./ ‘•/el 


c 

c 

c 

c 

r 

C 

C 

PURFCSt : 
FRIER AFTER 

CCFS THE PROCESS 
3LINC CALLEO BY 

INC FCR T(- 
0 I F C P . 

r C 1 FFER''-NT I AL CPEPA T CR 

fikjUlNAL PAGE IS 
DF POOR QUALITY 

L 

c 


PARAHE TER 

') EF IN 1 T 1 Ci\ 


c 

NAPE 

\ TYPE 

\ r.LASSV RANGE 

\ C'SCRIPTION 

L 

c 

I NF ILE 

VCHsMrJ 

\R 

\ 

\ INPUT filename 

c 

CUTFILE 

\Ch«13 

\R 

\ 

\CUTPUT FILENAME 

c 

ARRSIZE 

\ 1 

\R 

\ 

\NUKC.EIG OF elements If wINCCW 

c 

NCCLS 

\ I 

\R 

\ 

VCCLUKUS I'J INPUT IMAGE 

c 

NCUTRCW 

S 1 


\ 

\RCh'S IN OUTPUT IMAGE 

c 

W INCS I ZE 

\ 1 

\R 

\ 

WvINOCWSlZE 

c 

NCUTCCL 

\ I 

\R 

\ 

\CCLUMNS IN OUTPUT IMAGE 

c 

TTYCLT 

\ 1 

\R 

\ 

\cutplt tc terminal fileccoe 

c 

C 

\ I 

\R 

\ 

\C1RCLLAR CUfUE 

c 

APR VAL 

\ I 

\.^ 

\ 

\Nli\''CW ^RRAY 

c 

TeuF 

\ 1 

\U 

S 

\REAO IN aUFFEr* 

c 

CPNlf* 

\ I 

\R 

\ 

\CPERATCR NUMi’.cr 

c 

CLTTYPE 

\ I 

\R 

\ 

\Ai\GLE CR GRADIENT CUTPUT 

c 

r 


\ 

\ 

\ 

\ 

U 

c 

r 


NCN 

-LOCAL 

VARIABLES 


L 

c 


\ 

\ ■ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 

r 


\ 

\ 

\ 

\ 







c 


SUBROUTINES RECUIREC 


c 

r 

NAPE 

\ 

OESCR IPTICN 


c 

G6NP5K 

\GENERATES OIFF. CP MASKS 

c 

GETPNT 

\F INCS 

PCINT VALUE FOR CUTPOT 

c 

CPN 

\CPENS 

FILE AND ASSIGNS 

LCGICAI UNIT 

c 

UCLCSE 

\CLQSES 

F 1 1. E S 

OPENED W ITH CPN 

c 


\ 




c 


\ 



* 

c 

c 


\ 





C 

c 

SLOROUTINE CIFSU3 ( INFlLEt CUTFILE. ARRSIZE. NCOLSt N3LTPrw, 
i WINCSlZEt NCL'TCCL, TTYOUT, Q, ARRVALt TBUF. CPNUf'. CLITYPE) 
C 

LCGICAL ERR 
C 

INTEGER APRSIZEt NCOLSt NCUTCCLt NINCSIZE, CPNUR. CLTTYPE 
INTEGER It Ji Kt ERRNUKt STARTt ARRPCSt C'JTCCL, NPQRS 
INTEGER NUUTRCVtt TTYOUT t A RR V A L ( A R R S I Z t ) , TRLF(NCCLS) 


. L L r 


I 

1 




c 



I 



I 

I 



c 

c 


c 

c 

c 




c 

c 

c 

c 




I 


c 

c 

c 

c 




c 

c 

c 


IM6GER .. S I t ,Cc;i M t IShNu^. 

INTEG^>» MSK.M*}*. Gl, G2 

IMfCCR*--! lNFIl.6(UJ)t OLTFILt(lS) 

REAL ANGLE . I' I 

NCtlRCs = 0 
PI = <1 0 A T A N ( I . C ) 

CALL GE.\‘'SK( MSKi, MSK2t ARRSIZL. CPMjV) 

CALL CPN IINFNUK. INFILF. ’CLO’t 'CNF', EPRNU''. f.RR) 

I F (ERR » CCTC <<2 

CALL CPN (CfTNU»'. CUTFILCi ’NEw'f ‘UNF’. ERRNUPt ERR) 
IF (ERR) GCTC TA 

INITIALIZE GUF.UL 

DC 30 I = l.KINGSIZE 

REAC( INFNLP.ENC = 93) (T3UF(K), K = I, NCCLS) 

CC 20 J = l.NCOLS 
CM I . J ) » T6UF ( J > 

20 CONTINUE 
30 CCNTINL'E 


MAIN PRUCESSUG 

AO Continue 

CC 70 START = 1, NCUTCCL 
CC 60 I » 1 f WINOSI ZE 
DO 50 J - START, START ♦ WIKDSIZF. - I 
ARRPCS = W I NOS I Zt«( I - I ) ♦ J-START+1 
ARRVAL ( ARRPOS ) = C( 1 , J ) 

50 CONTINUE 

60 CONTINUE 

CALL GETPNTI KSKl, ARRSIZE, ARRVAL, Gl) 
CALL GETPNTI MSK2, ARRSIZE, ARRVAL, G2) 

IFICUTTYPE .EC. 1) THEN 
VALUE » SORT IREAL(G10#2) ♦ REAL ( GZ^'JZ ) ) 


ELSE 
IF ( Gl 

.EU. 0 

) 

THEN 

• 

IF( 

C2 .GT. 

0 ) 

ANGLE 

= P I /2 .0 

IF( 

ELSE 

G2 .LE. 

0 ) 

ANGLE 

a 3 <« PI / 2.0 


ANGLE a ATAN( REAL(G2)/ REAL(Gl) ) 
ENC IF 


ANGLE IS NCH BETWEEN -PI/2 AND PI/2, AND WE WANT IT BETWEEN 
C ANC 2PI. MLST CHEC); SIGN OF Gl ANO G2. 

IF( Gl .LT. 0) ANGLE * ANGLE ♦ PI 

IF(( G2 .LT. 0 ) .AND. ( Gl .GT. 0 )) ANGLE = ANGLE ♦ 2«PI 
IF( CPNUM .EO. I ) THEN 
ANGLE = P I /<» .0 ♦ ANGLE 
ENC IF 

VALUE = 255 « ANGLE / (2.C « PI) ♦ O.T 

FCR ROBERTS, T»HE ANGLES WERE RCTATEC BY PI/A . THIS COULC 

CAUSE PRCOLEMS WHEN THE ANGLE GETS ROTATED INTO THE FIRST GUADRANT 


LU 1 N*.i ' . 


OaHfTKE PAGE IS 
or. POOR QUALITY 


nnnoooo nr*o 


IF( value .G1. ) VALLt » VALLf - .1i‘> 

ENC 1 F 

CUTCCl = start 
TBUF(CUTCCL) = VALUE 
70 CONTINUE 
C 

NCUTitCK » NCLTROW ♦ 1 

WR I TE (CUT NUK ) (TiJUF(<0, K » I, NOUTCCL) 

UP0AT5 CUEUE 

cc ao I » i , wiNCsizE 
IF(1 ,LT. wINCSUE) then 
cc 85 J » 1. NCCL5 
C( 1 t J ) s 0( I M . J ) 

85 CONTINUE 
ELSE 

REAO( INFNUM.END = 98) (TBUFIK)# K » 1. NCOLS) 
CC 87 J = l.NCOLS 
C(I.J) « T3UFU) 

87 CONTINUE 

END IF 
80 CCNTINUE 
GC TC 

UE HAVE REACHED THE END CF THE INPUT PILE 
92 CALL FILERR (TTYOUT, INFILE. ERRNIJR) 

GCTC 99 

9A CALL FILERR (TTYOUT. DUTFILE. ERRNU>^) 

GCTC 99 


98 CCNTINUE 

CALL UCLCSE C INFNUM ) 
CALL UCLCSe(CUTNUM) 

99 CCNTINUE 
RETURN 
END 






ECF 

7 


w 

% 


c 

c 

r, 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


LNIVCRSITV CF KANSAS T E L E C CM K L M C A T I C N S ANi) I M i. S lAd 


PncCAAP SUTt : .\CISF FILTERS SFF. « : 


P9CCRAH NA^^E :CENRSK AI) THCJT : JEF F V.ATSC'I CAT*’ T 


PLP.PCSF. ; TC GtNLnATE THl lli.TfiF.HTS. >*rtfck.ITT, CP. SC'»cL 

ciffepfntial CPCPATOR RASK. 


PARAi'ETF.P TGFPJmCN 


NAPE 

\ 

TYPE \ f.LASS\ range 

\ CFSCRIPTluN 

PSK 1 

\ 1 

\W 

\ . 

NCIFFERINTIAL OPERAT'R PAS< 

MSK2 

\ I 

\A 

\ 

NCI FFE»EfniAL OPERA', SR PASK 

ARRSI ZE 

\ I 

\rt 

\ 

NMJPBEP. CF lilNCOV. ELEPENTS 

CPNLP 

\ I 

\R 

\ 

\rPERATCR NUPBFR 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ r> 

\ ^ 13 

N 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 



NON-LOCAL 

VARIABLES 



\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

N 


\ 

\ 

\ 

' \ 



SU3R0UT INES 

RECU IREC 


NAPE 

\ 

DESCRIPTION 



\ 





\ 





\ 





\ 





\ 





\ 





\ 





SUBROUTINE GEiMPSKC MSKl, PSK2. ARRSIZE. CPNUM) 


integer. ARRSUE. CPNUK. C 

INTEGER PS K 1 ( APRS I ZE ) . MS K2 ( ARRS I ZE ) 

IF (CPNUR .NE. 1) GC TO 10 
PSKl ( I ) = 0 

KSK I (2 ) = -1 
KSKl ( 3 ) = 1 

PS K 1 ( A ) = C 


<0 


o o n o 


r 


c 

►'SK2( I ) » 
^'SK2 « 2 ) » 
KI> K2 ( 3 ) * 
ksk: ( '. ) « 

cr tl’ :r 

PRFV. ITT ANC 
C. C=2 FCR 

IC lcntinle 
C » I 
IF (CPNUV 
C 

^ : K 1 { 1 ) = 

FSKl (? ) = 
J'SKl ( 3 ) = 
K S K M ^ ) « 
VSKl ( b ) » 
FS K 1 ( 6 ) « 
VSKl ( 7 ) = 
f'SKl (6 ) * 

KSKl ( 9 ) « 
C 

V5K2( 1 ) = 
HSK2(c) * 
>'SK2(3) = 
PS K2 ( <. ) » 
P3K2(5) = 
KSK2(6) = 
PS K2 ( 7 ) » 
PSK2(fe) = 
PSK2(9) « 
C 

20 CCNTIUUE 
RETURN 
END 

ECF.. 


- I 
C 
C 
I 


SCPEL DIFF. CFERATCRS HIFFEk C’ILY MY THc 
SC-’L-L. L»l FCR PREUITT. 


EC. 3) C 


1 
C 
- I 
C 

f 

W 

-c 

I 

0 

- 1 


origin Al* 

OF. POOR 


-I 

-c 
- 1 
c 
c 
0 
1 
c 

I 


( .SIAM 


P^GF IS 

quality 


c 

c 

r. 

c 

c 

c 

f. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


LNlVfR<nv C.F KANjAb T f L F L C f f' L M C A T I C N S AN) I N F C f< •' I I C N ''t .tNf.fb L A e> 


PPCfPflv SLITE : NCISC lllTlPS lUF. « : 


POCCRAy, NAMC-.Tl-CFNC A'nnC'>:J. Sf.CTI CM»ON>=P D M F : C <! / 2 7 / P 2 


PLIRFCSC : THIS 15 THE MAINLINE FCR Ti. TU'tSHOL') 

SEARCH RCUTINC WHICH USES A CUAHRATIC SEARCH 

TC FINC THE BEST THRESHCLC FCR FKCCL'CINC ORIGINAL PAGE IS 

AN ECGt MAI* BASEC JN THE CGCE FIGURE CF MERIT. OF POOR QUALITY 


PARAMETER OEFINITICN 


NAME 

\ TYPE 

\ 

CLASS\ RANGE 

\ 

CESCRIPTICN 


\ 

\ 

\ 

\ 



\ 

\ 

S 

\ 



S 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

S 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

’\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



NON-LOCAL VARIABLES 



\ \ \ 

\ 


s \ \ 

\ 


\ \ \ 



\ \ \ . 

S 


SUBROUTINES RECUIREC 


NAME 

\ OESCRIPTICN 


THCCUD 

\CCES THE PRCCESSING-- 
\ 

\ 

\ 

\ 

\ 

\ 

CALLED WITH VARIABLE PARAMETERS 

FRCGRAM 

ThCFNC 


INTEGER 

MAXSIZt HXWSIZ 


parameter 

( M A X S I Z = 5 1 Z ) 



parameter (MXWS12=15) 


C 


INTEGER RECS 

INTEGER TOUFF ( MAXSI2 ) , SIZE. QUTSIZ. WNCSIZ 


■ 






























lMrcni«l ICLFN(l'J», ACTLFN«ie)t llliJSfN(lB) 

C 

1 M Forw T T y I N , TT YCL T 
LCLICAL 

c 9^ boor 

ERF - . J A I 5 C . 

CAT A TlYir. . TTYCL'T /!', 16/ 

C 
C 

V^RlTf: (TTYCLT.6CI) 
bCl FCRFAT (IX, ‘ENTER ICF.AL IHACE FILFNAyE (PUST 5E AN CLC FILE)’) 
RE AC (TTY IK, 5 10) IDLFN 
510 FCRFAT ( I JAl ) 

C 

V,R I TE ( T T YfJbT , 70 1 ) ICLFM 
7Cl FCRFAT( IX, ISAl ) 

C 

ViRITE (TTYCUT, 602 ) 

602 FCRNAT (IX, ’ENTER ACTUAL IFAGE FILENAME (MIST QF AN CLC FILE)’) 
REAC (TTYIN,510) ACTlFN 
C 

V.R I TE ( TTYCUT , 70 I ) ACTLFN 
C 

WRITE ( TTYCLiT,bC3) 

6C3 FCRFAT( IX, ’F ILFNAME FGR THE EOGC ''AP (MUST BE A NEW FILE)’) 

REAC ( TTY IN, 5 *0. THPSFK 
C 

WR I TE ( TT YCUT , 70 I ) TFRSF N 
C 

WRITE (TTYCL'T.blO) 

blC FCRFAT (IX, ’ENTER Tf'E SIZE CF THE IOEAL' 1 IAGE ’) 

kEac (ttyin,<>) size 

V,RITE(TTYCUT,») SIZE 
C 

IF (SIZE .LE. HAXSIZ) GCTO 12 
WRITE ( TTYCLiT , 615 ) MAXSIZ 

615 FCRMAT (IX, ’« =>X«ERROR- - .THE MAXIMLlM SIZE = ’,15) 

GOTO ICIC 
C 

12 WRITE ( TTYCL)T,b27) 

627 FCRHAT (IX, ’ENTER THE SIZE CF THE FILTERED IMAGE’) 

REAC (TTYIN, <> CUTSIZ 
WRITE (TTYCLT, » ) CUTSIZ 
C 

IF (CUTSIZ .GT. SIZE) THEN 
WRITE ( TTYCLT, 6Al ) 

641 FORMAT ( 1 X , ’ F I L T E RED IMAGE MUST BE SMALLER THAN IDEAL IMAGE’) 
GOTO 1010 
END IF 
C 

C . 

C call SUBRCUTINE to DC THE WORK 
C 

c 

CALL THCCUO (IDLFN, ACTLFN, THRSFN, SIZE, REGS, OUTSIZ, 

C TTYCLT, TTYIN, T6UFF, £P.R) 

C 

C 


’ 1.S 

QUAUtx 


o* 


L 


630 


6 AC 

L 

C 

ICIC 

CF.. 


IM ) 6C TC lOlC 

V r 1 T E ( T I Y C L T . M 0 ) 

all r, CNF 6 6 6*) 

V.UITJ; ( rr ^ llt .6AO ) O'jtsi:, rfcs 

FCRFAT (IX, 'THE L'UTPOT IMAGE IS ' . 1 t ' WCtCS >Y '.15,' 


STCP 

6NC 


e^ALp. 

W4{'/ry 


CCPC 










(A 

or 


L fIVEPSITV CF Ka^•J..;, TfL( LCPPLMC.AT lC^S INFCRPATICK Si.ItM = S LSg 

C - - 

C PWCGRAP .>LITE ; Nri';L riLTtnS KFF. i : 

C 

C rwrCRAr NJPF :Thnt.UIJ AUTHCfl:JtrF x^TSC^ CA T F : ^ / I S /cl 

C - - - - — 

C PLt'FCSE : Tt-OCbC IS A CUACRATIC SFAPCH. MSINC A CUAC^'MC ipFPCX- 

C I>*A11CK. IT L!: TE«*'l A5S TmE ALXT THRCSHCLC TT. HAVE ITS FJGv.'Tc 

C CF KlRIT CAICLLAICO. I ASSLF’EO THAT A GRAPH OF EPH'i Vb. THRESHCLC 

C UCLIC PE APFUCX IKATFLY A 0UAC‘>ATIC. TmE SLARCh TPRMlNATtS V«HEN 
C THE tFP THAT J5, API’RCX I VA f t C •! T TMfc vjUAOkATlC IS WITHIN A S^-ALL 
r. A'-CLf.T (3FLTAI CF THfc AfTUAI. EF“ ANC THf^ SEARCH IS CvER A SHALL 
C RANGE OF THRESfiCLCS. DR WHLN A SET HaXIHLIH ITCRAMCNS CF TmF SEARCH 
C I S Rl ACHEC. 

C - 

C PARAHtTER DEFINITICN 

C HE \ type \ CLASSN RANGE \ C^SCRIPTION 


\IREAL ECCE hap 
\NCISY EDGE HAf» 

\0CST THRESHCLC ECCE HAP 
0-Sl? \NUH COLUHNS IN ICFLc 
C-512 SNUH RECS IN ACTFLE ANC THOFLE 
0-S12 \NUH COLS IN ACTFlF ANC ThDFLE 
16 \CUTPUT TC terminal FILECCOc 

15 MN'PUT FRCH TERHUAL FlLcCOCE 

\TEHPC«ARY BUFF-EP ARRAY 
\FLAG FCR FILE CPEMNC ERROR 
\ 

\ QUiuiKaE page is 

^ OB POOR QUALITY 

c 

C NON-LCJCAL VARIABLES 

C 

C \ 

C \ 

C \ 

C V 

C 

C SUBROUTINES BECUIREC 

C NAHE \ DESCRIPTION 

C 

C EFM 
r. THCRFI 
C THNEWX 
C THA2XO 
C THSPEC 
C 

C CPN 
C UCLCSr 
C FILERR 

C 

C 

subroutine THCOUOC I cl FL E t ACTFLEt TURFLCt NCCL t NCLTRCW. 

C NCUTCOL, TTYQUT, TTVIN, TB'JFt ERR) 

C 


\CALCULATES EDGE FIGURE OF HERIT FOR A GIVEN THRESHCLC 
\CREATES THRESHCLOED FILE 
\FINCS THE NEXT THRESHOLD 

\FINOS NEW INITIAL THRESHOLD AMD SEARCH STEPLENGTH 
\FCk A DIFFERENT IHAGE FINOS EFH ANC CREATES ANC EDGE 
\MAP FOR FIRST IMAGE'S BEST THRESHOLC 
\CPENS FILES 

\CLOSES FILLS OPENED wITH CPN 

\GIVES ERROR MESSAGE FCR FILE OPENING cRRCR 



c 

I CFLE 

\ 1 NT 

\Rf AC 

\ 

c 

ACTFLE 

\ I NT 

\KEAC 

\ 

c 

ThR f LE 

\ INT 

\WR 1 TE 

\ 

C 

NCCL 

\ INT 

\REaC 

\ 

C 

NCUTRCW 

\ I NT 

\hR I TE 

S 

C 

NCUTCCL 

\ INT 

\REAC 

\ 

c 

TTYCUT 

\ INT 

V<EAC 

\ 

C 

■ T Y I N 

\ INT 

\P. EAC 

\ 

C 

TBUF ( NCUTCCL ) 

\ INT 

\WR I TE 

\ 

r 

•-RP. 

\LCC 

\WR I TE 

\ 



S 

\ 

\ 



\ 

\ 

\ 

c 


\ 

\ 

\ 

c 


\ 

\ 

\ 




INTEGER NCCL. NJIJTROW, NOUTCCL. TTYCUT. TTYIM, MAXCCLNT 
integer ERRNLH. T (illF ( NCUT C CL ) . J. ICLN'.JM, ACTNUH, TH- JH 


IMfC* X, XltFI. N-^KiF-T. ALFSTAK. T. XM\|{Kt AlPHA 
If.TFCFP AlFjT, LCIAT, XO. Ai, A^' , ft b T A o , I N T »• , MAX7H, ANSWER 
C 

1 I !:i f L F ( 1 M ) . ACTFl«Hld), THiJUrcld) 

C 

LCGICAI P.ir.fl, frRU. CT PL"'. . IMCCUSC 

c 

kEAL CR‘!CAL. )firA, EFf'ARR < c : ?55 » . HAI.FSTVR, FI, F^, F2 

HFAL FFWAX. HALFST, HFp, hSTftP, FSTAR 

C 

PAKAyPlffO (VAxCCt'NT * 1C) * 

PARAi^tTER (PLLIA = .01) PACE K 

PARAMETER (KINTH » 0) QI/AMTV 

PARAME1ER(“AXTH=255) 

C 

xlgftcxvapk , alpha ) 

XR I OHT ( XH ARK , AL PHA ) 

ALFST(Fl,F2,F3,T) 

t 

HALFSK F1.)-2,F3) 
t 
C 

CALL CPM Tl'RNU^'. THRFLEt 'NEU*. 'LiNF*. ERRNLM, ERR) 

IF( ERR ) THEN 

CALL FILERR( TTYCUT, TliRFLEf ERKMJP ) 

GC TO 99 
ENC IF 

CALL CPA( ACTAUF't ACTFLEt 'OLD'. 'L'NF', ERRNLi'<f ERR) 

IF( ERR ) then 

CALL FILERKI TTYCUT, ACTFLE. ERRAUK ) 

GC Til 99 
ENC IF 

CALL CPA( ICLAUK. lOLFLE. ’CLD*. *UNF', EPRNU:< E«R) 

IF( ERR ) THEN 

CALL FILERRl TTYCUT. lOLFLE. ERRNUP ) 

GC TO 99 
ENC IF 

C. 

DC 10 J = RINTH, HAXTH 
EFPARR(J)»0.0 
10 CONTINUE 
C 

XC = 0 
A2 = 50 

X » XRIGKTI XO. A 2 ) 

C 

FI a EFR( ICLNUR. ACTNUR. NCCL. NCUTCCL. XC. TTYQIJT) 

F2 « EFR( ICLNUM. ACTNUM, NCCL. NCUTCCL. X, TTYCUT) 

EFM ARR ( XC ) a F 1 
EFMARR(X) a F2 
C 

CCLM a 2 

IF( F2 .CE. FI ) THEN 
RIGHT a .TRUE. 

ELSE 

R IGHT a .FALSE. 

XC a XRIGHK XO. A2 ) 

F3 = F2 
F2 = FI 
FI a F3 


a XRARK - ALPHA 
a XHARK ♦ ALPHA 
a I NT( (A0F2 - B-SFl - r?)itT/ 

( AOF 2 - 2 OF 3 - 2 OF 1 ) 1 
a FI - (AdS(A»F2 - 30F1 - F3))0«>2/ 
( 80( F 1 - 20F2 ♦ F3 n 



. V '•.i r tMv ir I'l.'jt i i A\r; A•>l;u^o he^f. 

t » . IF 

UM!Ul - -r Ai s‘ . . 

owg/nal pacf ,o 

u^-c, ,rc:< “‘^RQUautS 

?0 CCMIM,*' 

1C c r t. n N t F 

A 'J » <’ * <? 


C.UL PlGUr, xo. A3t MAXTh. MINTH, X. INECi'NC) 


1F( iKeCLVC ) then 

fF<tFMAiiR(X) .Fv. C.O ) Then 

Fjl = iFM( IDLNLK, ACTNUH, NCCL, '^CUTCOL, X, TTYCUT) 
EFMARI (XI = F3 

ELSL 

F3 = EFMARR( X ) 

F.NC IF 

CCLNT = CCL.'JT ♦ 1 

ELSE 

F3 * 0 
ENC I r 

IF( CEOUG ) then 

WR I TE ( TTYCUT , <,5 ) Fl,f2,F3 

<i5 FCRMAT( IX, *F1= •,F7.5,* F2= ',F7.5,* F3* *,^7.5) 

ENC IF 

IF(F3 .GT. F2I THEN 
A ? = A3 
F2 = F3 
GC TC ?C 
ENC IF 

1F((F1 .EC. F2) .AND. (F2 .EC. F3II THE''! 

WRITE! TTYOUT, 50 ) 

5C FCRMATI IX.'CRIGINAL OCUNOS OF SEARCH NEED TC 8E ChANCEC) 

GO TC 37 
END IF 

HSTAR = HALFST( F1,F2,F3 ) 

ASTAR = ALFST (F 1 ,F2 .F3 , A2 I 

|F( aSTAR .EQ. A2 ) ASTAR = ASTAR ♦ I 

CALL THNEWX( RIGHT, XO, ASTAR, MAXTH, MINTH, X, IN8CLND I 

IF( IN3CUND ) THEN 

IF( EF^'ARRC X ) .FC. 0.0 ) THEN 

FSTAR a EFM( IDLNUM, ACTNUM, NCOL, NOUTCOL, X, TTYCUT) 
EFMARR(X) = FSTAR 

ELSE 

FSTAR = EFMARR(X) 

ENC IF 

CCUNT a CCUNT ♦ 1 

ELSE 

FSTAP a 0.0 
END IF 




c 














ERff. M /.i, r. - ks'';.ri 

I r- ( ( ( L ►•• K I . o 1 . I ’ C L T /. ) . c R . i M . I". 1 . 3 ) I 
c .AM.. ( CCLM .LT. f-AAC CL\T )) T.(£\ 

r. 

c FINC MU '■.!£;» LcKT.TI' AN3 ■'CVL XO 

C 

C.UL ^H.^:AC( RIOHT, XO. ASTA3. a 2. FI. F?, FSIAS ) 


C 

Gc rr 3C 

c 

ENC I F 
C 

X = K I ;» T H 

EFMVAX = c.c 

CC 5 5 J = M fiT* , vaxTH 

IF( EFKAKR(J) ,GT. EFMMAX) THEN 
X = J 

EFRRAX = EFHARR(J) 

ENC I F 

55 CCMINUE 




V 


«» 


4 > 


V 




UR I TE ( T7YCUT ,70 ) X, EEMRAX 

70 FCRFAT( IX. •!?EST THKESHC LU , * . I 3 , ’ GIVES EDGE FIGURE CF HERIT,' 

C ,F8.5) 

C 

IF (CEEUC) then 

UR I TC ( TT YCUT , 75 ) CCUNT 
75 FCRFAT { IX , • 1 r TCO<*.13,' I T ER A T I C.N S * ) 

END : F 
C 

C CREATE TFRESHOLCED FILE 

C 

CALL THCRFM ACTNUH. THPNUR. T6UF, NCUTRCU, NiJUTCaL, 

C X, MINTH, MAXTH ) 

C 

57 CCNTINUE 
C 

KHILEI CEGUG ) 

REU INC( ACTNUV ) 

WR ITE( TTYCL'T. RO) 

90 FC9RAT ( :X ENTER SPECIFIC THRESHCLD (0-235). THRESHCLC • , / , 

C IX.'CUTSIOE RANGE TC STCP*) 

REACITTYIN, «) X 
KR 1TE( TTYCUT , >> ) X 

IF((X ,LT. 0) .09. (X .GT. 255)) GO TC 93 

FSTAR = EFP( lOLNUM, ACTNUM, NCCL, NCUTCOL, X, TTYCUT) 

WR IT6( TTYCUT ,95) X, FSTAR 
95 FOayAK IX, *FCP. THRESHOLD I 3 . ’ EFM = ',‘=7.5) 

ENC WHILE 

C. 

58 CCNTINUE 

• CALL UCLCSE ( THRNUM ) 

CALL UCLCSE(ACTNUK) 

C 

C 

C AT ABOVE BEST THRESHOLD FINC L'FH AND EOGE ■l/P FOR OIFFErENT IPAGES 
C 

ANSWER = I 

WHILE (ANSWER .FQ. 1) 


I 

r 

I 


c 

c 


C/Ll 1 H '■» t’ !: C ( X. lOlN'jy, A'JSwtfo, Nf.Cl. t I'vTurCGL. 

C M/xTM, TTY|,\, TTYGur > 


£NC ,»ML£ 

C/l. L IJCLCSt. ( I CL NU -f ) 
S'V CCNT I NCL 
RETUKN 
ESC 


ECF.. 

•> 








I ’Ik. 






I' 




i:. 


i:- 


j 


page js 

DF POOR QUALITY 




C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 




UM vr RS n Y OF K ANS/-. s ■» 

L u T 1 s r r. • . 

I M‘. L A 

frcgram 

S LI T t : 

Nd.Sr. FILTERS REF. 

1/ : 


PRCCR am 

NAME : 'F 

•' AUTHCRtJtPF 

A I S C N 

UA1E 


PUPPCCe : 


TM5 ALGCP.IT^•^' CALCULATES ThL ETCE TIGL^E CF FEPIT 
Cl AN ACIUAL EUGfc KAP AS PRCFCScO OV P«A7T, 

•CUAKTATIVL CFSIGN ANH FvALLATinN CF ENM A NC*" K C M / T H ^ E ShCL C 
L'CGE Ctlr(TGH5'. FHGCEEOnGS CF T!'£ IFF.E, MAY IS'/S. 

TFL CCGC FIGURE OF MERIT IS CrFIN^C BY 

F = I I .C/MAX( JCCCNT, ICCUNT ) ) « SU?I ( 1 . 0 / ( 1 . C ♦ ACCSO) 

THIS ririj^t aF KIRII CAN flC CSEC TC EVALUATE EITHER 
THE PERFCPKANCc CF ERCE CETECTCR5 C" THE PuLATIVE 
CUALITY CF VARiruS IMAGE EfHANCF.yFM ALGORITHMS 
hlTh RESPECT TCi ECGE CUALITY. htRE AN IDEAL ECGE MAR 
DEFINES THE EXACT LOCATION CF THL 60GE. IN THIS 
FRCGKAM I assume THAT THERE IS ONLY CNc tOGE AND IT IS 
VERTICAL. 


parameter OEFINITICN 


C 

r 

NAME 

\ TYPE 

\ CLAS 

S\ 

RANGE 

\ CSCRIPTIQN 

U 

c 

ICLNUM 

\ I 

\R 

\ 

17-AO 

\IOEAL FILE NUMBER 

c 

ACT NLM 

\ I 

\R 

\ 

17-AO 

\ACTUAL FILE NUMBER 

c 

NCCL 

\ I 

\R 

\ 

0-5 12 

\MJM COLUMNS IN ICEAL FILE 

c 

NCU' 

\ 1 

\R 

\ 

C-512 

\NUM COLUMNS IN ACTUAL FILE 

c 

TFi lC 

\ I 

\R 

\ 

C-255 

\ THRESHHCLO 

c 


\ 

\ 

\ 


\ 

c 


\ 

\ 

\ 


\ 

c 


\ 

\ 

\ 


\ 

c 


\ 

\ 

\ 


\ 

c 


\ 

\ 

\ 


\ 

c 


\ 

\ 

\ 


• \ 

c 


\ 

\ 

\ 


\ 

c 


\ 

\ 

\ 


\ 

c 


\ 

\ 

\ 


\ 

c 








C 

c 

c 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

w 

C 

c 


NON-LOCAL VARIABLES 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


FAME 


SUBROUTINES RECUIREC 
OESCP. !PT ICN 


] 


T Taun 




, LNCT 1 Of 

. actnum 


C 





1 NT EG CP 

AC T N U P . 



INI ECf R 

.G. 1*^1'. AEP 



j f. T F G £ R 

J . . I C C 1. M . 



REmL 

.NUP. A, OS 



LCG IC iL 

CEE'JG 



LATA 

FLAG/0/ 

Vir 

■ C 




c 

StT ER“CR 

ccnstam 


c 




c 

A = : .C /9 . 

0 


C 




c 

SET INITIAL CONDI T IONS 


c 

JCCLNT = n 


<# 


ICCLNT = 0 



1 ^ 

I' 

I' 

r 

I' 

I 

u 

I- 

r 

r 

I 


:i.u' 


c 

c 

c 

f 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 


Pacu 

OS foofi QUAinx 


csijy 

LEeUG 


0.0 

.TRUE. 


KAIN LOCP FCR PROCESSING UAT/S 

THIS LCCP PROCESSES ONE LINE OF THE EOCE IMAGES 
at a T1.VE 

10 CCNT I f. ;E 

R EAC V 1 CLNUH I END=<?e ) I 1 OKQVi ( I ) , 1 = l.NCOL) 

R EAC( ACTNUM . EN0 = 93 ) ( ACRCU { 1 ) f I = l.KOUTCCL) 

THIS LCOP FINOS THl- iCF.AL EDGE POINTS FOR THIS LINE 
(I ASSURQC ONLY ONE tCCE PCINT PER LINE IN THE IDEAL I PAGE ) 
ANC STORES THEIR LOCATION ANC COLNTS THEM 


FLAG) Then 
- NCUTCCD/2 
1 


20 

25 


CO 20 COL = l.NCOL 
IF I ICROiaCCL ) .EU. 
lEP = COL - (NCOL 
JCCLNT = JCOUNT ♦ 
GO TO 25 
ENC IF 
CONT INUE 
CCNT INUE 


THE FOLLOWING LOOP FINOS THE ACTUAL FCGE PCIMT (AEP) 
CCLNTS THEP. CALCULATES THE SQUARED DISTANCE BETWEEN 
IT AND THE lEP AND CALCULATES ITS CONTRIBUTION TO 
THE FIGURE OF MERIT. 


CO 30 COL = l.NOUTCOL 

IF (ACRCW(CCL) .GE. THRHCLO) 
AEP = COL 
ICCUNT = I CCUNT- ♦ 1 


THEN 


30 


CSC 
OSUH 
ENC IF 
CONTINUE 
GO TO 10 


* RFAL((IEP - AEPIOCIFP - AEP)) 

* OSUM ♦ 1. 0/1 1.0 + A#CSC) 


98 CCNTINUE 

CALCULATE THE EDGE FIGURE OF MERIT 


X^L^' » !^EAL • ICCL'M » 

1F( JCCLNT .CT. ICCLNT ) XNL'» « '<£AUJCOJNT) 
tFW » 


|F( CCeUC ) THEN 

V,R|TEITTYCUT,<,0» ThPMQLCt £FM 
FCKf-ATC IX, • TH(tEShnLC» EFK* 

END I F 

P.EH lwn( I CLNUK > 

RFW INC( ACTNUf' ) 

RETURN 

END 


• . FS . 3 ) 

dktgtnal page is 

OE POOR QL’AH'iV 




i 


•• 






•mt 






w 


3 * 


c 

r, 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

r. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 



UNIVERSITY GF KANSAS 

RiiMCTc SFNSlf.C LAc 

PRCCR AH 

SLITt : NCI Sc filters HEf 

. 1 : 

P R C G H A M 

NAMF:THCRFI AIJTHCR:JcFF 

WATSC V DA I F : 5 /? 5/83 

PURFC 5 C 

: CREATES THRCSHULOEC FILE 

original page is 



QK POOR QUALITY 



parameter OEFl.NIT ICN 


FAME 

\ type 

\ CLASSN RANGE 

\ C = 

SCRIPTIGN 

ACT NLM 

SINT 

\REAC \17-A0 

\LFN OF 

ACTUAL IMAGE 

THR MJM 

\ INT 

\RE4C \17-A0 

\lfn cf 

thresmolcec file 

TBC'F ( NCCTCCL ) 

\ INT 

\»JRIT£ \ 

\DIJFFET 

FCR PROCESSING 1 


\ 

\ \ 

SRECOkT 

AT A TIME 

NCUTRCW 

\ INT 

\HRITE \C-512 

\MJMMcR 

CF RCV<S IN ECGEMAP 

NC'TCCL 

\ I NT 

\REAC \0-5l2 

\ MIMHET 

OF CCLUMNS IN ACTUAL 


\ 

\ \ 

\ 1 mace 


THRESH 

\ INT 

\REAC \C-255 

\ THRF SMOLD 

M I N T H 

\ INT 

\REAC \0 

\ CLACK 


MAXTH 

\ INT 

\REAC \25S 

\HHITE 



\ 

\ \ 

\ 



\ 

\ \ 

\ 



\ 

\ \ 

\ 



\ 

\ \ 

\ 



NON- 

LOCAL VARIABLES 




\ 

\ \ 

\ 



\ 

\ \ 

\ 



\ 

N \ 

\ 



S 

\ \ 

\ 



SUOP.OUTIMES RECUIREC 



FAME 


OE SCR I PT I QN 


- 


\ 

\ 

\ 

\ 

\ 

\ 

\ 




SUBROUT I NE 

THCRF I ( 

ACTNUM, THRNUM, 

» TBUF, 

NCUTRGW, NCUTCCL, 

C 

THRESH, MINTH, MAXTH 

) 



INTEGER flCTNUM, THRNUHt ThRESHf NCUTRGW, NCUTCCL 
INTEGER MINTH, MAXTH, J, T3UFI NCUTCCL ) 


NCUTRCk, = 0 
10 CCNTINUF 

R EAC( ACTNUM , ENG = 30) (T3UFIJ), J = 1, NCUTCOL ) 

NCLiTRCw = NCUTHCH ♦ 1 


> 

t 


c 

V. 

v; 

v: 


OC 20 J * It UOUTCCL 
1F( TbUF(J) .GF. THRESH ) 
TBL'F(J) « MAXTH 
ELSE 

TBUFUJ » MInTH 
ENC IF 
20 CCNTINLE 

V,R 1 TE « THRNUF ) ( TBUF ( .1) t J 
GC TC iO 
C 

30 CCNTINOE 
RETURN 
ENC 

EOF . . 


THEN 


» I, NCtTCCL) 


TS 








0 


0 


0 





C - — 

C UNIVERSITY OF KANSAS REMOTE SENSING LAB 

C - - 

C FRCCRAK SLME : NCISE FILTERS REF. : 

C 

C PRCCRAf' NA^E:TFNEW>: AUTHCR:JEFF V.ATSON OArF:S/25/fc3 

C - 

C PLRPCSC : FINCS THE NEaT THRESHOLC TC HAVE ITS EFM CALCLLATEC 

C ANC CETERMINCS kvHl'THLR TtilS THRESHOLD IS WITHIN TMF ALLCwABL E 
C RANGE CF THRESHOLDS. WOGINAL PAGE IS 

C DF POOR QUALITY 

C 

c 

c parameter oefiniticn 


c 

r 

NAME 

\ TYPE 

\ CLASSV RANGE 

\ DESCRIPTION 

v« 

c 

R IGHT 

\LCG 

\REAC 

\ 

VPRESENT CIRECTICN CF SEARCH 

c 

XC 

\ 1 NT 

\REAC 

\0-255 

\INITIAL THRESHOLD 

r 

ASTEP 

\ I NT 

\REAC 

\C-25S 

\OI5TANCE FRCM XO TO 

c 


\ 

\ 

\ 

NTHRESHCLC aEING TESTED 

c 

M AX Th 

\ INT 

\REAC 

\255 

\MAXIHUM THRESHOLD 

c 

M. INTH 

\ INT 

\KEAC 

\0 

\MlNfM,HM THRESHOLD 

c 

X 

\ INT 

\WP I TE 

\-255-510\NEXT THRESHCLO TESTED 

c 

INOCUND 

\LCG 

\W 

\ 

\TRUE IF 0 < X < 255 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

r 

r 


\ 

\ 

\ 

\ 

L. 

c 

r 


NON 

-LOCAL 

VARIABLES 


L 

c 

• 

\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 

c 


\ 

\ 

\ 

\ 

c 


SUBROUTINES 

RECUIREO 


c 

c 

NAME 

\ 

OESCR I PT I ON 



c \ 

c \ 

c \ 

c \ 

c \ 

C \ 

c \ 

c 

SUBROUTINE THNEWX( RIGHT. XO. ASTEo. MAXTH, MIMH, X. IN0CUNO ) 
C 

INTEGER XO. ASTcP. MAXTH. MINTH, X 
LCCICAL RIGHT. INBCUND 
C 

1F( RIGHT ) THEN 

X = XO ♦ ASTFP 

ELSE 

X = XO - ASTEP 
END I F 


c 


ECF. 

t 


IF(( X .LT. »>INTm ) .CR. ( X .CT 
IN2CUNU = .FALSF. 

ELSE 

l^eCU^C » .TFUE . 

fcNU IF 
RETURN 
END 


. R A X T H ) » THEN 


j^GlNAL 

P® £qor quality 


n n 


c 

c 

r 



U M V L R S M Y 

OF KANSAS 

P ) ‘‘C’ 

t SENSING L At* 

r, 

r 

PPCGR Jr 

SLI T 

t : f. c ! :; *■ 

F 1 L T i P 

- ME F 

. •» : 


c 

P R C C R A M 

NAHF 

: T H A 2 X 0 

AUTHC> ; JEFF 

I.ATSCN UA It : *^/25/e3 

L 

c 

PUR PCSt 

• 

• 

Ti: FINj NEW STLP 

length ANO 

rcvr 

XO FOR CiJACRATIC SFARCF 

c 

c 






ORIGINAL PAGE IS 

c 

c 

c 

r 






D8 PCX)R QUALITY 

c 



PARAMfcTt P 

OEFINITICN 



c 

r 

NARE 


\ TYPE 

\ r.LA 

SS\ RANGE 

\ 

□ •^SCPI PT in.\ 

w 

C 

R IG»-T 


\LCG 

\MOD 

\ 

\C IRFC r 1 HN or SEARCH 

C 

XC 


\ INT 

\MOD 

SO-255 

\I.NIT1AL THRCSHOLC 

c 

ASTAR 


\ IM 

NRCAC 

\0-255 

\UISTANCF FkCH XC TC 

c 



\ 

\ 

\ 

\CUAD 

•S MEST GUESS 

c 

A2 


\1.NT 

\REAC 

\ 0-2 50 

\STEP 

LENGTH Cr SEARCH 

c 

F 1 


\ H c A L 

\MUO 

\0- 1 

\EFM 

AT XO 

c 

F2 


\RFAL 

\MDD 

\C- 1 

\ CFM 

AT A 2 FROM XC 

c 

FST AR 


\real 

\K EAC 

\C- 1 

\ EPM 

AT ASTAR PROM XC 

c 



\ 

\ 

\ 

\ 


c 



\ 

\ 

\ 

\ 


c 



\ 

\ 

S 

\ 


c 

c 



\ 

\ 

\ 

\ 


c 

r 



NCN 

-LUCAL 

VAR I AELt S 



c 



\ 

\ 

\ 

\ 


c 



\ 

\ 

\ 

\ 


c 



\ 

\ 

\ 

\ 


c 

r 



\ 

\ 

\ 

\ 


c 



SUaROLTINES RECUIRbC 



c 

r 

NAME 


\ 

OESCR I PT IGN 



c 



\ 





c 



\ 





c 



\ 





c 



\ 





c 



\ 





c 



\ 





c 

. 


\ 





c 









SLOROUTINE THA2XC( KIGHT, XO, ASTAR, A2, FI. F2. FSTAR ) 
INTEGER XC, ASTAR, A2. XRIGHT, XL£FT, XKARK, AU'HA 

lcgical right 

REAL FI. F2. FSTAR 


XRIGHTI XHARX. ALPHA ) = XKARK *■ ALPHA 

XLEFT( XHARK, ALPHA ) = XMARX - ALPHA 


IFIASTAR .G6. A2) THEN 

IF(F2 .Gr. FSTAR) THL'.N 


I 




I 

v_ 

L 

L 

h 

I" 

r 

L 

1 

K 

Ic 

K 

K 

TC 



L. 



C 

C 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

ECF .. 


PAX IS BcT,JEEN A1 ANC ASTAK 

A<» » A5TAK 
1 1 ( R I r> t ) 
h.FT 
XO 

f LSL 

.X IGHT 

XC 

ENC I F 
FI « F5TAR 

ELSE 

VaX I£ 

I F ( R ICHD 

IF(.NCT. RIGHT) XO » XLFFT(XOtA2) 

A2 * ASTAR - A2 
FI ■ F2 
F2 « FSTAR 
END IF 

ELSE 

ASTAR “cTWEEN A1 AND A2 

IF(F2 .GT. FSTAR) THEN 

PAX 3ETWFEN ASTAR AND A3 

IF(RIGHT) XO * XP IGhT( :<0. ASTAR) 
IFi.NOT. RIGHT) XO « X L E F T ( XO t A S T AR ) 
A2 » A2 - ASTAR 
FI a FSTAR 


- A2 
THEN 

» .FALS'’. 

» XRICllT(XOiASTAR) 
» .TKUr-. 

» X L L F T I X 0 » A s T A R ) 


QnQ|NA<: IS 

IB POOR QUALITY 


aF.TOEEM A2 AND A3 
XC « XR IGHTIXOi A2) 


ELSE 

KAX BETWEEN Al AND A2 


IF(RIGHT) THEN 



R IGHT 

S 

.FALSE. 

ELSE 

XO 

a 

XR IGHT( XO. A2 ) 


R IGHT 

a 

.TRUE. 

END 

XO 

IF 

a 

XLEFT ( XO ,A2 ) 

A2 = 

A2 - ASTAR 


FI a F2 
F2 a FSTAR 
END IF 

END IF * 

RETURN * 

END 


m 

J 


L 

J- 

I- 

r 

I 

I 

I 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 


UMVtRSITV OF KANS.'.' >’» 'oiL L A l' 


FOCCBAK lUTI ; NCIf-f. FlLTtOf, r F F . •» t 


FTCCRAF" ^A^'E : TI'SPCC /.MTHC^UE^F rAlSCN C A f F : *i / 2 1 / 8 3 


PLRFC5E : CRtMCS A T»-'RC SHOL C^.O FILF Fi*iCS Thf Tf*' FCB I ** A C 6 

CF Tl-C SAf^L' Size AS TMF. PREVIOUS 'ACTUAL FILE' . IT USES THE PRkVICLS 
FILE'S REST THREShOLJ AS ITS ThRfSHOLC. 


PAIUMETEP OEFIMTICN 

NAPE \ TYPE \ f.L/SS\ RANGE \ CC SCR I PT I OK 


\ 

\ 

\ 

\ 

V 

\ 

\ 

\ 

\ 

\ 

S 

\ 

\ 

\ 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


\ 

’\ 

• 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


BRWlNAr PAC.R ll 
OF POOR 01’ A MTV 


NON-LOCAL VARIABLES 


\ 

\ 

S 

\ 


\ 

\- 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

• \ 
\ 


SUBROUTINES RECUIREC 
NAPE \ DESCRIPTION 


THCRFl \CREATES THE THRESHOLOEC FILE 

CPN \CPENS FILES AND ASSIGNS LOGICAL FILE NJPHERS 

UCLCSE \CLCSES FILES OPENED BY CPN 

FILERR \REPORTS ERROR IN OPENING FILE V.ITH CPN 

EFP SFINOS THE EDGE FIGURE CF MERIT FOR A GIVEN ThRESbOLC 

\ 

\ 


• SUDROUTINS THSPECI THRESH, IDLNUP, ANSWER, NCOL, 

C NOUTCCL, TOUT, MINTh, MAXTh, TTYIN, TTYOUT ) 


INTEGER ICLNUX, ACTNUH, THRNLM. THRESH, 'ICUTC'Ji. 
INTEGER T3UF(N0JTCCL) , ANSWER, NCCL , H.lNTh, -1A/TM 
integer TTYIN, TTYOlJT, ERRNUM, KECS 

INTEGck#! ACTlF;^( i;n , THRSFMlfU 


iCtilCAL ‘•K? 











C 

c 

c 


c 

c 

c 


c 
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c 


«EAL 
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URIIFC TTYCUf, IJ ) 

10 FCk^ATl/,' IS THCRt A’IOTmER FlLTfcRtO IHAG- OF TFL SAI'f SlZc’t/t 
t .IXt'ThAT YCU WCULO like AN EP‘< A\C tOGF -1AF O'- FCS THIS • 

C • TIREShCLO?' ./ t U. • I » YCSt 0 « NO') 

R£AC( TTYIN. 0 ) ANS.vEK 
WftlUl TTYCLT. < ) ANSIiER 

IF I ANSV^ER .F.O. 0 I GC TO 6C 


•..RITE (TTYCUT, 20) 

2C FCRI'AT (IX, 'ENTER ACTUAL D-AGE FILENA'^E (MUST BE 4K CLC FILE)') 
F. EAC (TTYIN, 30 ACTLFN 
3C FCR^AT( ISU ) 

HRITECTTYOUT, h3) ACTLEI^ 

AO FCR)>ATI IX, 18A1 ) 

CALL CPN( ACTNUP, ACTLFN, 'CLO', 'LNF', ERRNU'^, ERR ) 

IF ( 5RP. ) THE.N 

CALL FILERR( TTYCUT, ACTLFIJ, cRRNUK ) 

GC TO 70 

END IF 


WRITE (TTYCUT, A‘J) 

A5 FCR f-AT ( I X, ' F I LE NAME FOP. THE EDGE HAP (HUST BE A NFw FILE)') 
REAC (TTYlfJ, 30) THRSFN 

WRl TEC TTYCUT, AO) THRSFN 

call CPN( THRNUH, THRSFN, 'NEW, *LNF', ERRNU'A, cRR ) 

IF ( EKR ) THE.N 

CALL F1L£RR( TTYOUT, THRSFN, ERRNUH ) 

SC TO 70 


C 

L c 

C UNO LFH FCR THE SPECIFIC THRESHOLD 

C 

C c 

F = EFH( ICLNU.M, ACTNUH, NCOL, NOUTCCL, THRESH, TTYCUT ) 

WRITE! TTYCUT, 50 I THRESH, F 

' 5C FCRFAT(1X, ' THRESHOLD, ', lA , ' GIVES EDGE FIGURE CF H E R I T , ' . F B . 5 ) 


C • 



c 

r 

CREATE THRESHCLDED 

FILE 

I 


call 

ThCP.FI (ACTNU.V, 

THRNUH, TbUF 


c 

C 

THRESH, MI nth, MAXTH) 

I 


CALL 

UCLTSE! ACTNUM 

) 



CALL 

UCLCSE( THR.NUM 

) 



RECS, 


NCUTCCL, 
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<^o ccntinlf 

70 LCMlfvLE 
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PARAMETER OCFINITICN 

\ type \ f.LAS.-iV RANGE \ CESCRIPTICN 
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SUBROUTINES REQUIRED 
\ DESCRIPTION 


ECUSLP \CCES THE PROCESSING — CALLED WITH VARIABLE PARAMETERS 

\ 

\ 

\ 

\ 

\ 
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■ PRCGRAP ECUFLT 

INTEGER HAXSI/.f MXViSIZt RUFCUE 


BLFCLE = KAXSI'E MAXW I NCCNS I Z5 
C 

parameter (vaxSIZ=8C0) 
parameter (MXWSIZ=15) 

PARAMETFR <BUFCUE = '^AXSI?'SMXWSIZ) 


V- 


C 






V- 










c 

C 

C/ 


c 

V/ 




c> 

c 


I K T E G f ►•■ •» I C S 

IMEGE;? 3IZC. CUTSlZf WNOSIZ 
INTEGER I N ( ^ A* S 1 Z ) f CUTj*/iXSIZ) 

C 

KEaL ^.ULUt ( liL'f '.'UO 

c 

IMcCEKfl INFr;>'(l«), CLTF^^'(l8) 

C 

1MEG6P TTYIN.TTYOUT, I NFCfOUTFC 
C 

CA7A TTYir., TTYOUT /15, 16/ 

C 
C 

KPMC (TTYCUT.oCl) 

6Cl FORMAT (IX, 'ENTER THE FILENAME FCR INPUT (KtST EE AN CLC FILE)') 
REAC (TTYIN,blU) INFNM 
KRITE (TTYCUT,700) INFNK 
7CO FLRFAT ( IX, ISAl ) 

5 10 FCR^AT ( 18AI ) 

C 

C 

WRITE ITTYCUT,6C2) 

6C2 FORMAT (lY, 'ENTER THE FILENAME FOR CUTHUT (MLST BE A NEW FILE)') 
REAC (TTYIN,510) C'JTFNH 
WRITE (TTYCUT,700) OUTFNM 
C 
C 

WRITE (TTYOUT. 610) 

610 FCRNAT (IX, 'ENTER THE SIZE OF THE INPUT IMAGE ') 

REAC (TTYIN,«) SIZE 
WRITE (TTYOUT, 705) SIZE 
705 FCRNAT ( IX, U ) 

C 

IF (SIZE .LC. MAXSIZ) GOTO 12 
WRITE (TTYCUT.615) MAXSIZ. 

615 FORMAT (IX, TERROR- - THE KAXIMUH SIZE = ’,15) 

GOTO lOin 
C 
C 

12 WRITE (TTYOUT, 027) 

627 FORMAT (IX, 'ENTER THE SIZE OF THE FILTER WINCQW) 

REAC(TTYIN, »)WNCSI2 
WRITE (TTYOUT, 705) WNDSIZ 
C 

c 

20 IF (WNOSIZ .LE. HXWSIZ) GOTO 13 

WRITE (TTYCUTf6Al) MXWSIZ 

6A1 FORMAT (IX, •« «« ERROR - - M«AX1MUM WINDOW SIZE = ',15) 
GOTO 1010 
C 

13 CUTSIZ = SIZE - WNDSIZ ♦ I 
0 

C 

C 

C CALL SUBROUTINE TO 00 THE WORK 
C 

CALL ECUSUB (INFNM, OUTFNH, GUEUE, 

C IN, CUT, SIZE, RECS, WNDSIZ, CUTSIZ, TTYOUT) 

C 
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C 

C 

V.aiTf ( TTYCL'T.n 20) 

6 30 F C A ► A T ( I << t ' ALL D C N E i:< 0 • ) 

WRITS ( 1 T VCUT I (,A0) RFC5 

6<i0 FCRMT (IXt’THC oUTPuT I^.i■•^G‘= 13 'flS,* wC7CS U '.IS,* RtCClCS') 

0 
C 

GCTC ICIC 


ICIC STCP 
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ECF.. 

7 


DRlGiNAL PAGE IS 

QB £Oor quality 


non 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

w 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


UNIVERSITY GF KANSAS 

TEL'=UCNML\I CAT ITNS AN) INFC-RHAIICN SCItNC‘=S LA0 

FRCGR Af' SU ! U : (.Cl S£ 

filters P. E F . : 

PRCCRAR NAME : FCU'SUI'. 

AIITHCiT:J. SCCTT GARDNER C A T E : 1 / -i / =» 2 
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PARAMETER 

OEF IMT ICN 



NAME 

\ TYPE 

\ CUASSS RANGE 

\ gescriptitn 


I NFNH 

\CHA1 J 

\R 

\ 

\INPIJT FILENAME 


OUTFNM 

\CH 


\l» 

\ 

\CUTPUT FlLENAMt 


CUEUE 

\R 


\R 

\ 

\ IMAGE CATA CUEUF 


IN 

\ 1 


\W 

\ 

\ INPUT eUFFER 


OUT 

\1 


\W 

\ 

\CUTPUT BUFFER 


SIZE 

\ I 


\R 

\ 

\IMAGE SIZE 


RECS 

\ 1 


\W 

\ 

\ NUMBER OF REC.TRCS 

IN imace 

WNCSI 1 

\ I 


\R 

\ 

\SIZE OF local area 

W 1 NCCv. 

OUTSIZ 

\ I 


\R 

\ 

\SIZE OF CUTPUT IMAGE 

TTYCUT 

\ i 


\r< 

\ 

\CUTPUT TC TERMINAL 

F ILECCOE 
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\ 

\ 

\ 
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\ 

\ 

\ 
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\ 
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\ 
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\ . 
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NON 

-UOCAL 

VAR I AdLES 
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SUBROUTINES PECUIREC 



NAME 

\ 


OESCR I PT I ON 



CPN 

\CPEN FIUE AND ASSIGN FIUECOOE 


UCLCSE 

\CLQSE 

FILES 

OPENED WITH CPN 


F ILERR 

\RE 

PORT 

TYPE 

OF FILE ERROR 



.SUBROUTINE ECUSUU (INFNM, CUTFNM. OtCU'^ t 
C INf OUT» SI7£t RECS. WNOSIZ. CUTSIZt TTYCUT) 


INTEGER SIZE. OUTSIZ. WNOSIZ 
C 

INTEGER ERRNUM 

INTEGER IMSIZE). OUT(CUTSIZ) 

INTEGER RcC. WRO. WORD. RECS. CREC. Tl^PC^C. START 



c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


30 

AC 

c 

c 

c 

c 


c 

50 

c 

c 

c 

c 


INTEGER SREC. WR'TRCt OVCRC 

REAL C'JEUt ( V.Nr;S I ; f S I Zr I 
REAL WNCRISt <NC'JGT, TCTALM, XEST^'T 

INTEGER^l INENKimi ClUlrNK(!‘n 

LCGICAL cRR 

INTEGER TTYGL'T , IN'FC.OUTFC 

CPEN FILES ANC CHECK FOR ERRORS 

CALL CPN ( INFCt INFNH, 'QLC’t 'UNF*, ERRNL'^'t ERR) 

IF (6RP. ) COTC 1009 

CALL CPN ( CUTFCt OLiTFNKt ‘NEW’t 'L'NIF', ERRNLMt ERR) 
IF (ERR) GCTO 2009 

WNDPTS = WNCS IZ « WNDS I Z 



CALCULATE THE WEIGHTING FACTOR 
WNCWGT = l.C / WNDPTS 
INITIAL!: the circular QUEUE 

CC AC RcC- 1 ,W.\DS I Z 

READ ( 1 NFC . ERR = A009 ) ( I N ( WRC ) » Wft C= 1 t S 1 Z E ) 
CC 3C WCRC=l,SI2f 
CUEUE(REC.WCRC ) = IN(WORC) 

CCNT 2 MJE 
CCN T I NUE 


BEGIN PROCESSING 

REGS = 0 
QREC = 1 

TRPQRC = QREC 

CC 90 START=l .nuTSIZ 

TCTALM = 0.0 


CC 2C kREC = I .s.NOS I Z 
kWCRD = I 

CC IC CkCRO=STA RT , S T AR T ♦ WND S 1 Z - I 
TCTALM = TCTALM ♦ QUE LE ( THPQRC . CkCRD ) 0 WNQWGT 
WWCRD = WWORC ♦ I 


L c ^ T I ix ■„ C 



T^'PC^C = “CC WKCSI!) ♦ 1 

CCNT INUE fS 

xestmt . rniALM 

PLT THE FILTEOEC VALUE IN THE OUTPUT liUFFE.< 

CLT(START) « INT (XESTMT) 

continue 

WRITE THIS RECQRO AND UPCAT? THE RCCCRC COUNTER 



WRITE (CUTFC) ( GU T ( WR 0 I . WP. C * 1 . CUT S I Z ) 
DON'T BOTHER TO CHC-CK FOR A WRITE ERROR 

RECS » RECS ♦ I 


REAC A NEV, RFCCRD INTO THE (JUEU6 AND UPDATE THE FRONT-ENC POINTER 



REAC ( I NFC I £KC = Z0C,lRR = 40C9 ) ( I N ( WR C ) t WRQ » I , S I Zc ) 

CO ICO WC»C=l,SI7t 

QUEUE ( OREO .WORD) = IN(wnRC) 

CONTI NU6 

CREC = HOD (CRcCt WNOSIZ) ♦ 1 
GOTO 50 


WE'VE REACHED AN ENO-CF-FILE SC THE OUTPUT NUMBER OF 
RECCRCS WILL BE WNDSIZ-l LESS THAN THE NUMBER INPUT 
WRAP IT ALL UP AND QUIT 


CALL UCLCSE ( INFO ) 
CALL UCLQSE (CUTFC) 


GOTO ICIC 


u 

ioo<; 

I 

[ 


CALL FILERR (1TY0UT, INFNM, ERHNUM) 
GCTC IClC 


zee*; 

c 

c 

*1009 

660 

c 

1 OIC 
ECF .. 
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CALL riL£XK (TTYClJTf CUTFMP, CKRMJf') 
GCTC ICIO 


WPITc (TTvCUT.663) 

FCPf-AT (IX, 0 « ER^'.CR IN SfAClNC INPUT IPAG; f » <*» 


RETLRN 

END 
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PLRFCSC : THIS FUCCPaH CWFATES AN CUTPCiT IHAGE CF MECUNS CF 

VARIA«1LS SIJEC WINUCHS CN THE INPUT 1KAC;E 


hAHE 


PARAMETER OEFINITICN 
\ type \ CLASSV RANGE \ C5SCRIPTICN 
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NCN-LOCAL VARIAELES 
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\ \ 

\ \ 

\ \ 

\ \ 
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NAME 

\ 

SUBROUTINES RECUIREC 
OESCR IPTICN 



MEDIAN 

\FROCESS6S THfc IMAGE 
\ 

\ 

\ 

\ 

\ 

\ 




PROGRAM MECFLT 

INTEGER MAXSIZ. MXKSIZt SLFGUEt BLIFHNC 


eUFCUE = KAXSIZE * MAXH I MCCWS I ZE 
SUFHNC = MAXNI NOCWSIZE <= M A X V. 1 N OCX S I Z E 


parameter ( KAXS IZ=512 ) 



. tF'vLH «M/'XS I I Z » 

PiPA-*CU3 ldLFHKJ»MX.*';iZ«-^XWSIZ) 

IMEC'"'; T'UjF r^AXi I z ) 

IMECE-^ SIZE. CuTSU. Wh.tlSIZ. AftOSIZ 
IMECEk OlEUF.« ' iUFCUi. ) . k I NCCW ( (^UF kNC ) 



IMEGEROl 1AFN>M1R). CUTFN'f<(lS) 

INTEGER TTYIN. TTYCUT 

DATA TTYIU. TTYCUT /15. 16/ 


kRITE (TTYCUT. oCl) 

FGFNAT (IX. 'ENTER THE FILENAME FOR INPUT (NLST BE AN CL3 FILE)') 
READ (TTYIN. 5101 INFWK 
FCRRAT (13AI) 


WRITE (TTYCUT, oCiO 

FCRRAT (IX, 'ENTER THE FILENAME FOR OUTPUT (MUST BE A NEk FILE)') 
REAC (TTYIN, 510) QUTFNW 


WRITE (TTY CUT, 610) 

FORMAT (IX, 'ENTER THE SIZE CF THE INPUT IMAGE ') 

REAC (TTYIN, «) SIZE 

IF (SIZE .GT. MAXSIZ) THEN 
WRITE (TTYCUT, 615) MAXSIZ 

FORMAT (1X,'=» 06ERR0R- - THE MAXIMUM SIZE » '.15) 
GCTC ICIC 
END IF 


WRITE ( TTYCUT ,f>27 ) 

FORMAT (IX, 'ENTER THE SIZE CF THE FILTER WINCCW) 

REAC (TTYIN, » ) WNCSIZ 

IF (WNCSIZ .GT. MXWSIZ) THEN 
WRITE (TTYCUT, 62n) MXWSIZ 

FORMAT (IX, '« O* ERROR MAXIMUM WINCCW SIZE * ',15) 

GCTC IClC 
END IF 

CLTSIZ * Size - WNCSIZ + I 
ARRSIZ s WNCSIZ « WNOSI Z 
REGS * C 


CALL SUBRCUTINE TO 00 THE WORK 


CALL MECIAM INF.;M, UUTFNK, ARRSIZ, 
C WNCSIZ, CUTSIZ. TTYCUT, CUEUE, 


SIZF, RPCS. 
WiNOCw, TlUF) 


630 


W 9 1 T ( I I 'i t I t <• i<J ) 

fCH^ATtlx.'^OT;! ^LL CCNE 6 6 i»M 
h«ITC CTTYCLT OUTSIZ, HECS 

6«C FORMAT llY.,*Tht OUTPUT IH6GC IS HC^CS TY 'fl5i' 

C 

ICIC STOP 
ENC 


PAf»r re 

Of fOOH 


RCCC-JCS • ) 


EOF 
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C 

C PncCkAP SUTc : NCISE rlLTf;0'j TEF. « : 

C - - 

C FftCCRAf' NAkE:vtCIA\ AUTHC-CJtFF «AISC\ C A I t 2 


c 

c 

c 

c 

c 

c 

r 

PCRFCSF ; 

callec ey 

CCES THE PRCCLSSINC FCR TmE ''F.OlAf. FILTER AFTi.R TEINC 
PRCGR/.y <ECrLT. 

c 


PARAMcTbP OEF INIT ICN 

c 

r 

NARc 

\ TYPE 

\ CLASSY range 

N OESCRIPTI^'N 

c 

INFILE 

\CH*ie 

NR \ 

\ INPUT filename 

c 

OUTFILE 


NR N 

NCUTPUT FILENA^^E 

c 

ARKSl ZE 

\I 

NR N 

NNUMtiER OF ELEMEMS IN WINOCW 

c 

NCCLS 

\ 1 

NR N 

NCCLUMVS IN INPUT IMAGE 

c 

ncutrcm 

\ 1 

\W \ 

NROWS IN CUTPUT IMAGE 

c 

WIN'CS IZE 

\I 

NR \ 

NWINOCWS I ZE 

c 

NCUTCCL 

\ 1 

\R N 

NCCLUMNS IN CUTPUT IMAGE 

c 

TTYCUT 

\ I 

NR N 

NCUTPUT TC TERMINAL FUECCOE 

c 


\ 

N N 

N 

c 


\ 

N \ 

N 

c 


\ 

N \ 

N 

c 


\ 

N \ 

N 

c 

r 


\ 

N N 

\ 
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c 

r 


NON 

-LOCAL variables 


1# 

c 


\ 

N N 

N 

c 


\ 

\ N 

N 

c 


\ 

\ \ 

N 

c 

r 


\ 

\ - \ 

N 

c 


SUaROUTINES RECUIREC 


c 

f 

NAPE 

\ 

OESCR 1 pticn 


c 

heofnc 

\F INDS 

MEDIAN OF LCCAL 

V.INOCW 

c 

CPN 

\CPL NS 

file and ASSIGNS 

LCGICAL UNIT 

c 

UCLCSE 

\CLOSES 

FILES OPENED WITH CPN 

c 


\ 



c 


\ 



c 


\ 



c 

c 


\ 




C 

c 


c 

c 


• SUenCUTINF ‘^ECIA.N ( INFILCf CUTFILE. AKPSI2E, NCCLS. NQUT«CW, 
C WINCS12Et NOUTCCLt TTYCUT, Q, A.'t‘»VAL. TOUF ) 

LCGICAL CRR 


INTEGER ARRSIZEt NCOLSt NCUTCCL . WINCSIZE 
INTEGER l« Jf K. ERRNUH, START, ARRPCS, CUTCCL, NRCW'S 
INTEGER NCL'TRCW, TTYOUT , ARRVAL ( ARRS I ZE 1 , TRLF(NCCLS) 
INTEGER 0 ( V, I NOS I ZE fNCCJl S ) , INFNUK, QUTNU'', l'€U* “EOPCS 
INTEGER*! lNFlLc(l8», CIJTF1LE(13) 


oon o o o o nnnono ooo 


C 

KCUMC.' « 0 

c 

CMl O'N (!\rK','«. INFIIE, •r3L'*/'t E 1 < » 

GAIL CP^- nuu;v. CUTFlie, •KFW'i *UNF'. FR'^Mjyt 

IMlIALl^t i.’JfcUt; 


Cr 10 ! •« I . V. 1 ACS 1 Zc 

Bt Atll I NFNGVt t'iC » «;3) ( T Tl. F ( K ) . 

CC ?C J « IfN'CCLS 
C ( i . J ) * T JUF ( J I 
2C CCNTINL? 

3C CCNT 1 NLF 


yAlN PBCCESSINC 


I . NCCl S > 




PEDPC5 « ARRSIZE/2 
AO CCNTINUE 

CC 70 START » I, NOUTCCL 

CC 6 C I » I . M I NOS I ze 

CC fO J « START, start ♦ WINOSIZF - I 
ARHPCS * W I NOS I ZC «( I - I » ♦ J-START + l 
ARBVAL ( ARP POS ) » C ( 1 , J ) 

50 CCNTINUE 

60 CCNT INUE 

call MECFNC( ARRVAL ,y6D. ARRS 1ZE» KECPCS) 
COTCOL = START 
T»UF( CLITCCL ) = (-'EC 
70 CCNTINUE 


NCUTRCW = NCUTRCW f 1 

ViRl TE ICUTNUH ) (TBWFIK), K = 1, NOUTCCL) 
UPDATE CLEUE 


CC 80 1*1, 'rtINCSIZE 

I F ( I .LT . V. I NOS I ZE ) THEN 
CC 65 J = I, NCCLS 
C(I,J) * Q(H-l,J) 

85 CCNTINUE 
ELSE 

REACi INFAUH.ENC = 98) (TBUF(K), K = 1, NCOLS) 
DC 87 J = 1, NCOLS 
C ( I f J ) = T3UF ( J ) 

87 CCNIINUE 

ENC IE 
80 CCNTINUE 
GC TC 40 

HE HAVE REACHED THE END CF THE INPUT FILE 

98 CCNTINUE 

call UCLCSE ( INF NU'3 ) 

CALL UCLCSE ( CUT ) 

RETURN' 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


L 1 V F R S 1 T Y C f K 4 N S A •> 

I FLEC CPPLM CAT H'nS 

A NO I^FC^^'AI1L•^ ECIINC'S LAd 

r'»cGRAH SL ! Tt : NCI 

r ILTtPj PFF. 

•t : 

F P C C K A f' f. A P r : ' i C r N C 

AUThCR : JEFF 

IwU S C '< fj A T E : ? / S / 5 3 

PURFCSL : Flf.C^ T 

HERE TttE ARRAY 13 Tl-f 

HL .irniAN value cf 
LCCAL area WINCCV, 

TPc ARRAY PASSct: T C IT. 

DFK3INAE PAGE IS 
OF POOR QUALITY 




P AP A'^cTE R 

'JEFIMITICN 


NAPE 

\ 

Type \ r.LA3 

S\ RANGE 

\ cescription 

ARR 

\ I 

\R 

\ 

NAR^AY OF CLCPENTS IN WINCCK 

PEC 

\ 1 


\0-255 

\PEDIAN CF ARRAY 

S.ZE 

\ I 

\R 

\ l-AOC 

\NUPRE'’. OF ELEHENTS IN ARRAY 

MCCPCS 

MNT \RCAC 

\ 1-20C 

\PCSITICN OF MECIAN VALLE IN 


\ 

\ 

\ 

\CRCEREC array of PIXEL VALLES 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 



NON-LOCAL 

VAR I ABLES 



\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ ’ 



SUBROUTINES RECUIREC 


NAPE 

\ 

Dt SCR I PT I CN 



\ 




• 

\ 





\ 





\ 





\ 





\ 





\ 





SL'URCUTINC KEDFnC(/ 1RR. MFC. SlZEt MtCROS) 
IMECEP SIZE, 1, MEC, MEDPCS 

INTEGER ARRISiZE), P I • ARR ( C: 255 ). C.C;.NT 

ECR 1 = 0, 255 

PIXAPRl 1 ) = 0 

END FCR 


1 


C 


/ 


I 

I 

I 


c 


c 

c 


ECF . . 

7 


f C I - ^ • '> I ^ f 

( 1 ) ) * P1XA“R( I 

EN'J KCk 

I « -1 
C C L M = 0 
CC 

1 = 1*1 

CoLNT = CCLM ♦ PIXAHP.l 1 ) 

LMIL ( CCLNT .GT. VcDPCS ) j 

PEC = 1 I 

RFTlJP^ 

ESC 


I 



) * I 


DRTOINAE PAGE IS 
DE POOR QUALITY 


I 


o o o n n n 


DKIUINAI! Pii(CE IS 
OF POOR QUALITY 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c' 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


L.\IVERSITV Cr K j T ‘ L E L C ► N L M C A T I C \ s INFC^KATICN SCJENC-S LA2 


PPCGRAM i,LlU : MLTl.PS RtF. « : 


PRCGRAf' NA^E:ltEflT AlMHC'<:J. SCCTT GARTNFu L A I E : C 2 / 1 2 / A 3 


PfRPCSC : THIS INCCRPCRATES LCfS R LI L T 1 P L I C A T I w C 

NCISC yCCCL INTC HIS rCGC- PHESCRVINC ADAPTIVE FILTCP. 


RAPE 


PAR A PETER OEFINITICN 

\ TYPE \ CLASSV RANGE \ CESCRIFTICN 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

S 

\ 

\ 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


NON-LOCAL VARIABLES 


\ \ . \ 
\ \ \ 

\ \ - \ 
\ \ \ 


\ 

\ 

\ 

\ 


SUdROUTINES RECUIREC 
RAPE \ DESCRIPTICN 


LEESL'6 \CALLED WITH VARIABLE PARAMETERS 

S 

\ 

\ 

\ 

\ 

\ 

\ 


PftCGRAP LEEFLT 

INTEGER PAXSI2, MXwSIZ. OUFCUE. 5LFViN'C» TUFP'wN 


BLFCUt = PAXSI2E * ’■1 A XVJ I NCOW S I 2 E 

BL'KPkN = 8 « PAXW HmOCWS I 2 E P.AXV. I NCCVIS I 2E 

BLFV.NC = PAXW I NCOWS I 2E •> P AX V, I NPCWS 1 2 E 

parameter (PAX S 12 = 512) 


IF (w MSI.’ .i . .wSIJ) CJ1C li 
V.RITC ( TTYCUT , o.' n ■•XwSIZ 

6?8 FCRKAl (iXt'^ £KRC» - - •«A'(r<C’* tlNOCh SU^ = ’,15) 

CCTC ICIC 


CUT 3 1^ = SI Zb - v-\CS I / ♦ 1 


CALL SLMr.CLT I NE TC CC THE WHUK 


OWGJNAL PAGE IS 

POOR QUALITY 


CALL LLESLiiv (INFj^, f)UTFf,M, CUEUE . WINCC.^, ^'SK^.^C^ IN, CLT, 
C SIZE, REGS, WNJS17, OUTSIZ, THShLr., NUhLKS, TTYTUT) 


UPITc <TIYCUT,63C> 

63C FCRFAT ( I < , • « «« ALL CCNE <■« ■> • ) 

WRITE ITTYOUT,OAO) OUTSIZ, RECS 

6<iC FORMAT (1X,'TFE OUTPUT IKACE IS MIS,' WCRCS «Y RECCRCS') 


CCTC lOlC 


IClC STOP 
END 

ECF .. 




DRIGTNAE page 0 
OF POOR QU^T.TTY 




L’MvfRSlTY r.F KiJr.'iAS T( LI?CCMMLM r,.>T AM' I N F i; K A I I : N 


FRCG.lil^ 5L1TF. : NC 13L FILTfPS 


KEF. " : 


C PROGRAM r.AKf :LE ESlu; 


ADTML-MJ. SCCTT GAWnN-iP 0 A T ^ : C 2 / 1 2 / " 3 


PtRPCSt : This SUHROUTIM COGS THE PRCCESSINC 

FCR LEE'S ECGC FILTER AFTER JEING CALLET 3 
BY FRCGRAK LCEFLT VI IT H VARIABLE ARRAY C I K^^N'S I ONS . 


FAKE- 


PA->A‘*ETEH OEFINITICN 
TYPE \ CLASS\ RANGE \ 


D 6 SCR I PT IGN 


C INFKP 
C OUTFNV 
C CUELE 
C ViINCCU 
C MSKVNC 
C. IN 
C’ CUT 
C SIZE 
C REGS 
C ViNCSlZ 
C CUT SI 7 
C ThShLC 
C ^U^'LKS 
C TTYCUT 


\CH« 10 
\GH« 13 
\R 


\R 

\R 

SR 

\R/V( 

SVI 

\W 

\W 

\R 

\W 

\R 

\R 

SR 

SR 

SR 


S INPUT FILENAME 
SCUTPUT FILENAME 
SlhAGE CATA CU^OE 
SLCCAL AREA KIMCCV* 

SEDGE TEMPLATES 
S INPUT buffer 
SCUTPUT RUFFET 
SIMAGC SIZE 

SMJMBER Cr REGHRCS IN IMAGE 

SSIZF OF LOCAL APEA V, 1 NCCV* 

SSIZE CF CUTPUT IMAGE 

SEDGE ThRcShCLD 

SNUM'JtR CF ICUKS 

SCUTPUT TC TER'IINAL FILECCOE 


NON-LOCAL VARIABLES 


NAME 


SUOROUTINES BECUIRED 
S DESCRIPTION 


C GNMSKS 
C LCSTAT 
C SUBMSK 
C FNCEDG 
:C EGSTAT 
C ESTMTE 
C CPN 
C UCLCSE 


SGENERATES MASKS 

SCCMPUTES LOCAL STATISTICS CF V^INCQW 
SGETS A 3 X 3 LOCAL MEAN FRCM WIMCCV^ 

SGETS EDGE ORIENTATION AND SELECTS A MASK 
SCCMPUTES LOCAL STATS CLTSICE THE EDGE 
SFINDS AN ESTIMATE FOR THE SIGNAL 
NCPcNS FILE AND ASSIGNS LCGIGAL UNIT 
SCLOSES FILES OPENED WITH CPN 


SUBROUTINE LELSU3 (INFNM, OUTFNMt QUEUEt V<INOUVtt MSKV»NC. 

C INt OUT, SIZE, RECS, WNDSIZ, OUTSIZ, THSHLO, NUMLKS, TTYCUT) 


INTEGER SIZE, OUTSlZ, WNOSIZ 


T 


INTEGER MSK 3 X 3 C 3 , 3 ,A) , SBARtA( 3 , 3 ), MSKl,\C ( . 1 , WNOS I Z , V>NC S 1 Z ) 
INTEGER IN(SIZE). CUTICUTSIZ) 

INTEGER RFC, WRD, V,ORD, RECS, CRFC, TMPQRC, START, lEKR 


nno on n n n o no nnnon non noon nnnoo 


C 


IMrJcK V*; I, J. K, 


4 


c 

c 

c 


REAL 
REAL 
R*: AL 
REAL 


CUEUt ( V»\PS I ? I S 1 7E ) , W I NCCW ( liNCS 
SUinMS. I, /MEAN, VAR/, 
IHKSMX, 7MHS.11, TTLTMN, TTLTSO, 
IGTMKS, NMliPTS 


! , hNC S I / ) 
WUML^‘;. lhSHl'7. 
n L'tMN 


/O^iTy T 


iNTtjEK<-i I.^F^y^a), CLirNMiitj) 
LCGICAL FR‘! 


DraCINTAE PACTE IS 
OF POOR QUALITY 


integer TTYCUT , IN'EC,CLj7rC 


IMTlflLl/E TmE 3*3 GRADIENT RIN'TCWS 


DATA ( ( (VjK3X3( 1 , J,.<) , 1=1,3), J=l,3), K=l,A) 

C /l,l,l,C,C,C,“l,“l,“l, 0,“l,“l,l»C,“l,l,l,C, 
C l,C,”l,l,0,~l,l,0,”l, lfl,0,l,C,“'l,0,~l,~l/ 


initialize ZK.EANOZKEAN/VARZ MIN AND MAX ANT TOTALS 


THRSMI 
T H R S K X 
TTLTMN 
TTLTSD 
TTLHMN 
NMHPT S 


lOCOO'J 

- I CCCOC 

C.C 

O.C 

C.C 

O.C 


CETER>'IM: THE ND.-IBE'! CF PCINTS IN IMAGE (ASSUMED SO'JARE) 
NUKPTS = OUTSIZ <> OUTSI Z 


OPEN FILES ANC CHECK FOR ERRORS 

CALL CPN ( INFC, INFNH, ‘nLC*, ’UNF*, ERRNLH, ERR) 
IF (ERR) GUTC 1009 


CALL CPN ( CUTFC, CUTFNM, 'NEW, 'CNF', ERRNUM, ERR) 
IF (ERR) GOTO 20C9 


FIRST, GENERATE THE EDGE TEMPLATES 
CALL GNMSKS (MSKWND, KNDSIZ) 


WNDPTS = WNCSIZ 0 WNOSIZ 
SUePTS = (WNDSIZ/2*1) « WNCSIZ 


initialize the CIRCULAR CUE'JE 



CC AO REC=l, WNCSIZ 

READ ( INFC ,ERR = 3009) ( I .N ( WRO ) , WRC= 1 , S I ZE ) 



/ 


cc 3c :zt 

CffUE ( pf C ' ) « INC'CKL) 

3 0 CHM I N'JE 

<»C CCMINlE 

C 
C 

C BEGIN PACCES:iNG 

c 

RCCS * C 

CBtC « I 

C 

5C TPPCRC * CR5C 

C 

DC <;C STAPTsl.CuTSi: 

C 

C GET THE LCCAL STATISTICS FCR Tl-F. AREA CEFINFC BY THE VilNCCV^ 

C ANG FILL THE W 1 NOOW ARRAY TO QE USLC LATER IF AN EGGE IS F CLINC 

C 

CALL LCSTAT (CUEUEt V^MCUW, TPPCRCt START, WNOSW. WNCPTS, 

C SIZE, ZKEAn, VARZ. Z) 

C 

C 

C FIND THE EDGE THRESHQLO VALUE 
C 

EGTHRS * ZMEAN«ZMFAN/VARZ 

c‘ 

C 

C UPCATE RLNMNC SU'<S FCR CETERPINING THRESHCLC 
C 

TTLTMN » TTLTKN ♦ E G TH R S/ NU P P T S 
TTLTSC * TTLTSO ♦ EG TH R S^E GTH R S/ NU PP T S 
THRSPI » PIN (THRSHI, EGTHRSI 
THRSMX = PAX (THRSMX, EGTHRS) 

C 

C CETERPINE IF A.\ EDGE EXISTS BASED ON THE LOCAL STATISTICS. 

C THE NUP3ER CF LOOKS AND THE USER SPECIFIED THRESHOLD. 

C 

IF (EGTHRS .LE. NUMLKS- THS HLD ) GOTO AA 
TTLHPN=TTLHPN ♦ EGTHRS 
NPHPTS * NPHPTS ♦ l.O 
GOTO 60 
C 

c 

C WE HAVE AN EDGE, SC PRCCEED WITH THE EDGE FILTERING 
C 

C GET THE 3X3 SUBAREA LOCAL MEAN 
C 

4A CALL SUePSK (WINDOW, SSAREA, WNOSIZ) 

c 

c 

C FIND THE EDCr OR I EN TAT I ON AMD DETERPINE WHICH EDGE 
C TEMPLATE TC USE IN CALCULATING OUR NEW LOCAL MEAN 
C 

CALL FNCFDG (SHAKEA, PSK3X3, MSKNUP) 

C 

C 

C CCPPUTE A NEW LOCAL MEAN AM) VARIANCE USING THE APPROPRIATE 

c template 

C 

CALL EGSTAT (WINDOW, HSKWND, WNOSIZ, MSKNUP, 




r 


vA«z 


c 

c 

c 

c 

6C 

c 

c 

c 

c 

c 

c 

90 

c 

c 

c 

c 

c 

c 

c 

c 

c. 

c 

c 

c 


ICO 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

2CC 

c 

635 


C 

637 

638 
C 


/ ) 


FINO AN cSTH^AfP ru,? T!|c SIGNAL. XtSTyT 
CALL tST'Ml; VAKZ, 2. NUKLKS. 

IF (<EST^'7 .CT. 2 •^5.0) X»:STJ^T«255.C 

PLT THE rlLftPLC VMLE IN THE CUTPOT GUFF 
C L T ( 5 7 A P T ) « I •; T ( X r S T M T ) 


V LSI U ) 


original page is 

OE POOR QUALITY 


CCNT 1 NUE 


WRITE THIS RECCPI) AND UPDATE THE RECCRC CDUNTErt 
WRITE (CUTFC) (I'UT (WHO) ,URC = l.CUTSI2 ) 

DON'T BOTHER TO LHcCK FOK A WRITE ERROR 
RECS « RECS ♦ 1 

REAC A NEW RECORD INTO ThE DUEUE AND UPDATE THE FRCNT-ENC POINTER 

READ ( INFO . EN‘C = 20C.cRR = 30C9 ) ( I N ( WRC ) . WRO= I . S I ZE ) 

DC 100 WOR 0=1. SIZE 

OLEUF I CRFC . WORD » = UJlunaO) 

CCNT I NUE 

CREC = yCD (QREC. WNOSIZ) 4 1 
GOTO 50 


WE'VE REACHED AN ENO-Cr-FILE SC THE OUTPUT NUMBER OF 
RECCRCS WILL BE WNDSlZ.-l LESS THAN THE NUMBER INPUT 
WRAP IT ALL UP AND QUIT 


CALL UCLCSE ( INFC ) 

CALL UCLCSE (CUTFC) 

WRITE (TTYCUT.635) THR S M 1 . THR S M X , T T L T MN . T T L T SO-T TL T MN«T T L TMN 
FCRFAT (lX.'ECGE THRESHCLO HIN = '.F9.3. 
t /lX.'CCGE THRESHOLD MAX = '.F9.3, 

C /IX. 'EDGE THRESHOLD MEAN = '.FO.S. 

C /lX.'FCGE THRESHOLD STANDARD OEVIATICN = '.r-9.31 

WRITE ( T TYCUT . o37 ) T T L HMK/ N MhP T S 

FCRMAT ( IX . ' HCMOGENECUS AREA MEAN = '.F0.3) 

WRITE (TTYCUT.63B) NHHP TS« I CO . 0 / NUMP T S 

FCRHAT <iy.F9.3.' PERCENT CF ThE IMaGE WAS HCMOGENECL S ' ) 


GCTC ICIC 


c 


c 

c 

1 cc^ 

Cil L 

FILERT (TTYCUT, ISi'NK, FRRNUP) 

OWcinal Her 

c 

2 cos 

GCTC 

IC 1 C 

CAL 1 

riLcRP iTTVU’UTt CuTF^:H, ERRMJR) 


#• 

GC I C 

ICIC 


L 

3COS 

WRITE 

( T T YCUT . <ibO ) 


b6C 

FCKf-AT 0 0 tPKCP IN REACING INPUT 

I RAGE >> * <!' ) 

c 

ICIC 

RETURN 



ENC 



6CF.. 

•> 





c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


L N I V r 

ITY OF <AN,AS 

TFlECCMMUMCATICNS AND INFC-MAlll.N 

SC I 1 NU S L A i> 

PRC CP. A- 

SUTL : NCI SI 

FILTLRS RFf-. ‘t 


program 

N A M E : r S T y T .• 

AUTHCH:J. SCOTT GAKQN'^R lJ.ATE 

: 0 2 / 1 A / S 3 


PL'PrCOe : TM:5 ‘■)URI>!)UT i nf lstipatps the signal fkcm the 

LCCAL HEAN AM’ VAfUA'.'Cl. . 

ORIGINAL PAGE IS 
D£ POOR QUALITY. 


NAME 

\ 

PARA 

TYPc 

''ETtR DEFINITION 

\ r.LASSX range \ DESCRIPTION 

ZME an 

\R 


\R 

\ 

\LCCAL mean 

VARZ 

\R 


\It 

\ 

\L1CAL variance 

Z 

\R 


\R 

\ 

WALUt OF CENTFH PIXEL 

XES TMT 

\R 


\W 

\ 

\EST1MA7E FOR SIGNAL 


\ 


\ 

\ 

V 


\ 


\ 

\ 

\ 


\ 


\ 

\ 

\ 


\ 


\ 

\ 

\ 


\ 


\ 

\ 

\ 


\ 


\ 

\ 

\ 


\ 


\ 

\ 

\ 


\ 


\ 

\ 

\ 


\ 


\ 

\ 

\ 


\ 


\ 

\ 

\ 



NON- 

L JCAL 

VAR I ABLE S 


\ 


\ 

\ 

\ 


\ 


\ 

\ 

\ 


\ 


\ 

\ 

\ 


\ 


\ 


\ 


SUBROUTINES REGUIREC 
^AME \ OESCRIPTICN 


\ 

\ 

\ 

\ 

\ 

\ 

\ 


SUBROUTINE ESTMTE (ZMEAN, VARZ. Z. NUHLKSt XESTMT) 


REAL ZMEAN, VARVf VARZt Zt XESTMT, XM£AN, VARX. K 
REAL NUHLKS 


SOURCE - - SPtCxLH ANALYSIS AND SMOOTHING IF SY,S'THF‘'IC 
APERTURE RAOAH IMAGES 

JCNG-SEN LEE CC.''PUTER GRAPHICS AND IMAGE f’f ESSING 1 7 t 2 A- 3 2 ( I S 6 I ) 


oooooonon 


C jlCNAl MEAN * IVAnF. MEAN / NCISF «EA.\ 

C 4SSLVE NCISF HFAN * I 
C 

XHEan » ZPPAf; 

c 
c 

VARV » l.C / NJMLKS 
C 

c 

VARX » (VARZ ♦ ZMEANO/.MEAN J / (VARV ♦ I.O) - X K F A X E A N 


I CAN ae LINEARIZFC 9Y THE Fl-^ST OkCER TAYLTR SiRIGS 
EXPANSION AROUT (XMCAN, VPEAN) 

2 * VHEANOX ♦ XMtAN * (V - VHEAM 

FRCM THIS AN ESTIMATE FOR X IS DEVELOPED 


K a VAPX / (Xf'CAN « XMtAN « VARV ♦ VARX) 
C 

. . XESTKT * XMEAN ♦ K # (2 - XMEAN) 

C 

c 

RETURN 

END 

EOF .. 

7 


'•Aog n 
POOR QVALny 


V. 


n o o o n n 


C PUCCRAM 3L1TC : KCI SC FILTERS 


LNIVERSITY CF KANSAS T t L E C CM R L M C A T 1 C N S ANO InFCRMAII 


A fc 


P CF . 


PRCCRAM NAMFiADPFLT 


AUTHCR:J. SCCTT CVUnN"R OATC:C2' 


PURPCSE : this IS THE MAINLINE FCR THE ADAPTIVE FILTER 

P.CUTINC XHICH USES A FILTER 

WINCCH WITH AN EXPONENTIAL WEIGHTINC FLNCTION. 


DRIUINAC PA 09 la 
DK POOR QUAUTY 


NAPE 

\ 

PARAMETER OEFINITICN 
TYPE V CLASSV RANGE \ CFSCRIPllON 


\ 


\ 

V 

\ 


\ 


V 

V 

V 


\ 


\ 

\ 

V 


\ 


V 

\ 

V 


\ 


V 

V 

\ 


\ 


\ 

V 

V 


\ 


V 

\ 

\ 


\ 


V 

V 

\ 


\ 


V 

V 

V 


\ 


\ 

V 

\ 


\ 


V 

V 

V 


\ 


\ 

\ 

V 


\ 


V 

V 

V 


\ 


V 

V 

V 



NON 

-LOCAL VARIABLES 


\ 


\ 

\ 

\ 


\ 


V 

V 

V 


\ 


\ 

V 

V 

•• 

\ 


\ 

V 

V 



SUBROUTINES RECUIREC 

NAPE 

\ 



OE SCR I PT ICN 

' 

AOPSUB 

VOCES the 

PROCESS ING — 

CALLEO WITH VARIABLE PARAMETERS 


V 






V 






V 






V 






V 






V 



• 



prcgrah acpflt 


INTEGER MAXSI2, MXWSIZt 3UFCUEt BUFKNOt MXFLTSi RCFFWN 


BLFCUE * MAXSIZE « M A X W I UCOH S I Z E 
BUFWNC = MAXW 1 NCCHS I 7.E * PAX W I NOCHS I Z E 
OL'FFWN IS THE BUFFER FOR TriE FILTERS ARRAY 


PARAMETER (PAXSIZ= 512 ) 


c 

c 

c 

c 

c 

6C1 

5 10 

700 

C 

C 

602 

C 

603 


C 

C 

610 


705 

C 


615 

C ' 

1 2 
625 


710 

C 

637 


reu 
:rR 
s 1 r R 


( X S S I Z « l 5 ) 

(vxFLTS»25) 

( GLF'jljF 7«MXHSI 2 » 

( = I 2<‘HXU<i I Z ) 

I 8UrriiN«KXFL TS c KXWSIZ 


0 “X\iSIZ) 


we quautv 


IMFvlEr' RECS. FLTNUM 

IMEGE'». INJ^'AxSiZ). CUKMAXSI Z ) . SlZEt CUlSIZt wNOSlZ 
IMEOEK FLThSl l./XFLTS) 


REAL tUEllE ( ilLFC’JS ) . W 1 NDH W ( 6UF WNC J i ^ l tr S ( ^UF F WM 
REAL NURLKS. NRAX, NMIA, C 1 S T ( 3U F WN C » . EOLiRfcS 

IMECEr.Ol lf»F.<*M( I a ) . CUTFNM(l8>t CUlZFN(ld) 


INTEGER TTY |Ki TTYUUT t I NFC fOuTFCi CUT2FC 


DATA TTYIM, TTYOUT /15. 16/ 


WRITS ( TTYCliTtCOl ) 

FCRFAT (IX, 'ENTER TFt FILENAHE FCR INPUT (HLST 6E AN CL C FILE)') 
READ (TTYIN.5I0) INFNM 
FCRXAT (13AI) 

WRITE (TTYOUT, 7C0) INFNH 
FCRFAT ( IX, IdAl ) 


WRITE (TTYOUT, 602) 

FCRRAT (IX, 'ENTER THE FILENAME FOR OUTPUT (MUST 6G A NEW FILE)') 
READ (TTYIN,510) OUTFNM 
WRITE (TTYCUT,700) UUTFNM 

WRITE (TTYOUT, o03) 

FORMAT ( IX, 'FILENAME FUR FILTER OUTPUT (MUST 3E A NEW FILE)') 
READ (TTYIN, 510) CUT2FN • 

WRITE (TTYCUT,700) 0UT2FN 


WRITE (TTYCUT,619) 

FORMAT (IX, 'ENTER THE SIZE OF THE INPUT IMAGE ') 

READ (TTYIN,*) SIZE 
WRITE (TTYOUT, 705) SIZE 
FORMAT ( IX, U ) 

IF (SIZE .LE. MAXSIZ) GOTO 12 
WRITE (TTYOUT, 615) MAXSIZ 

FORMAT (IX, '« **ERROR-- THE MAXIMUM SIZE » '.15) 
GOTO lOlO 

WRITE (TTYOUT, 625) 

FORMAT (IX, 'ENTER THE NUMBER CF LOCKS ') 

READ (TTYIN,*) NUMLKS 
WRITS (TTYOUT, 710) NUfaxS 
FORMAT (1:<,F6.3) 

WRITE (TTYOUT, 637) 

FORMAT (IX, 'ENTER THE ECUI VALENT RESOLUTION ') 

READ (TTYIN,*) EOUPES 
WRITE (TIYCUT.TIO) ECURtS 


ii. >iLUTt6^j) 

FCf'f'M (IXt’ENira THt CF L*:CKS •) 

RE&C JHYIN.O) ;m| n 
WPIT fc (TTVUUT.MO) NM|\ 

C 

W?M0 ( nVJLT.o^.?) 

647 FCHFAT (ir, 'ENTER THE MflXlKUK ^UM'JE9 CF LOCKS •) 

R E A C ( T 1 Y I \ t <:• ) 'i « A V 
WRITE ( TTYCL'T,7 10» NMAX 
C 

WRITE ( TTYCUT,64‘,) 

649 FCRf'AT nx, 'ENTER THE NL'MOER OF FILTERS *> 

REAC (TTYIN.Ol FLTNUH 
WRITE (TTYCUT.705I FLIKUM 
C 

c 

WRITE (TTYCLT.62M 

627 FCRHAT IIX, 'ENTER THE SIZE OF THE FILTER HINCOVC. 
£ /lOX.'CThlS PARAMETER MUST QE OOC ) ' > 

REAR (TTYINt «) WNCSIZ 
WRITE (TTY0UT.705) WNCSIZ 
C 




IF (MCC( WNCSIZ, 2) .NE. CJ GOTO 20 
WNDSIZ « WNCSIZ ♦ I 
WRITE (TTY0CT.632) WNCSIZ 

632 FCPHAT (IX, 'THAT IS NOT AM ODD NUMBCK. 1 WILL USE ',12,' INSTFAC) 
C 

2C IF (WNCSIZ .L5. MXWSIZ) GOTO 13 
WRITF. (TTYCUT,641» MXWSIZ 

641 FCPYAT (1X,'9 ERROR - - MAXIMUM WINOCW SIZE ■ *,15) 

GCTC lOlC 
C 

13 OLTSIZ»SIZE- WNDSIZ ♦ 1 
C 

c 

c 

C CALL SUURCUTINE TO 00 THE WORK 
C 

CALL ACPSUli (INFNM, JUTFMM, CUT2FN, CUEUE, WIMCCW, FLTRS, OIST, 

C IN, CUT, FLTHST, size, FLTNUM, RECS, WNDSIZ, CUTSIZ, NLMLKS, 

C NMIN, NMAX, ECURES, TTYOUT) 

C 

c 

-C 

WRITE (TTYOUT, 6301 

'.6 30 FORMAT (IXt’RO'X ALL CCNE 0 « «•) 

WRITE (TTYOUT, 6401 OUTSIZ, RECS 

'*640 FORMAT (IX, 'THE OUTPUT IMAGE IS 'tlS,' WORCS BY ’.IS,* RECCRCSM 
C 

C . 

COTC lOlC 
C 

c 

c 

ICIC STOP 
END 

EOF.. 

? 


*< a 



Meg n 

fOOR QUAUTY 


c 

c 

r 

LN IVFR S 1 T Y 

cr 

K At.*: a:; 

T FL rr C mmuM CAT 1 ONS ANO INFCRMATICN SCIENCES LAB 

L 

c 

FRCGR AM 

SU 

TE : 

NCI SL 

filters 

REF 

. i : 

c 

r 

PRLGRAM 

NAME : AOPSU". 

Al 

UTHC'C.I. SCCTT GSRONER OATE:T/A/“3 

U 

c 

PWRFCSE 

• 

• 

THl 

S IS TME SUBROUTINE TO PERECtm 

c 

THE ACTL'AL 

FRCCESSI NG 

FC.t Tht 

APART IVE 

WEI NER 

c 

c 

r 

F ILTLP 

RUUT 

1 NE. 





w 

c 




PAP AMETfc R 

DEFINITION 


c 

f 

NAFE 



\ TYPE 

S CL AS 

SS RANGE 

S CESCRIPTICN 

c 

I NF NM 



\CH«11 

SR 

S 

S INPUT filename 

c 

CUTFNM 



\CM«IU 

SP. 

S 
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S INPUT eUFFER 
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SFILTER USAGE HISTOGRAM 
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SIZE 
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SR 
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FLT NUM 



\ I 
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SMJMTER OF FILTERS 

c 

REGS 
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SW 
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SMJMBER OF RECORCS IN IMAGE 

c 

WNOSl Z 



\ I 

SP 
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SSIZS OF LOCAL AREA WINCCW 

c 

CUTSIZ 



SI 

SR 

s 

SSIZ'" OF CUTPUT IMAGE 

c 

NUMLKS 



\« 

SR 

s 

SNUHiE*. OF LCUK3 

c 

NM I N 



SR 

SR 

s 

SMIMMUM NUMBER CF LOCKS 

c 

N M A > 



SR 

SR 
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SMAXIMUM NUMBER CF LCCKS 

c 

ECUPE5 



SR 

SR 
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SECUIVSLENT RF.SCLUTION 
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TTYCUT 



SI 

SR 
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NAME 
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GENFLT 
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FILTERS AND PRINT THEM 

c 

LCSTAT 



SGETS 

LOCAL STATS AND FILLS A WINDOW ARRAY FRCM CUEUE 

c 

F ILTER 



SAPPLY 

THE PROPER FILTEF 

1 TO The LOCAL AREA 

c 

CPN 



SCPEN 

FILE AND 

ASSIGN F ILECODE 

c 

UCLCSE 



SCLCSt 

FILES OPENED WITF CPN 

c 

c 

F ILERR 



SREPORT TYPE OF FIIE ERROR 


SLORUUTINE ADPSUB (INFNPt CUTFNM, OLT2Fh, QUEUr, V.INCCW, FLIRS. 
C OIST, INf Ct'Tf FLTHST» SIZE. FLTMIM, RtCS. WNDSlZi CLiTSIZ. 

C NUPLKSt NMIN, NMAX, ECURES, TTYOUT ) 

C 

c 

c 

INTEGER SIZE. OUTSIZ. WNDSIZ 
C 

INTEGER FLTNL'M, FLT. ERRNUM 

INTEGER IN(SIZE). CUT(CUTSIZ). FLTHST (FLTNL’M 
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Pf.AI riT'<S JFlTUU-*, Wf-Clil/. '.<NCSIZ>. J I i” ( WN3S 1^: . V.NCS 1 2 ) 

RCAL WM:PT5. f /f*rAf , VAf’, .\UHlKSt IhSHLr.t XTSTl'T 

KFAL Mil.--., JUTSi), j-LTAN, Ny|-<. .WAX 

RFAL fCThPS. •.U'^PISt {CURFS 


IME^ERX*! INrN,«/(lrf). JUTr\»'(ia)t CUT?F\(15> 

LClilCAL 3RRt RN.:?PR 

IMEGFR TlYUUT.IN'FC, OUTFC , CUT2FC 




IMTIA;. IZE ZPEAN«7PEAr,/VARZ VIN ANC "AX AND TOTALS 


ThRS"I » ICCCCO 
THRSMX * -ICCCOO 
TTLTMN = 3.C 
TTLTSO = O.C 


RNGERR = .FALSE. 


ZERO hi ST ARRAY 

CC 5 FLTsI.FLTNUR 
Fl.THST(FLT) » 0 
CONTINUE 


CETERNINE THE NUMBER OF POINTS IN IMAGE (ASSoMEO SOUAPE) 
NUMPTS = OUTSIZ O OUTSIZ 

calculate the incremental NUMBER OF LOCKS 

TO BE USED IN INDEXING THE FILTERS 

• • 

OELTAN « (NMAX-N'MIN) / FLO A T < F L TNUM I 

OPEN FILES ANC CHECK FOR ERRCRS 

CALL CPN { INFCt INFNH# ‘OLCt *UNF*. tRRNUMf ERR) 

IF (ERR) GOTC 1009 

CALL CPN ( OUTFCf OUTFNM. 'N£W*» *UNF't ERRNUMt ERR) 

IF (ERR) GOTC 20C9 

CALL CPN ( 0UT2FCt CUT2FN, 'NEV.'* ‘FCR*. FRf>^lJ^/, .jkR) 
IF’ (ERR ) GOTC 3009 

HNOPTS * WM05IZ <» WNDSIZ 


v*)r 


) 


C oE-.tKAT^ 7h' FllTjr^S 
C 

CALL GCNFLT (ri. T-; 5. JIST, CLTNU^, WNC'; 

t LCLScS. C'.I?rC . ITYCU1 ) 

c 

C initialize Ti-E CIRL’.JLAK CL'l-JE 
C 

CC AC REC= I t W'iCS I £ 

c 

RcAO ( I NFC lEP < =A0U9 ) ( IM URD ) t WRC= 1 t S 
CC 3C 1<CRC= I t S I £f 
OUE'Jc ( REC I iCrC ) = IMVCRU) 

3C CCNTIMJt 

AC CCNTINLE 

C 
C 

C 8ECIN PRCCESSIUC 
C 

RECS * C 

CREC = 1 

C 

50 TMPCRC = ORcC 

C 

CC 90 STAPT=l,CUTSI£ 

C 

C ■ 'GET THE LCCAL STATISTICf. FCR THE AREA CEFIN'EC EY THE WINCCK 
C ANC FILL THE WINDOW ARRAY 
C 

CALL LCSTAT (CUE'JEt WlMCCW, TMPCRCt START, VJN3SIZ, WNCPTS, 

C SIZE, ZMEAN, VARZ, Z) 

C 

C 

EGTHRS a ZMEAN * ZMEAM t VARZ 
C 

C CHECK TC SEE IF LOCAL NUH0ER CF LCCKS IS OUT OF THE USER RANGE 
C 

IF (EGTHRS .LT. NMIN ,QR. EGTHRS .GT. NMAX) RNGERR = .TRUE. 
C 
C 
C 

C UPDATE RUNNING SUMS ECR DETERMINING THRESHOLD 
C 

TTLTMN = TTLTMN 4- EGTHR S/ N'UMPT S 
TTLTSD = TTLTSD ^ E GTHR S*E G TH R S / NUMP T S 
THRSMI = min (THRSMI, EGTHRS) 

THRSMX » MAX (THRSMX, EGTHRS) 

C 

C 

C - CALCULATE WHICH FILTER TO USE 
C 

FLT = FLTI,UM - I M T ( ( EGTHRS -NM I N ) / OELTAN) 

IF (FLT .LT. 1) FLT * I 
IF(FLT .GT. FLTNUM) FLT = FLTNUM 
C 

C UPDATE HISTOGRAM 
C 

FLTHST (FLT) = FLTHST (FLT) ♦ 1 
C 
C 

C PERFORM THE FILTERING 


1 , CEL tan. NUMIk*-, . . 


I Z c ) 




CPT S , 


c 


65 

r 

call riLTER (MMirw, FITTS. WNCSIZ. FlTNJM. FlT 
C ZMEAN. VARZ, Z> 

t 

c 

XESTMT » ’MTAf; 

pack a 
OE POOR QUAUTV 

c 


c 

c 

PLT TFE FILTCRL'D V..LUF 1/. THE CUTPjT '^UFFLI 
IF (X6STVT .C.T. 25-3.0) XESTmT = ?5*^.C 


CLT(SlflRT) « !M (XESTK7) 

C 

C 

90 CCNTIM;£ 

C 

C 

C WRITE This RECCRD AND UPDATE THE 9ECC''.C CCLNTEK 
C 

c 

WRITE (CUTFC) ( CU T ( 0 ) . WR C = I . CUT S 1 7. 1 

C DCN’T aCTHER TC CHECK FCK A WRITE ERRCk 
C 

RECS » RECS ♦ I 
C 

c 

C READ A NEW RECCRO I.'ITC THE QUEUE ANC UPDATE THE FRCNT-ENC PCINTEB 
C 

READ ( I.\FC t EN0 = 20CtERP.»A0C9 ) ( I N ( WR C » t '<R0 « I t S 1 ZE ) 

C 

DC IOC W0KC»ltSI7.E 

CL£UC( CRECtV.'JRO) = IMWORC) 

100 CCNTINLE 

c 

CREC « KC'J (CRECf WNOSIZ) ♦ I 
C 

CCTC 50 
C 

c 

c 

C WE’VE REACHED AN E.ND- CF -F I L E SC THE CUTPUT NUHUER CF 
C RECCPCS'WILL OE WN0SI2-1 LESS THAN THE NUMBER INPUT 
C WRAP IT ALL UP AND CUIT 
C 
C 

2CC CALL UCLCSE ( INFC ) 

CALL UCLCSE (OUTFC) 

C 

WRITE ( TTY0UTt'635» THRS H I . THRS H X , T T L TKN . TT LT SO- T T L TM N«T T L TMN 
635 FORMAT (IX.’XINlMUM NUMBER CF LCOKS * 'tF9.3, 
t /IX. ’MAXIMUM NUMBER CF LCOKS = ’.FS.3. 

C /IX. ’MEAN number CF LCOKS * ’.F9.3. 

C /IX ,’ STANDARD OEVIATICN CF NUMBER CF LCTKS = ’.F9.3) 

C 

IF (RNGERR) WRITE (TTY0UT.6A5) 

6<i5 FCRNAT (IX. ’» O*ERR0R — ENCCUNTEREO LCCAL NUMBER’ 

C /5X,’ CF LCCK5 rtHICri WERE OUTSIDE THE USER 3PEC1FIFC f.ANGc.’ 

C /5X.’ The FIRST DR LAST FILTERS WERE USED IN THESE AREAS’I 
C 
C 
C 


FILTER US Af. F 


P. * • 


VvPITL ir^ i-Ct 6SC) 

650 FCk fat J /• X . ’ - - - 
C USAOEM 

C 

CC 150 FLT a l,FLTNii« 

ViRITc ( CUT2FC t 70”) » F L T » F L C A T ( F L TH ST ( F L T ) 
7C0 FQRFAT ( I X t 1 3 t 6X , F 6 . 3 ) 

15C CONTINUE 

C 

C 

CCTC lOlC 


C 

C 

C 


ICC'! 

r 

CALL 

CCTC 

FILER’ 

ICIG 

I ttycut 

, INFNPt 

ERRNUK I 

C 

2005 

CAL L 
GOTO 

FI LERR 
lOlC 

( TTYOUT . 

CUTFNF t 

ERRNUF ) 

3C0<J 

r 

CALL 

GOTO 

FI L ERR 
1010 

( TTYOUT . 

CUTZFNf 

ERRNUH I 

5C09 

WRITE 

( 7TY0UT.66U) 




660 FCkFAT C1X,'« » 0 ERROR IN REACING INPUT 

c • 

ICIC RETURN 
END 

EOF.. 
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EADSUD \CCES THE PRCCESSING— CALLED WITH VARIABLE PARAMETERS 

\ 

\ 
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program EACFir 

INTEGER MAXSIZ. MXV»SlZf BUFCUCt QUFV<N0 , MXFLTSt OUFFNN 
INTEGER BUFHWV 


BLFCUC * MAXSIZE => H AX W U.CCW S I Z E 
BLFWNC = MAVU I NCUKS I ZE « MAXW I NOOWS I ZE 
BLFFWN IS the buffer FUR THE FILTERS ARRAY 



P At* A 

PAPi'-’ ** 
PARAiETER 
PAR A METT h 
PARAMETER 
P Ak AMET r fl 

parameter 


( VA> I ; = j 1 ,! ) 

« ”X'.»S W = 9) 

( ''XFLTS = 20J 

( ?-i;F.:ijL*MAaSI Zf.MXWSI Z ) 

( L‘l;K Wl JsMXV.'S I Z «MXWS I Z ) 

( bi;f-Kk>''»T;'f',xv,:>i zoMxkSi z ) 

( dUFr.JN=MxFLTS => MXkr.lZ « 


•'XWS I Z ) 


DWCINal pxg. « 


INTEGER RfCj. FLTNUM, MS K WNC ( EUF ) 

INTEGER SIZEt PujTSIZt VtMTSIZ. OU T ( M AX:, I Z ) 

INTEGER IMMAXSIZ) 

INTEGER FLTHSl CUFLT5) 

REAL CU£U£( SL'FCoE ) . W I NDflk ( eiJFWNC ) . FL TRS ( .1 UF F V.N ) 
PEAL MJMLKS. NMAX, .^MIN, C I ST ( 3UFViN'C » , EQLPEj 
REAL TMSHLD 

INTEGER*! INFNM(IO), OUTFNM(ia)t CUTZFfJdS) 
INTEGER T TV IN'. TTYOUT , I NFC t CUTFC. CUT2FC 
DATA TTY If., Tlfi'jT /1 5, 16/ 


kPITE (TTYCUT.6C1) 

FCRMAT IIX, 'ENTER THE FILENAME FOR INPUT (MUST BE AN CLD FILE)') 
RE AC (TTY IN, 510) INFfjM 
FCRMAY (13A1) 

WRITE (TTYOUT,7CO) INFNM 
FCRMAT (IX.iaAl) 


WRITE (TTYCUT,oC2) 

FCRMAT (IX, 'ENTER THE FILENAME FOR OUTPUT (MUST BE A NEW FILE)') 
REAC (TTYIN,510) CUTFMri 
WRITE (TTYUUT.7CJ) OUTFNM 

WRITE (TTYCUT,603) 

FCRMAT ( IX, 'F ILENAME FOR FILTER OUTPUT. (MUST BE A NEW FILE)') 
REAC (TTYIN.5’0) QUT2FN 
WRITE (TTYCUT,700) 0UT2FN 


WRITE (TTYCUT,6I0) 

FCRMAT (IX, 'ENTER THE SIZE OF THE INPUT IMAGE ') 

REAC (TTYIN,*) SIZE 
WRITE (TTYCUT,710) SIZE 
FCRMAT ( IX , I 5 ) 

IF (SIZE .LE. MAXSIZ) GOTO 12 
WRITE (TTYCUT,ei5) MAXSIZ 

FCRMAT (IX,'* ** ERROR- - THE MAXIMUM, SIZE = ',15) 

GOTO 1010 

WRITE (TTY0UT,625) 

FCRMAT (IX, 'ENTER THE NUMBER ’CF LOCKS ') 

REAC (TTYIN,*) NUMLKS 
WRITE (TTYCUT.720) NUMLKS 
FCRMAT (1X,F10.1) 


.ITj: ( 1TYt:UT,c-3( : 

537 hCRI-ftT THU tCUIVALtM P£SCH>Tr:N •) 

kft.Z ( T'^Yl? .»» • CUr’uS 

WRIT? ( TTY^LT f 72'JI cOURbS 
C 

WPITt; (TTYi.Lr.o2a) 

625 FCR'-AT (lX.*h'M*’P. THi EnCC THPcSHCLO VALLi*. ') 

P. EAC (TTY1‘<.«) THSHLf; 

WRnb (TTYnUT,720) THSHLU 
C 

WRITE ( TTYCUT . t,^b) 

6«S FCRFAT (l<, 'ENTER THE HIMKUK NUHBEP. CF LICKS ') 
REAC (TTYIN.O) NMIN 
WRITE (TTYCUT.720) mMIN 
C 

WRITE ( TTYCUT.^*./) 

5<i7 FCRPAT (IX, 'ENTER THE MAXIHUM NUMEER OF L'TCKS ') 
R EAC ( TT Y I \ ,« ) N -(AX 
WRITE (TTYCUT,720) NMAX 




^fJTv 


C 

WRITE (TTYLLT.oAtn 

549 FCRHAT (lY, 'ENTER THE NUMOtP. CF FILTERS ') 

REAC (TTY1N,«) FLTNUM 
WRITE (TTYCUT,710) FLTMJH. 

C . 

C 

WRITE (TTYGUT,527) 

527 FCRPAT (IX, 'ENTER THE SIZE CF THE FILTER WINDOW', 

C /10X,'(Th1S parameter must OE OCO'I 
READ ( fTYiN, #) WUnSIZ 
WRITE (TTYCUT,710) WNLSIl 
C 

IF (KCC(WNDSIZ,2) .NE. 0) GOTC 20 
HNOSIZ = HNCSIZ ♦ I 
WRITE (TTY0UT,532) WNOSIZ 

632 FCRKAT (IX, 'THAT IS NOT AM ODO NUMBER. I WILL USE 'tI2,' INSTEAC) 
C 

2C IF (WNCSI2 .LE. MXWSIZ) GOTO 13 
WRITE (TTYCLT,6A1) HXWSI7 

6A1 FORMAT ( I X , ' « «« ERROR - - MAX.IMUM WINCCW SIZE « ',I5» 

GCTO lOiO 
C 

13 CUT SI Z * SIZE - WNDSIZ ♦ 1 
C 

c 

c 

C CALL SUBROUTINE TO DO THE WORK 

CALL EACSUB (INFNM, JUTFNM, 0UT2FN, CUEUE, WlMOOW, MSKWNC, FLTRS, 
t OIST, IN, CUT, FLTHST, SIZE, FLTNUM, PECS, WNCSIZ, 

E CJUTSIZ. NUMLKS, THSHLO, NMIN, NMAX, E'JURES, TTYOUT) 

C 

C 

C 

WRITE (TTYCUT,o30) 

630 FCRMAT (IX, '« ALL DCNE «') 

WRITE (TTYCUT,540) OUTSIZ, REGS 

640 HCRHAT (IX, 'THE OUTPUT IMAGE IS ',15,' VICRCS BY ',15,' RECCRCS'I 
C 

c 


n o o 


I 

k 1 1 c 


ICIC STOP 
END 
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c 

KSKViNH 
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FLIPS 
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\R/W 
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\FILTERS ARRAY 
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c 
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FLThST 
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\ 
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c 

RECS 
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c 
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\R 
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\SIZr OF CUTFUT IMAGE 

c 

NUKLKS 
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\ 
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c 

ThSHLD 

\k 

\R 
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\ECGE THRESHOLD 

c 

NP I N 

\R 

\R 

\ 

\MIMMUM MJM.RER CF LCCKS 

c 

NPA X 

\R 

\R 

\ 

\MAX1MUM MJM.btR CF LCCKS 

c 

ECURES 

VR 

\R 

\ 

\EQU1VALENT RESCLLTICN 

c 

r 

TTY CUT 

\1 

\R 

\ 

\CUTPUT TO TERMINAL F ILCCCCE 

w 

c 

r 


NON-LOCAL 

VAR 1 AELES 

c 

r 


\ 

\ 

\ 

\ 

c 


SUDROUTINES RECUIREC 

c 

r 

NARC 

\ 

OESCR I PT ION 


c 

GENFLT 

\GE.MCRATE THE 

FILTERS 

AND PRINT THEM 

c 

CNHSKS 

\GENERATES MASKS 


c 

LCSTAT . 

\GETS 

LCCAL STATS AND 

FILLS A UINOQW ARRAY FRCM CLEUE 

c 

SUORS K 

\GETS 

A 3X3 LOCAL MEAN 

FROM wi icnv. 

c 

FNCEOG 

\GETS 

EDGE UPIFNTATION 

AND SFL'CTS A MIASK 

c 

FILTER 

\APP LY 

THE PROPER FILTER TO THE LOCAL AREA 

c 

egstat 

\GET THE cDGF 

STATISTICS 

c 

ECGFLT 

\CCMPUTLS LCCAL STATS 

CtTSIDF the edge 

c 

CPN 

\CPLN 

FILE AND ASSIGN 

F ILECGOE 

c 

UCLCSE 

\CLCS£ 

F I LES 

OPENED WITH CPN 

c 

c 

F ILERR 

\REPO.TT TYPE 

OF FILE ERROR 


SUORQUTINE EAOSun (INFNP. CUTFNM, OUT^TN. OLC-UC. *<INCCh, HSKWNCt 


S I i ' I r L I N 1 .' ■* t .» I' L ^ t ^ ^ u S I J 
t^UfieS. TTYCUT) 



t. f n?'n ... CjT, I L 1M1.T , 

C CL'TSI.’. KS. TMSHI D. N^IN, 


1 N T E Vi TP •: 1 .’ L . ['[} r S U . XN') S W 
INTEut-rv If.(SUE). :UT (CUTS I/. ) . N«0 

IMfulP Vi K3x:i ( ;j. 3 . ^ , SiiAMe A ( 3 . 3 ) , “S K w'<C ( o . WNO S 1 Z t V.KH S I Z » 


I M E Jcl- 
I M ECEL 
INTEGER 
INTSGEk 


FLTNU.^. FIT, CP.rNCiM 
FtTHST (FlTMJ'^) 

PtC, WCRD. :<LtS. CRtC. TMPOKC. START. 
>'SK\L'f., I, J, K 


ORIGINAL PAGE H 

1£,;0F poor quality 


RFAL -UEUE ( WM S 1 Z . S 1 ZE ) , W 1 M'CW ( GNDS I Z . V.NUS 1 ’ ) 

REAL FLIPS (FLTNJ.«1, WNCSIZ. WXCSIZ), D I ST ( WN3S I Z . GNCS 1 Z ) 
REAL GNCFTS. Z, ZMtAM, VARZ, NUKLKS. THSHLD, .XESTkT 
REAL THRSMX, THRSHI, TTLTPN, TTLTSO. OELTaN, MFIV, NKAX 
REAL EGTHBS, NUMPTS, tCURES* TTLhPN, N'AHPTS 

INTECEROl lNF'vM<l3), OLTFNy(l3). CUT2FN(ie» 

lcgical err. RNGERR 

INTEGER TTYOiJl . I NFC . QU TF G , CUT 2 FC 


initialize the 3X3 GRADIENT WINTCWS 


DATA ( ( ( MS K3X3( 1 . J ,Y. ) , I.l,3>, J«l,3), KM. A) 

C /l.l.l.O.O.O.-l.-l.-l, O.-l.-l.l.O.-l.i.l.C. 
C I.C.”l.l.C.~l.l.O.“l. I.l.C.I.C.'*’ItO.'*I.~l/ 


initialize 2XE ANOZMEAN/ VARZ KIN ANC MAX ANT TOTALS 


THRSMI « lOCCOO 
THRSMX ■ -ICCOOO 
TTLTHN « O.C 
TTLTSC » 0.0 
TTLHMN * 0.0 
NKHPTS » O.C 


RNGERR » .FALSE. 


zerc hist array 


CC 5 FLTal.FLTNUM 
FLTHST(FLT) s 0 
CONT I NUE 


DETERMINE THE NUM>3FP CF POINTS IN IMAGE (ASSUMED SQUARE) 
NUMPTS * nuTSIZ o OUTSI? 

CALCULATE THE INCREMENTAL NUMBER CF LOCKS 


T C i i. *J - w 


IKJtM.vC IML t^lLfthS 



CKT/.N » (^' ’MN) ! FI 'J.\T(PLTMJK) 
CF6N FILLS »'NC r.M'Cs E^RCKb 


P 

* 


call CPN < INFC. If.r-.M, •'!Lr*. ‘UNF*, l‘9.^1 

IF (C«P. I r.'jTC ino") 


CALL CPN I CUTPCt ODTFi.P, ‘NlW’i 'LNF', tHRNj-t, c«K) 
IF (c»:l) GJTC 2JCV 

CALL CPK ( ClT?FCi CJT2FN, 'PCK*. CRP) 

IF (c’RI GliTC jJO 

FIRST. CeNcHATE THE EO^F. TFKPLATES 
call GNPSKS (MSK'lNO. WNOSIZ) 


WNOPTS « WNCSIZ >> WN0SI2 


• GENERATE THE FILTERS 

CALL GENFLT IFLTRS. OIST. FLTNU*-. WNCSI.T. DFLTAN. NUFLKS. WNCPTS. 
t ECLRES. CUT2FC. TTYCUT ) 

IMTIALIZG the CIRC'JLAR CUE'JE 

DC ^0 REC»l.WrJDSIZ 

READ ( INFC . EM)»<.009, FRR«<.C09 ) ( I N ( WRC ) . hP D » 1 . S I Z E J 
DC 30 V.ORO-I.SIZE 
QUEUE! REC.wnR0)3iN( ^ORC) 

CONT I NUE 
CCNTINUE 


BEGIN PRCCESSING 

RECS = C 
OREC = 1 

THPCRC » CREC 

DC 90 STABT»l.OUTSIZ 

GET THE LCCAL STATISTICS FOR THE AREA DEFINED EY THE V.INDCU 
AND FILL THE WINDOW ARRAY 

CALL LCSTAT ( QUEUE . WINDOW. TMPQRC. START. WNOSIZ. WNDPTS. 
C SIZE. Zy.EAN. VARZ. Z) 


IF (VARZ .N£. 0.0) THEN 
EGTHRS » ZRCAN » ZH6AN / VARZ 
ELSE 



ukjuinal pace is 
rciMij;; ' or POOR quauty 

£NC If- 

c 

C C^•EC^; TC SE? IE irCAL MineL'i GF LGUKS is gut GF the L'SL« 'U'.'GF. 

c 

IF (LSn-^S .LT. ( .Oh. rCTHRS .GI. 

c I NMA''* .cl >MJ 1L KS ) ) «NGeo.l « .TRIJt. 

c 

c 

c 

C fPOaiF PCNMNC SU^'S FC'< OETr.BMINlNC THRESHOU; 

C 

TILTMN ■ TTLT^'N ♦ r. CTHRS/NU^'PTS 
TlLTiR « TUTSu ♦ E G TM R S»:>F C T H R S / NU hP T S 
TI-KSKI » PIN EGTMBSI 

THRSPX ■ :<AX IfHRSMX, EGTMRS) 

C 

C 

C DETERPINE IF AN EDGE EXISTS BASED ON THE LOCAL 
C NUHHE’. CF LOCKS 
C 

IF (EGTHRS .LE. NUHLKS - ThSHLD) GOTO AA 
TTLHPN « TTLH.<»N ♦ EGTHRS 
NMHPTS » NPHPTS l.O 
CCTQ 6C 

C • 

c 

C WE HAVE AN EDGE. SO PROCEED WITH THE EDGE FILTERING 
C 

C GET THE 3X3 SUBAREA LOCAL MEAN 
C 

tfi, CALL SU3MSK (WINDOW. SBAKEA. WNOSIZ) 

C 

C 

C FIND THE EDGE CP.IENTATION AND DETERMINE WHICH EDGE 
C TEMPLATE TO USE IN CALCULATING OUR NEW LOCAL MEAN 
C 

CALL FNCEDG (SBAREA. MSK3X3. msknuM) 

C 

C 

C CCMPUTE A NEW LOCAL MEAN AND VARIANCE USING THE APPROPRIATE 
C TEMPLATE 
C 

■ CALL EGSTAT (WINDOW. MSKWNO. WNOSIZ. MSKNUM, 

. .. C ZMEAN, VARZ , Z ) 

C 

IF (VARZ .NE. 0.0) THEN 
EGTHRS « ZMEAN * ZMEAN / VARZ 
. ELSE 

■ EGTHRS » NMAX 

END IF 
C 
C 
C 

C CALCULATE WHICH FILTERS TC USF. 

C 

FLT = INTI ( EGTHRS-NMIN) / DELTAN) ♦ I 
IF (FLT .LT. 1) FLT » 1 
IF (FLT .GT. FLTNUM) FLT » FLTNUH 
FLT a FLTNUM - FLT ♦ 1 



Ic 

c 



i; 

lie 


fPCAlt; 

fL THST ( FL1 > f L THSi ( FL ^ ) 

PFRFC^K the FCGF FMTCilINC 

CALL ECCFLT (v.UOOn. FLIPS. 
C FLT, . Sr NLIK, 2VCA ,, VAP Z , 


♦ I 


^ q[}m .. • 


'^SKWMC. Wf.'lSI Z . FLTNMM. 
Z ) 


GCTC 63 


THIS AREA IS HCHCGENECOS 

calculate which filter to use 

FLT « INTI ( ECTHRS-KHIN) / CELTAN) ♦ I 
IF (FLT .LT. l» FLT » I 
IF (FLT ,GT. FlTNUM) FLT » FLTNUM 
FLT ■ FLTNU.** - FLT ♦ I 

UPDATE HISTOGJUV 

FLThST (FLT) « FLThST (FLT) ♦ I 


PERFCRH THE F ILTt RI NG 

CALL FILTER (WINOCM. FLTRS. WNCSIZ. FLTN'jy. FLT, 
C ZMEAN, VARZ, Z) 


XESTHT « ZMEAN 

PUT THE FILTERtC VALUE IN THE OUTPUT SUFFER 

IF (XESTHT .CT. ZS5.0) XFSTP.r « 255.0 . 
OUT(START) » INT (XESTMT) 


CCNT I NUE 

WRITE THIS RECORD AND UPDATE THE RECCRC COUNTER 
WRITE (OUTFC) ( OU T ( WRO ) ,WR C» 1 . CUT S I Z ) 

RECS a RECS ♦ I 


READ A NEW RECORD INTO THE QUEUE AND UPDATE THE FRCNT-ENC FCINTER 
READ ( I NFC . ENUaZOO.ERRaAOCR ) ( I N ( WRC ) . WRO* I , S I ZE ) 

CO iio wcrd=i.si:e 

CUFUE(CREC.WORn) » IN(WCRC) 

CCNT I NUE 


♦ i 


r 


T- 


c 


635 

• • 

4 • 


•’ c 

* ^ 

,-637 

• 638 
C 

: 645 
f * 

! c 
c 
. c 

• 650 


C 


t-7C0 
■ 15C 
• C 
C 

V 

c 

c 


catc » ^'CC { .'r-r, 

CCTC 50 








wE'Vt REAChzr a;; f'jp-rr-riL" sc thl cui-'ct cr 

RECC7CS WILL 8 6 .iNUSlZ-l LrSS THAN THf. IfJPUT 

WRAP IT ALL UP A IC. CLIT 

CALL L'CLCSL ( Iflf C » ‘ 

call UCLCSf (CUT?C) 

WRlTfc- ( TTYCUT.633) THRSPItTHRSMX.TTLT'^N. TTLTST-TTLT’^N'STTLTyN 
FCKNAT (lX,*yiMMUM NUHBER Cr LCGKS = • . F ■? . ? t 
f. /IX, • VAXIRUF NUMiiEP. CF LCOiCS = ',F5.3, 

£ /IX.'yeAN NUMBER CF LCOKS - 'fF9.3. 

£ / IX standard deviation CF NUMBER CF LCIKS = ',F9.3) 


V.PITE (TTYCUT.037) TT L HHN/ N MhP T S 

FORMAT ( IX, •HCMITGEnEO'JS AREA MEAN « *,F9.3) 

WRITE CTTYOUT,o38) MMHP T S* 1 00 . 0 / NUMP T S 

FORMAT (IX,F9.3,* PEROENT CF THE IMAGE WAS HCMOGf N EC L S * ) 

IF (RNGERR) WRITE (TTY0UT,t45) 

FOFMAT (IX,'R <1 r p C R -- ENCOUNTEPEO LCCAL NUMBER' 

£ /5X,* OF LOCKS WHICH WERE C'LTSIDE THE USER SPECIF I EC RANGE.' 
£ /5X,' THE FIRST OR LAST FILTERS WERE USED IN THESE AREAS') 


WRITE (CUT2FC, 650) 

format </5X,'- - - FIL TER USAGE - - - ' , / 1 X , ' F I L T £ R A' 
1 . 5X,'i USAGE') 

CC 150 FLT * UFLTNUM 

WRITE (QUT2FC,70C) F L T , F L 0 A T ( F L TH ST ( F L T ) ) « 1 CO . C/ Nu M P T S 
format < IX , I 3, 6X , F6 . 3 ) 

CCNTINLE 


CCTC lOlC 


vC 






i 1009 

CALL 

F 1 LlfRR 

(TTYOUT, 

I N'FMM, 

ERRNUM) 


GCTC 

10 10 










1 2CC9 

call 

FILERR 

( TT YO'JT , 

CUTFN’Mi , 

C-RR NUM ) 

1 . 

r 

GCTC 

ICl C 




f*3CC9 

call 

FI LEPR 

( T fYOUT , 

O'JTZFN , 

ERRNUM 1 

1 .. 

GGTC 

IC 1 C 





i .<rOC9 WRITS I TTYCUT.oA.T) 

560 FORMAT ( 1 X,'« O .HRCF I r. RLACInG INPUT IMAO': * !.'• «•) 

■c 
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c 

c 

c 


c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

r 

c 

r 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


L.\ IVfk I f Y C 

iS «• a •• 1 

1 tlECC 

p L 

ICATICNS A,\’) INFCRPATICN ' 

^L 1 cNC'" i LA e 

pajCnA." ^lJTE : 

r.f i .V 

l‘ 1 1 Tt F 

«' 

J 

RtF. *» : 


P F C E R A M N A K t : .* L.' 'TIT 


-l‘iTHCR:.l. 5CCTT < ‘RPNtK 

Z2/ 

FLRFCGL : THIS E-jR" 

.'k;t j ' E 

COPP 

ltec ti-e local TM1:.TICG 

FCR 

A K INlTCV, If YLR 

APPL Y 5 

’.C A f ILTLR 

WHICH PLST 


paEVlCLSLY HAVE 

J f. 1 

N'l.R PALI 

; c . 




PAR A.'IETE R 

OEF 1\1T ICN 


NAPE 

\ TYPt 

\ CLASSY 

range \ CiSCRIPTIGN 


W INCCU 

\R 

\R 

\ 

\LCCAL AREA NlNCCki 


FLTR5 

\K 

\K 

\ 

YFILTtRS AFRAY 


PSKWN': 

\ I 

Vti 

\ 

\£2GE TEPPLATES 


N.N0SI 7 

\ I 

\R 

\ 

\SIZE )F V«I.\CijW 


nlpflt 

\ I 

\P 

\ 

\mipte: cf filters 


FLT 

\ 1 

\R 

\ 

VMIP'JL’ CF filter 

TC LiSc 

PSKNJP 

\ i 

\K 

\ 

\MJP1ER. CF PASK TC 

USE 

i >' L A M 

\F. 

\.i 

\ 

\LCCYL PFaW 


VAR Z 

\R 

\w 

\ 

\L"!CAL VARIaNC*" 


Z 

\9 

\W 

\ 

WAulL CF CENTER PIXEL 


\ 

\ 

\ 

\ 



\ 

S 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

. \ 

\ 



\ 

\ 

\ 

\ 



NON-Lt.T.A 1 

. V A R I A G 1 E S 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



\ 

\ 

\ 

\ 



SL'-i^lOUT If.ES RECUIREC 
NAK£ \ OESCP. IPTION 


\ 

\ 

\ 

\ 

SL'ERlUT 1 t.c ECCriT (,<INCCiW, FLTRSt i-'SKy4’';j. WN'lSl't 
C NLHFLT, FLT, MSKNUM. ZMEANt VAp.Z, L) 

IMEoEit UN'CSIZ. FLT, RFC, UC-IC, NtiRFLT 

INTEGER .-ISKjJhO ( H , WNOS I Z . KNC S 1 Z ) , MSKNUK. FlTNHX 

REAL FLTTS(Mj>'FLT,WNDSIZ,WNCSIZ) 

REAL V, I NOCV, ( kflDS I Z , WNCS 1 Z > ■ 

real ZREAN, VAPZr Z, TCTALy, TCTALV, FLT7AL 



TCTiL*' * n.C 

TCTALV » .r 

uc 2C -Jtc * I . .. . V. i; 

cc 10 = ll■^\C.>Iz 

I F I MS K' ( t-SKM. ' , f \.0^. C ) . 'EO.O ) TH'^N 

fl T..t :< «= mjmi l t - fl t » i 

fcL it 

FLlf.CA = rLT 
ENC IP 

FLTVAL » WllOTi'K F.L fViCk'' > O FLTKS ( F L T N )X t ftrC . 0 ) 

TCTALM ■* TCItL.* ♦ CLIVAL 
TCT'-LV » TCIMV ♦ HTVAL « FLIVAL 
CCM 1 1.MC 
ccr.T inlf: 

ZVEA'x s TOTALM 

VARZ = TCTALV - ’ /!; AN* 7.MF A tx 

RETURN 

ENC 




1 

I 

1 

1 

1 

1 

T 

I 

r 

i: 

i. 





Lf%IVL'ASlT^ c? Kd»»>J^, T r L f(. f'.-'f'l. M r. A T J ^N*i /.r. I A r' C> K A T I N iLICNL'C LAE 


PPZCtiAf'. '■.LITc : I.Cl^i riLTtf»;, f . •< : 


FRCGUA»' NAPfc : *•. I CFL T Ai)TI'C'’:J. SCLTI CMinKc.it (; *. T F : / 1 7 / » 3 


PLRFCSC : THIS I'.Ct.TINC PjRHG.TKj LFE’G G I G ’A FlLTcR. 


P 43 AKGTE F AEF INI n CN 

NAKC \ TYPJ \ CL/SS\ RANGE \ CGSCRIPTIGN 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


NON-LOCAL VARIACLcS 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


SUBROUTINES RCCUIRED 
NAVE \ OtSCRIPTICN 


SIGSue VOCES THE PROCESSING 

V 
\ 

V 
\ 

V 

V 


FRCGRAV SIGFLT 


CALLED WITH VARIAdLE PARAMETERS 


INTEGER MAX31Z. MXwSIZt PUFQIJE 


BLFCUE = VAXSIZC « VAXW I NCCW S I ZF 


t 

« 


PAkA-^ETF.R (VAXSIZ = 512) 
PARAMETER (VXRSIZ=15) 


1 


( W ^ • u . • 


1 

I 




I 


720 

C 

c 


528 


C 

c 


627 


t 

r 




•'* * i I w • • ' t .X l ( j I /. ) 


1 N T : c E R =< I c j 

IMECEK if.c-’AASi ; ) f \.7i ^ ) , sij;. L'ji'i’, ..\'c;i; 

»<r;.L a.:’.'- ( ) 

|’*:AL Na-LK*;. f 

I M P OE - I 1 N r ; M ( J J ) t L T F t. V( 1 ,i ) 

IM=:£^ TTYIS,TTY.;UT.ltFC.':L!TFC 

DAT A 7 TV U- , TTV^'JT / 1 «i , 16/ 


■;v^> 


:i 


'Arf. 


/ 6 - 


WRIT£ (TTvCLT.tCU 

FCPf'AT (U.'LMFP THi: FILrNAyF FC^ INPbT (l-'LiT Pf CLC FILE)') 
READ (TTYIN.510) I.>IF,M 
FCR XAT ( I 3A I ) 

WRITE ( TTYCLiT. 70J) IMFNK 
FCR KAT I I 1 ) 


l.P I T£ ( T TYCL’T.oCZ ) 

FCHFAT (lX,‘tMER THE FILEAAHE FCR CLTPJT (--^LST dE A NES FILE)’) 
REa: (TTYlNfSlO) OUTFNfX 
WRITE (TTYCLT.7Cb) CUTF.'JM 


WRITE (TTYCt'T.MO) 

FCRHAT (IX.'t'XTER THE SIZE CF THE INPUT IXAGt ’) 

READ (T7YIiv,») SIZE 
writ:- (TTYCUT,710) 51/E 
FCRHAT ( IX, 10) 

IF (SIZE .L':. .iAXSIZ) GCTO 12 
WRITE ( TTYCuT, 610 ) fXAXSIZ 

FCRHAT (IX, •» X>-Y ERROR- - THE VAXMLH SIZE = ’,10) 
GOTO 1010 

WRITE (TTYOUT.625) 

FCRHAT (IX, 'ENTER THE NUMBER CF LCCKS ’) 

READ (TTYIN,«) :4UMLKS 
WRITE (TTYCUT,720) NUKLKS 
FCRHAT ( 1X,F10. 1) 


WRITE (TTYCUT,623) 

FCRHAT (IX,'ENTFR TH.: SIGMA THREShCLD ') 
READ ( T T Y I N , « ) K 
WRITE ( TTYL'UT, 720) X 


WRITE (TTYCLT,6Z7) 

FCRHAT (IX, ’ENTER THE SIZE OF THE FILTER -.MNCOW’, 
C ‘/lOX,MThIS I'ARAMETER must ‘JC ODD)’) 

RE AC (TTY IN, *) WNCilZ 
WRITE (TTYCUT.7I0) VJNCSIZ 



IF ( .' L. J 1 h > 1 S 1 Z . 

2) .Ri . 

C) cal': 2c DWUiv3c 


wr.osiz •» •'^CSI.’ 

♦ 1 

POOR 


WRIT; ( T TYCLT . . 

32) W'iOS 

1 z 

612 

FCRFAT ( IX. • IHA 

I 16 .Cl 

.*.M CU3 M.MJCR. I lILL uSf. 

c 

2C 

IP (t.\.ySIZ *L;. 

I'c.is; 1 ) 

(.OTC 13 


WRITE 1 ITYCLT, 

6<.i) (y... 

SI 2 

6<. 1 
r 

F C 1 ► A T ( 1 / , • 
GCIC lOlC 

•> 0 ^ P 

P C K - - MAxra' MMCS 

U 

1 3 

CLTSI Z » SIZE - 

^\ns 1 Z 

♦ 1 


c 

u 

C CALL SCTiPCUTINd Tf*. O'J THF 
C 

CALL SiCS'Jc (I\FjM, OUTFM*', CUEUCt INi ?LiT, 

C SIZE. necs. V.-SJOSIZ. CL'TSIZ. NUVLKS. K. TTYCUT) 
C 
C 
C 

V.RITE ( TTYOUT.f>3'J) 

6 30 FCR ►AT ( I X , • « 0 •> ALL C C N E « » «•) 

WRITE (ITYCUT.tAOl JUTSIZ. REC3 
6A0 FCRHAT (IX, ‘THE OUTPUT IMAGE IS *,IS,' WCTCS lY 
C • 

C 

GCTC 1010 
C 
C 

c 

ICIC 5TC? 

ENC 

EOF.. 

*> 


" PAojf n 
QI^AUTV 

• . I ? , • I N ‘i 1 P A C • I 


b I F * ’.15) 


15,' RECCPCS • ) 


n n o o o o o o n n r» o o o o n o r* t'l o ,•* r-v 


f. 

C 

c 

r 


C 

c 

r 

C 

C 

c 

r 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


dtoginai: page is 

DE POOR QUALITY 

lmversity ce 

K/.NS.*.-, TrLcGf-.MiaMCAT IONS A.\T |NrC<V*.nON SCIENL'S LAi 

PRCCRAP slit? : 

\ris. I- 

I L Tt Ri 

K r-F 

. •» : 

PRCGf'A.y NAi-E:S! 

SSL i 

A'JTHfrJ. S 

COTI gsRHNER 'JATF:.V^/K2 

PURPOSE ; Tin 

S IS T n. 

SU.IRIJUT 1 NE TO P 

ERPC- y 

THE actual PKrCFSSIVG F 

O LGC* 

S SIGKA FILTtR. 


r A { 

: R TEFISITICN 


NA^E 

\ typ: 

S CLASS 

S RANGE 

S C'SCRIPTICN 

INFNf 

NU*!' 1 J 

SR 

S 

sinp‘jt filename 

llUTF.'.y 

\CH« 1 J 

Si* 

S 

SCUTPUr «=lLtNAMu 

CLELE 

\R 

SK 

S 

S IMAGE CATA GUELE 

I N 

\ 1 

\M 

S 

sirgpuT eiiF«=E® 

CUT 

\ I 

SiJ 

s 

SCUTPUT *^UFFER 

SIZE 

\ 1 

SK 

s 

SI PAGE SIZE 

REGS 

\ I 

S.J 

s 

SN'IHiER Cr RECTK':S IN IMAG‘D 

WNCSIZ 

\ 1 

SK 

s 

SSlZr. OF LP.UAL ARcA KINLCk 

CUT SIZ 

\ 1 

SR 

s 

SSIZE nr- CUTPjT image 

NCR LKS 

SR 

sy 

s 

SNUMBE" -OF LCCNS 

K 

SR 

SR 

s 

SSIGMa THREShDLO 

TTYCUT 

S 1 

SK 

s 

SCUTPUT TC TERMINAL FILcCCOE 


S 

S 

s 

S 


S 

S 

s 

s 


S 

S 

s 

s 


S 

S 

s 

s 


S 

\ 

s 

s 


NON'- 

L.'JCAL VARIABLES 



S 

S 

s 

\ 


SUiJROUT INES 

RECU I REC 


NAKC 

S 

CESCR 1 PT ICN 


CPN 

SCP: S FILE AND 

ASSIGN F ILECCOc 

UCLCSE 

SCLUSt f 

ILES OPENED WITH CPN 

F ILERR 

SRE.'CRT 

TYPE OF 

FILE ERROR 

SUHR3UT I NE 

!j 4 G S J (I 

1 1 NFNM, 

CUTFNM. 

CUSUE. IN. CUT, 


SIZLt Rr.St WNDSIZ. CUTSIZf ^U^'LKSt K, TTYCUT) 


IMGO^R TUTSU. W\DSIZ 

INTEGER LRRNJP 

INTEGER IMSIZEI. CuT(CL'TSIZ) 

INTEGER REC. -iCRC, RECS. CREC, KICRtC. T'PCRC. START. IE»R 
INTEGER CVtCRC, MIGKilC. URG. WREC. HXTORC. INTCX 


real OUEUE ( V.NDS 1 Z . S I Zb ) 

REAL NURLKS, XESTRT, R, WSRCNT 


J' *- ♦. i L 't - L 1 M t, ) , 

VAL'Jc. S')?.*’ 


L-’h^K, UPPl’, L!.'Kt l,V ^ 
ToTalk, cent'ii, usft n 


1 N T J i /• * ! I F '' i 1 ) . i ‘ l; I r r '^ ( ! d J 

0/; 

T 1 V .Ll . ! '-PC , ''Ij T» c Qs 

imtiil?. uf pc a aj: l }v.'La t. arrays 
data .UI'* /.«?•.. .121. .<>07. .273,. 367. .',c/ 

DATA UPRU*^ '’•'>'»»^‘*7«5,2.A08,2.1<»l,l.<j/,5,l.7Cw/ 
T*-= FCLLlV,I-jC data ‘TArrKEMS A<^ 6 FGR A CAl.SSlA.M 

AS' ,v,>T ic\ 

c 4 " A L w ^ L I « / c . 0 . C . 0 , C . 'J , 0 . tj , . I )J , . 3 7 / 

C4TA UPRM-'I /3.C,2.^i,2.15,2.C,l.R2,1.63/ 

I A 3 ‘ ' NLi''t K s ♦ 0 . ; 

^ 1 1 • . c;< .G A ) tmsn 
F ( ^- CcX .lE. d) THEN 
• O'"* s 5 

- ^ -• 

t ° C?}' = c 
i ■ . If 


^^cri^AL 



c 

EAl! if 


• 


LCViCR = 

L WRI. I R 

1 

4 - 

C 

vjPPER = 

•JPRL I.-' 



CSRCf.T 

= C.C 

1 

c 



« 

c 




CPtA fills anc check FJP ERRCRS 


!■ ")LC*, *U.\F', ifRRNt.H, ERR) 

IF ( I R{. ) GO T C lOOy 


CALL CRN ( CL'TFC, OL'TFIJR, 'NCW*. *CNF*, ERRMj'^, ERR) 
IF (cRR ) GOTO 20C9 


imtialicf The circular CLLUE 
DC AO REC=1 ,V,.gCSI2 

READ (INFC,ERK=-ir!09) (Ii\(wRl;),ViP.D=l,Si;£) 
DC 3C WCRD=1,SI2E 
DULUc ( RFC ,i.‘CRG ) = IM WCRD) 

COM INUF 
CCK’T I KL'C 


fio oooofior* ooono ooiu 


PtCC'S.ll.L 


C 

c 


c 

‘ c 

c 

c 

c 


c 


c 

c 


7C 

c 


0 


R r c ■ 

C‘>t:C . 

t 


M C R • r. 

i 

. / j • : 

T “ P C ' C 

» 0 ; :c 


cr T3 

' T A R T . 1 

.OCT 31 ; 

TLT AL V 

» c. : 


SL’dPTS 

a 0.? 


M ICWRO 

a s TART Ml 

iPS 1 /./: 

CENTER 

a OLEUF 

( .>^1 OR iC . WI0V,RD ) 

LWR a 

LChER » 

CEMTEO 

L'PR * 

upper »:• 

CF'lTl R 

CC 90 

WREC a 

1 .WHCSI 2 

CC 7 0 

CWCP.C 

a START , STAR T-^WNDS I Z 


VMUL = CLcUi:(T <i>CftC.CWCKC) 

IF (VftLUc.GE.LHR V A L Lf. . L E . R ) 

TCTALK = TCTAL'1 ♦ VALUE 
SJ3PT3 = 5LiPT3 1 
EMJ IF 
COM I MIC 

T VPCRC a MuU ( T MpQ„c # '<\nS I 2 ) I 
CCM 1 \lE 



J 




IF (SUUPTS .LE. Ki THEN 
totals = CENTL=R ♦ OUELF ( M DR EC . K I OV.R C ♦ i ) 

NXTDRC = f'CD ( fM CREC . WNOS I Z ) ♦ I 

TCTALV * TCTALM QUEUE ( N X 7 C RC t K I CHkO ) ♦ CU tUE ( f,X T C RC i '■< 1 DkR C ♦ I ) 

XtSTyT a TCTALM / 4.0 


USECM a USOC.JT «■ 1.0 


ELSE 


XEST»-T a TCTAL‘' / SU'^PTS 
END IF 


PUT THE FILTERCC VALUE IN THE OUTPUT 8UFF£f< 


IF (XESTPT .GT. 255.0) XESTH7 * 255.0 
CUT(START) = INT (XESTKT) 


CONTI NUE 


WRITE THIS RECORD AMO UPCAT? THE RcCCRC CCLMTF.R 
WRITE (CUTFC) (CUTIC.IC) .W«C = l fOUTSIZ ) 

CCN’T fJCTHcR TC CHEC< FCP A WHITE ERROR 




p r L ' 


c RtAC i f.‘ i iM-: iiM Ar.c tP':.Ui. rcr.ips 

r 

sf/.c ( 1 NF^. . t;.:= jc . 1 1 *. ( whd » . i . : i :n 


L'C 10 ''. «cf:* i . ' ! ; j 

CLFiiC ( C •lit,.':' ) » 

ICO CCNTINl*; 

c 

CPLC » HCC ICR£C. * l 

Kii;5£C « *'cc (••io.<tC. viNcr. JZ) ♦ i 
r 

GCTC «'J 


^^GlhtAL Pa.'v 
POOR QUAuVy 


I C WE'VE REACHED dfv E;0-CF-F I LC SC THE CUTPlT CF 

C RECCF.es ViILL SE /NO'ilZ.-l LESS THAN THE NJ’10:K INPLT 
• C WRAP IT ALL JP A',0 OUT 


2CC CALL L'CLCSE ( I\FC ) 
call uclcse <cuirc) 


V.RITE ITTYCCT. »>50» 100. 0 » LSECM / L'.' A T ( OU T S I Z ••‘tJU T S 1 2 ) 

S5C FCRFAT (/5/.'T!iF SIG'.1A TiiRESmClC WAS APPLI^^O TC '.F^.2, 

C / lx . • P£RC~:iT OF THi IHi-.CE') 


CCTC ICIC 


IlCO^ CALL FILERR ITTYOUT, INFNF. ERR NUF) 

GCTC 1010 

I t 

2CC9 call FILCRR (TTYCUT, CuTF.NM, EURNUWl 
CCTC ICIO 
C 

|aC 09 WRITE ( TTYCUT. SSOl 

66C FCRFA7 (IX, •» O o LRRCR IN READING INPUT I^AGE «*•>•) 

IlClO RETURN 
• ENC 

SCF.. 


LNlVE*'Snv CF KANSAS T « I M CM »» 1. M C A T I T N S AN:> |NFL:!mATICN SCIENC'^S L A u 


pa c gram SLMt : KCl^JL MltLF*, 


P F F . -I : 


FRCCRAM NAHE:C^.PSK5 


AUTMC«:J. SCCTT CAwnNFR CA f F ; C 2 / H / “ 3 


PL'flPCSE : THIS IS A SL'JRCLTIAF TC liFNFkATF fl.F 

ECGt TEMPLATES FCP. AfJ AJAPTIVE COCE FILTEP. 


^^OOlhfAL PArr r* 

* '«OR QI'autv 


kAME 


PARAMETER OEF IN I T I CN 
\ TYPE \ CLASS\ RANGE \ 


DESCRIPTION 


NAME 


\ARRAY CF 3X3 MASKS 
\ARRAY CF EDGE TEMPLAIFS 
\SIZE OF SQUARE MlNCCkS 
\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


NON-LOCAL VARIABLES 


SUBROUTINES RECUIREO 
\ OESCRIPTICN 




SUBROUTINE CNMSKS IHSKwNC* UNCSI2) 


INTEGER WNCSIZ. LINE. CLM, START 
INTEGER MSKWNUIO.WNUSI ZfUNCSIZ) 


FILL THE EDGE TEMPLATES 



c 


1C 

c 

I 5 


2C 

C 

2S 


30 

C 

35 


^0 

c 

45 

• 

50 

C 

55 


60 

C 

65 


7C 

C 

75 


8C 

C 

• l‘?0 

'200 

C 

c 


EOF . 
? 




ST APT » I 
CC ?00 L1KE»1 . 

cc ISC CL“»i .’.jNOSi; 

IFJCLK .LT. GCir, 1C 

P$KV.N!:( I .L l\r . CL'^J « 1 

CCTt 15 

y S K i* N C « I I L I E I C L >' ) » 0 

IF (CLM .LF. WNUSIZ/2»l» CCTC 20 
HSKUNRI 5 . LI KE. CLM » ■ 0 
CCTC 25 

ySKliNC ( 5 tL INEt CLK I » I 

IF (CL.“ .GE. START) GOTO 30 
ySKWNOI 2tL I NL.CL •<) « 0 
CCTC 35 

y$KViNC( 2 iLlNE. CLP ) • 1 

IF (CLM ,GT. START) GOTO 40 
.MSKWNOC 6tL INEt CLM) « I 
GCIC 45 

ySKWNC(6tL INE, CLM) « 0 

IF (LINE .LF. WNDSIZ/2tl) GCTC 5C 
MSKV.NO( 3.L INS, CLM) « 0 
CCTC 55 

MSKV.NC( 3.L INE.CLM) » I 

IF (LINE .Lc. WNDSIZ/2) GQTC 6C 
MSKKNO( 7 ,L I NE. CLM) » I 
GCTC 65 

ySKWNCI 7 ,L INE, CLM) s 0 

IF (CLM .LE. WNOS 1 2-START + l ) GCTC 70 
MSKHNC( 4 ,L IKE. CLM) * 0 
CCTC 75 

MSKWNC(4 ,L INE. CLK) « I 

• 

IF (CLM .LF. MNDS I Z-START ) GOTC 80 
yi5KWNC( 3.LINE. CLM) ■ 1 
GCTO 19C 

MSK),NC( a.L INE. CLM) « 0 
CONTINUE 

START a START ♦ I 
CONTINUE 




AL 




RETURN 

END 


c 

c 

r 
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i A\0 1^FC•1^AIIC^ SLILNC'S 1, A d 

c 

r 

FRCGRAi- SLME : JvCI^c 

filters ref 

. K : 

L 

c 

c 

prccf.am nape 

:LCST AT 

Ai,lTHCR:J. SCCTl GAitDN'R 0ATF:C2/ 13/32 

K- 

c 

PURPOSE ; 

this SUeUOUTlNb CQPPUTES THE LOCAL STATISTICS FCK A 

c 

wiNCCVi IN A 

C IRCULAR 

CUEUE AND FILLS A 

WINOOW ARRAY 

c 

FOR LATER pnCCESSi.NC. 


ORIGINAL PAGE 13 

c 

c 

c 

r 




QJE PCXDR QUALITY 

w 

c 


PARAMETER definition 


c 

r 

AAPE 

\ type 

\ r.LASSN RANGE 

\ C^-SCRIPTIOn 

w 

c 

CLEUE 

\R 

\R \ 

\DATA DUFUE 

c 

W INCCR 

\R 

\W \ 

\WIN...'< FILLEC FP.CF CUEUE 

c 

TMPCPC 

\ I 

\R/W \ 

\FRCNT OF OUEUF. 

c 

START 

S I 

\K \ 

\START CF WINDOW IN GLEUE 

c 

WSCSIZ 

\I 

\R \ 

\SIZE OF WINDOW 

c 

KNCPT S 

\R 

\R \ 

\NUP3E3 CF PCINTS IN WINPCW 

c 

SIZE 

\ 1 

\R \ 

\SIZF OF CUEUE 

c 

Z»'E AN 

\R 

\W \ 

\LOCAL mean 

c 

VAR 2 

\R 

\W \ 

\LCCAL VARIANC' 

c 

Z 

\P 

\W \ 

WALUE CF CENTER FUEL 

c 


S 

\ \ 

\ 

c 


\ 

\ \ 

\ 

c 


\ 

V \ 

\ 

c 

r 


\ 

\ \ 

\ 

c 

r 


NON-LOCAI VARIABLES 


V 

c 


\ 

S \ 

\ . 

c 


\ 

\ \ 

\ 

c 


\ 

\ ■ \ 

\ 

c 

r 


\ 

\ \ 

\ * 

c 


SUJIROUTINES REGUIREO 


c 

r 

NAHE 

\ 

OE SCR 1 PT ION 


c 


\ 



c 


\ 



c 

10 

\ 



c 


\ 



c 


\ 



c 


\ 



c 

r 


\ 




SUaROUTINc LCSTAT lOLlEUE. WINDOW, 

. THPQRC, START, WNCSIZ, 

r 

£ -WNDPTS, SIZE, 

ZMEAN, VARZ, Z) 


w 

INTEGER 

WREC, WWORCt OWCRC, START, 

1 WNCSIZ, SIZE, TMFCRC 


REAL Ct'EUE(KNOSl7.SI7fc»t WNCPTSt ZMEAN, VARZf Z 
REAL TCT4LH, T'lTALV, W I NO C W ( WNC S I Z t ViNCS 1 Z ) 


C 

C 


TOTALS' » 0.0 


onn oofNi nnnn 


TCTALV * 


C.C 


C 


1 c 


CC 2C ViiJCC* I . rtMJ*; I ? 
wwC(»3 * I 

CL 10 C..C!<J»ST \RT , ST ART ♦WNOS 1 Z- i 

tctalv = TCT.uM ♦ OULU. ( TKPCRC iCwcxrn 

TCTALV = TCTALV ♦ OUf: UE ( T P POR C t 0 CR C ) * « ’ 
h I NCLU ( riREC . WWCRU ) » CUE LE ( T MP UR C » CW CR O 
H.CRC «= WWCRC ♦ 1 

CCNT I Nue 



TVPCRC = '<CC (TMPORC, WNDSIZ) ♦ I 
CCfiT I NUE 


ZPEAN * TOTALP / W.NCPTS 
VAR2 » TCTALV / WNOPTS - ZKtAN«ZMEAN 

FINC the center PIXEL 

Z * W INCCWI WNCS I Z/ 2+ I . WNi)SIZ/2^1) 

C 

RETURN 

END 

ECF.. 

o 


DRic;r?^Ar pa(;e rs 

OF POOR QUALITY 


• ' / 


(. 

r. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


•C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 


c 


c 

c 

c 
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REF 

. It : 


PRCCt' AM 

NAME 

; SUBmSk 

/ 

M 

llTHCRiJ. Sr.CTT CARONEk OA I E : 

C2/ 1 3/ T3 

PURPCSR 

J 

THIS IS 

A SUBROUTINE TC GENERML A 3X3 UINCCR 

c= lccal 

KEANS FROM A 

UINUCU 

or CaTa PCINTS. 






parameter 

OtF INIT ICN 

• 

NAME . 


\ TYPE \ CLAS 

S\ RANGE 

\ description 


W INCCV, 


\R 

\R 

\ 

VUINRCf OF AREA OF 

INTEREST 

S3ARE A 


\ I 

\W 

\ 

\3X3 LOCAL KEANS 


WNCSI Z 


\ I 

\R 

\ 

\5IZE OF WINDCM 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 





NON-LOCAL 

VARIABLES 





\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ 




\ 

\ 

\ 

\ • 




\- 

\ 

\ 

\ 




SUBROUTINES RECUIRED 



NAKC- 


\ 

OESCR I PT ICN 


' 



\ 







\ 







\ 







\ 







\ 







\ 







\ 






SUBROUTINE SUDMSK (WINCOVJt S^AREA, V.NCSIZ) 


INTEGER WNCSIZ. SrtWSIZt OVHLAP, SUAR»;A(3f2) 
INTEGER STARTR, SVARTW. MR EC# MWCKC, KEC. WCRU 

REAL I, I NOCV, ( V,NrS I Z fWt^DS I Z ) 

REAL NCHPTS» TOTAL 


CALCULATE THt SUiUREA SUE >NC PVTLAP SIZE 


1 


c 


CVPL/»P » •^CC C.JNOSI 7 , 3 ) 

SPWSi; « WNC5l<!/3 ♦ L:VKLAP 



^LMPTS » S!1V,5I/ » SR.vSIJ 
STARTS = 1 
STARTW « 1 

CC <.0 Msec » 1.3 
CO 3C MWCRC » 1,3 
TCTAL « 0.0 

DC 2C SEC * STARTS, ST A R TR ♦ S 3 W S I Z - 1 
CC 1C WCRC = STARTU, S T A R T « ♦ S 3 R S 1 2 - I 
TCTAL « TCTAL ♦ Ulf.nuv^C RfC.WCSC) 

COM !NUE 
CCNT I NCR 

S RARE A ( PREC ,Hl,CRO ) » TCTAL / NUMPTS 
STaPTw s start W ♦ S'^^‘S1Z - CVRLAR 
CCNT INUE 
STARTW » 1 

STARTS a starts ♦ SUWSIZ - CVRLAP 
CCNT I NUE 


RETURN 

END 


ORIGINAL PAOK It 
DE £00R QUAUTY 




n n 


DRIGINAI. 

OF POOF QUALITY 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


L M V t K 5 

ITV OF 

TH. r '..OIMLM CAT ITN AN’' INFlTHATICN SUlNCES LAi 

PROGRAM 

SLITJ : nC!‘. 

f 1 1. T t R S ntF. i : 

PRCCR am 

name : F N C E C G 

AUTHCRrJ. SCOTT GARCN-R Q A T E : C 2 / 1 A / a 3 


PURPCSfe : THIS SU'’P.*JUT I fjf USfS A 3 X 1 r.rtAIjlc'.'T HASK TO 

CETE.'^^'I^E ThL tCGC OR 1 u M A 1 I ON AND THEN THRt SHLL L*S 
TC FIND WHICH S U)F OF TH£ C'TGC THF CLMER PlX-'L 
LltS. FRCV THIS A MASK NJMCEK IS SELLCIEC FOR SUhSECUCM 
CCGE FILTERING. 


PARAHETLP DEFIMTIGN 


C 

r 

NAPE 

\ TYPE 

\ 

CLASS\ RANGE 

N description 

L 

c 

SBAREA 

\ 1 

\R 

\ 

\ 3 XJ S-JbAREA MEANS 

C 

MSK 3 X 3 

\ I 

\R 

\ 

VGRADI^NT MASKS 

C 

PISKNL.P 

\ I 

\W 

\ 

XN'IPRER CF P.ASK TC USE 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 


\ 

\ 

\ 

\ 

c 



— 




C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


KON-LJCAL VARIABLES 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


\ 

\ 

\ 

\ 


NAPE 


SUfjROUT I.N 6 S RECUIREC 
DE SCR I PT I ON 


\ 

\ 

\ 

\ 

\ 

\ 

\ 


SUBROUTINE FNCEOG (SUAREAt MSK3X3» MSKNUK) 

INTEGER SBAREA(3t3). P SK3 X 3 ( 3 . 3 . A ) . HSKMJPt MAXVAl 
INTEGER CENTER# MSKNOX. ROW. CLM. SU^* 

INTEGER THRSHlf THRSH2 


CENTER = S8AREA(2t2) 
MAXVAL = -ICCCO 


c 


ORIGINAL PAGE IS 
DB £QOR QUALITY 


CC 50 HSKNCX « I . 

5LK « 0 

CC IC RCW = l.J 
CC IC U ^ = 1 , 1 
SU'^ » Sum ♦ SuARCA( kCw, CLM) ■> y S K 3X 3 ( R T I, , CL M . M S K ,\ C X ) 
10 CCMINUe 

2C CCMlMJc 

IF (A&S(SUM> .LT. kAXVAD CHTC 5C 
MAXVAL ’ ASSIS'JMJ 

PSKNC^' » *'SK.\CX 

50 CCNTINUF 

C 

c 

IF (MSKNC'*' ..Nf. 1) GOTC (.C 
ThSSHl = S0ARLA(2,3) - CENTER 
THRSH2 s S3ARcA(2,l) - CENTER 
GOTO 90 
C 

60 IF (MSKNUM .NE. 2) GQTC 7C 

THPSHl = S0AREA(l,3) - CENTER 
THR3H2 - SCAREAO.l) - CENTER 
GOTC 90 
C 

70 IF (MSKNUM .NE. 3) GOTC 80 

THPSHl = S6AREA(l,2) - CENTER 
THRSH2 = SbAPEA(3.2) - CENTER 
GOTC 90 
C 

80 THRSHl = SDARrAd.l) - CENTER 
THRSH2 = SBAREA(3.3) - CENTER 
C 
C 

90 IF (AflS(THRShl) .GT. A6S(THRSH2)) MSKNUM = MSKNUM ♦ A 
C 

RETURN 

ENC 

EOF.. 


tf' { 


n n 
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C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

t 

c 

L 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


LMVFRSITV CF KA\‘J 

V, 1 FLFCCMPL'JI CAT ICNS AnT INFC.IPATICN SCI(NL‘‘3 LA? 

PunCRAM 

SLiir ; i.n St: fu teps 

RE^ 

. n : 

PRCCR AM 

NAPE :FCSTAt 

AIIThC'TrJ. SCCTT G'.RON'R C A T C : C Z / 1 <> / ’ 2 

PLR POSE 

THIS SJ 

.'’PiJllTINE CCPPLTLS T f- 

r LOCAL STATISTICS ICR 

A ECCE 

UINOCW AFTER 

AOPL Y INC 

A PASK WHICH POST 

PPCVICLI 

5LY HAVE SEEN NORPAL IZ 

CC. 




PAR APcTC R 

OEFINITICN 

NAPE 

\ TYPE \ CLASS\ RANGE 

\ DESCRIPTION 

W INCCU 

\R 

\R 

\ 

\ LOCAL AREA NINCCV. 

MSKVtNC 

\R 

\R 

\ 

\EOGC- TEMPLATES 

WNCSIZ 

\ 1 

\R 

\ 

\SIZE IF WINCCU 

PSK NUP 

\ I 

\R 

\ 
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Abstract 

Classification of synthetic aperture radar (SAR) imaqes has important 
applications in geology, agriculture and the military. A statistical model 
for SAR imaqes is reviewed and a naximum likelihood classification algorithm 
developed for the classification of agricultural fields based on the model. 

It is first assumed that the target feature information is known a priori . 

The performance of the algorithm is then evaluated in terms of the probability 
of incorrect classification. A technique is also presented to extract the 
needed feature information from a SAR image; then both the feature extraction 
and the maximum likelihood classification algorithms are tested on a SBASAT-A 
SAR image. 
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1 . 0 Tntr<~i1m-ti on 


Synthetic Apertun? Rnrtar (SAP) imAninn systoms )wivf been iisei to obtain 
imaqes of the Earth's surface. This pafnT finals with the prohl* n of antnnatic 
classification of SAP imaaery. 

Accordina to t)ie definition in (1), classification of an object in the 
imaae is the decision reqardinq the catnoory to which the object helonqs. For 
example, if the imaqe is of an agricultural area contalnino a number of fields 
and there are three cateqories, or classes of fields, such as earn, wheat, and 
soybeans, classification is the process of assigning a cateqor'" to each field. 

Each object (tarqet) class is characterized by a unique set of measurable 
properties, or features. Itie classification is done by assiqninq the tarqet 
to a class based on how closely the observed set of target feature:, .matches 
the set of features for that particular class. The assignment of an object to 
one of several classes can be done on the basis of a maximum liltelihood 
criteria, that is, a classification decision is made such that the probability 
of incorrect classification is minimized. The problem of maximum li)telihood 
classification has been presented in general form in [91, and solved for the 
specific case where the observed feature vector is distributed js multivariate 
normal random variable. In our case (SAR imaqes), the classification is based 
on one observed feature, gray level, which is not normally distributed. The 
major problem with processing SAR images is that coherent nature of the 
microwave illumination gives rise to a phenomenon called speclcle. Specltle 
noise seriously degrades the quality of an image [12]. It is signal- 
dependent. Because of the dependency of the noise on the signal, processing 
techniques designed for additive Gaussian noise fail. Thus other statistical 

models are needed, and one which has been successfully used in the past 

[15,16,20] will be introcuced in the next section. 

As has been mentioned above, a tarqet in a SAR imaqe can be described by a 

single feature - gray level mean. This is especially true for agricultural 
targets, where textural features do not offer much discrimination. Here, a 
maximum lllcelihood classifier based on gray level mean is designed. This 
classifier assigns every pixel of an image to one of several tarqet classes. 

It is originally assumed that the gray level means are Icnown a priori for all 
target classes. A method is then presented to estimate those values 
automatically. To evaluate the performance of the classifier, the 



probability of cla*5«tf Irntion rrror is dorlv^d, an«l thf? classifier performance 
is tested usinq radar imaqe sinuilation and SFASAT-A SAR imaqes. 

2 , 0 Statistical Model for Intenaities of Pixels in Synthetic Aperture 
R adar Images 

It has been confirmed [iC] that heuristic imaqe processing techniques, 
such as gradient edge detection algorithms do not work well for SAR images. 
Thus a statistical model is needed, one which would incorporate the physics of 
the SAR image formation process. This model is the basis for image processing 
techniques developed later. 

The main ch >• acteristic of SAR images is that the noise variance of 
the pixel Intensity X is proportional to » the square of the mean of this 
pixel. The true intensity is proportional to power return from the pixel 
n the absence of noise. The ramifications of proportionality of mean and 
variance is that targets with higher intensity have higher noise variance, and 
one can conclude that the noise power is the function of the signal power. 

This precludes the use of an additive white Gaussian noise (AWGN) model, where 
signal power and noise power are assumed to be independent. 

Synthetic aperture radar images belong to the category of speckled 
images. The statistical characteristics of speckle noise have been developed 
in ( 13 ] and applied to synthetic aperture radar in ( 14 , 15 ). The complete 
noise model was developed in ( 16 , 201 . According to this model, let R be a 
random variable representing observed pixel inte;isity from the SAR image, it 
has been shown that the probability density function (pdf) of R is exponential 
with the mean of p where p is the true intensity of the target to which this 
pixel belongs. It has been shown in (15) that pixel intensities of 
neighboring pixels can be assumed to be uncorrelated. Here we also assume 
that they are independent. In raost SAR imaging systems more than one 
independent sample of detected power is obtained for each pixel. Tbe number 
of independent samples, a is often referred to as "number of looks per 
pixel". Now pixel intensity X is obtained by averaging these independent 
samples: 

1 ° 

- a "^i , 

The probability density function of X can thus be expressed by the formula for 
OAmma distribution: 


2 


(1.3) 


f^(x) 


r(a) 


where, 

B .■ U/a 

r(a) ■ (a-1)! for a an inteqer. 

Mean and variance of pixel intensity X can be calculated usinq the 
properties of qamma distribution. They are equal to u and v /a, 
re*^-Ctively . The ratio of siqnal power to noise power (siqnal-to-noise 
ratio) [4] can be calculated as 

S _ E^[Xl 

N var(x) ^ 


Thus it has been shown that noise power is proportional to siqnal power with 
coefficient of proportionality of 1/a. 

This relates directly to the observation about dependency of noise on 
siqnal, made at the beqinninq of this chapter. Statistical noise model for 
SAR imagery can thus be best characcerized as "multiplicative noise model", 
where siqnal and noise power are proportional to each other: 

X - S • N (1.7) 

where X is defined as before, S is noise-free signal, and N is random 
variable, representing multiplicative noise. One should note that signal S is 
a continuous random variable and its probability density function is 
continuous. 

There are restrictions to the model, worth mentioning: 

(1) The statistical model for pixel intensities introduced above is only 

valid for a type of region that can be defined as homogeneous. A homogeneous 
region is an area where all pixels belong to the same target, best 
characterized by a single feature: mean y. If more than one target is 

present in the region the distribution of pixel intensities does not follow 
gamma distribution, and more sophisticated model is needed [16]. Such regions 
are called edge regions, and since the model is not valid for them, they 
should be excluded in any analysis. 

(2) There are several conditions to be imposed on the target for 
multiplicative noise model to be valid for SAR imagery. These conditions were 
outlii ed in [13]. One of them is that the surface roughness of the target 
should be large, compared with wavelength of radar signal. There is a variety 


3 


of tarqet classes that would meet such conditions. An example is a tarqet 
that contains an aqricultural area, such as corn field. 

In summary, the statistical model that incorporates properties of SAR 
imaqe formation process was introduced, and its limitations mentioned. This 
model will he used throughout this paper. 


3.0 Maxlnuw t.ikplihotyi Classification of P«»nion8 in tht: SAH Iwanos 


This section presents a solution to the problem of the classification of 
pixels in a SAP inane, based on qrav level. Also, an expression for 
probability of incorrect classification is derived, and a test to exclude edqe 
reqions is presented. 

In this section it is assumed that the imaqe contains M tarqets. The true 
tarqet means are known a priori . (An alqorithm to extract tarqet means from 
the SAP imaqe will be presented later). A slidinq KxK processinq window is 
applied to the imaae, and sample mean is calculated for each window 
position. *lhen, based on the outcome of the test to be described, the center 
pixel of the window is assiqned the value of the tarqet mean, to which the 
region, defined by the window, is most likely to belonq. 


A. Maximum Likelihood Approach to Classification 


For a given window position, let X^, ...» Xjj (N»KxK) be the gray 

level intensities of the pixels within a window. According to our model, we 
can assume that , . . . , X^ are independent and uncorrelated. 

1 ” 

Let X ■ — X, be the observed mean, and X. , X. , . . . , X. be true means of 
” iTi ^ 

all M targets in the imaqe, known a priori , as previously assumed, it is 
assumed that all the pix»»lc in the window oclonq to the same homogeneous 
region and that the region is characterized b/ the mean . Peqions that 
contain a mixture of tarqets (edge reqions) will be excluded by a test to be 
described later. Now consider the following set of hypotheses: 




X. (region belongs to Target 1) 


p ■ X, (region belongs. to Target 2) 


(3.1 ) 


u » Xj^ (region belongs to Target M) 

Tbe likelihood function under hypothesis is defined as joint probabi- 
lity density function of all samples under the hypothesis 
H^, 1 < i < M, multiplied by a priori probability of hypothesis being 
true, P(H^). Using the SAP imaqe model previously described, the likelihood 
function can he written as 


L(X. 


X,, X^, 


•*' ’'n' 


P( 


( M a-l\ -oNX/X 
n X le 


(T(a)) 



(3.2) 


5 


region 


. The test is carried out this wavt 



less than — 
2 


, the region 


is an edoe region, '.tireshold of — has been chosen heuristically . To identify 

2 

edge region, the center pixel of the window is assigned zero intensity, and no 
further testing described in previous section is needed. This test is a quick 
way to check whether the region defined by a processing window at a given 
position is homogeneous. However, it is not necessarily the optimal test. It 
may miss some of the edge regiont . A better test (Likelihood Ratio Test) has 
been derived in (16). Likelihood Ratio test would define the edge regions 
better at the expense of computational efficiency. 


C, Performance Analysis of Maximum Likelihood Classification Algorithm 

Section 3. A considered the problem of assigning center pixel of processing 
window in a SAR image to one of M targets. Tliere is probability that the 
pixel is assigned to the wrong target. Hie purpose of this section is to 
derive an expression for this error probability. In [231, a performance 
analysis criteria were derived for two targets in the image and with Gaussian 
noise model assumed. Hie following analysis uses the multiplicative noise 
model and M targets in the image for derivation. Let 


Y. ■ — — + in X . , 1 < i < M, be our decision criterion. Let P be the 
i Li e, 

probability of assigning the pixel to a target other than target 1, when 

hypothesis is true (that is, pixel belongs to Target 1). 


Identically, p^ 


1 < i < M is the probability of assigning the pixel to 


target other than correct target with mean X,,. P^ can also be expressed in 

IrfX 6 £ 

terms of probability of correct decision: 


P. - 1 - P„ (3.6) 

*i *=i 



hypothesis is true. Since all targets occur with equal probability, the 
average probability of misclassif ication can be calculated as: 


M 


-E 


P /M 
®i 


(3.7) 


Thus the task of performance analysis boils down to finding P , 1 < i < M. 

^ i 
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Consider the followinq set of statistics: 


Y, - + tn X. 


1 < i < M 


(3.8) 


Let the true tarqet means be arranqed in ascendinq order: < Xj^j < ... < 

Xlm* probability of correctly assiqninq the pixel to Tarqet 1, when 

hypothesis is true can be calculated as; 

Pc, “ P(<Y, < Yj, Y, < Yj Y, < Y^)|hJ (3.9) 

where the comma in probability expression is defined as loqical "and”. Next, 
calculate the difference of random variables Y^ and Yj where 1 < i < M and 1 < 
j < M: 


’'i - ’) ■ 4: * ‘"Xu - 


*LJ - hi . hi 
^ Tx x“ 

Lj Li Lj 


(3.10) 


Hie inequality Yj^ - Yj < 0 can thus be expressed as: 


5 < (tn ^ 

Li Lj • Li 


(3.11) 


By inspection, one can conclude that 


’^Lj^^Li . ’'Lj 

Let . be equal to — — In — — , 

' ’ ’'l j ■’'li ’'li 

Zj^^j > 0 and Zj £ ■ 55^ j • Rewriting the probability of correctly assiqninq 
pixel to target 1 , we obtain: 


■■■ P-, -'P((X < Z, X <' 2 / 3 ', ..., X < z, „)|hJ 


(3.12) 


This joint probability can be expressed as the product of the following 
conditional probabilities: 


**C1 - " ^,3>^ " ’'l,2’l”J 

P(X < z, J^|x < Z, 2' •••' ^ < ^1,M^I”l^ 


(3.1 3) 


But Xj^, < Xj ^2 < ... < therefore Z, ^2 < ^ 1,3 < •••• Z, < Z, so 

all conditional probabilities in the expression above are equal to unity i.e.. 


p( (5? < 3 IS < 


- 1 


•<'* ‘ ' ^,2' * < '-1,2 


f •••pX^Z 


l,M-.> "i> ■ ' 


and 


(3.14) 


By analogy the probability of correct classification for the largest a priori 
n.ean is: 

P(!i > z, P((8 > Zj,„|!I > Z,,„>|H„) 

P(S > X > •••' ® ^ 


Again, all conditional probabilities are equal to unity, and 
PcM - P(X > Z,^„|h„) 


(3.15) 


For any 1 < i < M, the expression becomes: 

Pci ■ "* ^i,l' ^ ^ ^i,2' •••' * ^ ^i,i-1' ^ ^ ^i,i+1' *•*' ^ ^ ^IpmI^U 

Or: 

Pci - P<5c > Zi,i_,|Hi) P(X < > z. i_,)|Hi) 

It can be easily shown that the only tarqets that the pixel can he incorrectly 

assigned to when it belongs to target i, are tarqets with true means X. and 

^i-1 

X. , directly below and above true mean X, 2 . The probability of correct 

i+1 _ 

decision can be expressed as: 

Pci - Pt’' > Zi,i-ll»i) P(5' < Zi,i+ll5^ > Zi,i-l)l»i) 

■ P*^'i,i-1 < ^ < 2l,i+il”i^* (3.16) 

In summary, probabilities of error Pg^ ■ ^"Pci' calculated as: 

P„ -P(S > Z,,2lH,l 




(3.17) 


PeM - P(X < 

Since the processing window size is large, one can apply central limit 
theorem to fi:id the probability density function 55 under any hypothesis 1 < 

i < M. Expected value of X under hypothesis is: 

N N 


FixiHlI --5 E/lXil - S 2 ■ *Li 
1-1 i-1 


(3.18) 


9 


And variance of X under the aamo hypothesia iai 


Var (S|Hj^) 


1 




Var(X^) 



No 




(3.10) 


Thus, according to the central limit theorem, it can he assumed that Jc has an 

approximate normal distribution under arbitrary hypothesis Hj^, with mean of 

2 

Xj^£, and variance of Xj^£ /oN. Let be equal tot 
X - X . 

R, - ( — 5 ) /on (3.20) 

Li. •, . 

Under hypothesis R^ is approximately normal with zero mean and unit 
variance. Now we can express error probabilities in terms of R^t 

P - p(r > (-1^1 /^|h ) 

.1 


- P(- 


" --H) /^< Pi < /Sn)|h ) 

Li Li 


(3.21 ) 


P - p(R < 

"m ” *LM 


^i.j - ’'Li 

Tbe threshold ( — ) /aN is equal to 


Li 


'i.j 

where T 


_ (tn t.l - Tj_j) 


i.l 


— is the relative target contrast of Target i with respect to 
Lj 

Target j, 1< I<M, can also be expressed in decibels. Now 

error probabilities can be expressed in terms of Q function, where the p 
function is defined as: 


1 

rs -i 


(3.22) 


The error probabilities can then be calculated as follows: 


i) /^)) 


p 


el 


of( 


ln(T, ,) + 1 - T 

T - 1 

1,2 


* T. . . + 1 - T. . . InT. . , 1 “ T. . 

%i- '-«<( 4 ^ - «(( 


i,i+1 - 1 




(3.23) 
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1-P(( 




1 - T, 


M,M-1 


T - 1 

M,M-1 
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D. Conclu8ion3 

Inspecting the expressions, derived in Section C one can malce some 
important conclusions: 

The maximum likelihood classifier can incorrectly assign the center pixel 
of processing window of size N only to targets with true means either above or 
below the mean of the target to which the reoion belongs. The error probabi- 
lity depends on two relative target contrasts Tj^ i-l* "^i i+i • experiments 

to be discussed in a later section are consistent with these conclusions. 

Hie performance of maximum likelihood classifier depends on the product of 
number of independent, samples averaged by the imaging system {.lumber of looks 
per pixel), a, and number of pixels in the processing window, N. Since the 
error does not depend on a and N individually, but rather on their product, 
there is a tradeoff between a and N. For instance, performance of maximum 
likelihood classifier on an image with 4 looks per pixel (a>=4) and 25 pixels 
in the processing window (N>25) is equivalent to performance on the image with 
10 looks per pixel and 10 pixels in the processing window. Thus we can call 
aN total number of looks or total number of samples . 

Hie classification decision, derived in this section will yield the best 
error performance for a type of SAR image consistent with our statistical 
model. It has been shown in [9] that maximum likelihood classifier will 
minimize error probabilities. Hierefore, this technique results in optimum 
performance for the given statistical model. 
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4.0 Targpt M«an Fxtraction T>achniqu«» 


The maximum likelihood claesif icatior. alq^^rithm, derived previously, 
assumes that true tarqet means are known a priori . In many cases these moans 
may not be available. 'Itiis section deals with the problem of estimatinq 
target means, when the estimation process is complete, the resultant tarqet 
means can be supplied to maximum likelihood classification algorithm. A test 
to compensate for imperfect extraction is also described here. 

According to our statistical model, the ;^ixel intensity of a homogeneous 
area is. a. . gamma distributed random variable. Therefore the probability 
density function* (pdf ) "of the whole" imaoe containing a number of targets can 
be characterized as pdf of a mixture of a number of gamma distributed random 
variables. The problem of estimating target means is equivalent to estimating 
the parameters of each mixing gamma pdf. Such estimation is possible if the 
mixture is identifiable (for definition of identif lability see [28]). There 
are methods to estimate parameters in identifiable mixtures (9). However, the 
mixture of gamma variates is identif idOle only if the random variable 
representing tarqet means takes on only discrete values; otherwise the mixture 
is generally not identifiable [291. Unfortunately, in our case the random 
variable representing ♦■arqet means has continuous pdf, and the mixture 
parameters cannot be estimated. Even if we assume that pixel intensities have 
approximately normal distribution, the mixture is still unidentifiable, 
because mixtures of normal pdf's are identifiable only if all normal variates 
used in the mixture have equal variances [28]. In our case the variances and 
the means for each area will be unequal. Therefore a different approach 
towards estimation of tarqet means is needed. One of the approaches 
considered here is based on jelecting the area that is most likely to be 
homogeneous and estimating target mean from such area. 

A. Automatic Extraction of Target Means 

The purpose of this section is to describe an automatic (unsupervised) 
target mean extraction procedure. The basic idea of the procedure is the 
following: 

(a) The homogeneous areas of the image are identified, e.g. the ones that 
belong to the same tarqet. 

(b) These areas are combined into groups on the basis of likelihood of 
belonging to the same tarqet. Therefore each tarqet is identified by a 
respective group. 


(c) All pixel IntensitieB within e/ich qroup <ir« avoraqi'd to obtain 
eRtimates for tarnet meana. 

In RumiTO y, pixels from homogeneous areas defined by the algorithm serve 
as a basis for estimating target means. Therefore the test to identify such 
regions has to be very stringent. The probability of selectinn an edge area 
(e.q. area containing miltiple tatgets) as homogeneous must be small. 
Conversely, the probability of missing (rejecting) a homogeneous region is not 
important. For example, if the target contains SOO local areas and only one 
of SOO is selected, this is still enough to provide a good estimate for the 
target mean." However, if a non-homogeneous area is selected as homogeneous, 
the estimate based on such area is wrong. Therefore, trying to minimize the 
probability of selecting edge areas as homogeneous is an important 
consideration. 

The complete flow chart of the procedure is shown in Figure 1. The 
following is the description of flow chart block by block. 

(1) The first step is to apply a .tliding processing window (typical size 
13x13) to the image. Then, for a given window position, the mean and variance 
of pixel Intensities within the window is calculated. 

(2) The variance test is a quick way to check whether the neighborhood 
(defined by the processing window) is homogeneous. It is basically a moment 
matching technique that tells how closely the observed variance of the 
neighborhood conforms to maximum likelihood estimate of that variance. 

The maximum likelihood estimate of the variance is predicted from the mean and 
is equal to: 

^2 X . . 

0 - - (4.1 ) 

a 

a is a number of looks per pixel and is always known for a particular SAR 

I 2 ^2 1 

imaging system [16] . If |S - o | < P where P is a threshold, the region may 
be homogeneous and chi ;uare test is invoked to impose tighter constraints. 

If this inequality fails, the region is not homogeneous and no further testing 
is needed. The sliding window is then moved to another neighborhood. 

(2) The next test for homogeneity uses a chi-square goodness -of -fit 
test. This test is a more stringent check of homogeneity. It provides a 
qualitative measure of how closely the gray level distribution measured from 
the neighborhood fits the gamma distribution, predicted by the model. 
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Figure I. Flowchart of the target mean extraction algorithm 









Sp»*cifically, l«t , Xj Xjj l>e pixel intensities from the neighbor- 
hood. Aqaln, accordlnq to the model, we cr.n assume that random variables X^ , 
.... Xjj are independent and uncorrelated. A hypothesis test of the followinq 
form is then used: 


H t f(X. ) - 
o 1 

f(Xj ) * 




1(a) 

'b-"v 

T7TT 




1 < i < N 
1 < i < N 


(4.2) 


where flX^) is the probability density function of a given pixel X^; a is 
number of loo)cs per pixel and 0 is u/a, as defined before. Since true mean u 
is unlcnown, it can be replaced by its maximum likelihood estimate. 


X 


1 

^ i-1 ^ 


(4.3) 


This is equivalent to the following hypothesis test: 

HqI neighborhood is homogeneous 
t neighborhood is not homogeneous 

Dividing the range of pixel intensities into K intervals, the probability of 
rando'< variable X^ being within given interval (a,bj is; 


P(a < X^ < b) 



a-lg-Oe-Xi/B 

TTaT 


d - F(b) 
*i 


F(a) 


(4.4) 


Where >'(/) is the distribution function of the pixel intensities. Thus 
expected number of samples from a neiohborhood to fall within an interval 
(a,b) is N(F(b) - F(a)). The observed number of samples can be calculated 
from pixel intensity histogram for a given neighborhood. Define expected 
number of samples to fall within ith interval as Nexp^, and actual number of 
samples as Nobs^, respectively. Then define test statistic, x*' as follows: 


K (Nobs, - Nexp, )‘ 

Z 


i-1 


Nexp, 


(4.5) 


where K is the number of intervals 


X has an approximate chi-squtire distribution with (K-l ) denrees of 
freedom, x defined as our measure of "deqree of homooenei ty" . Ttie oreater 
X is, the more non-homoijenoous the neiqhborhood is likely to be. Thus we 
define homooeneous neiqhborhoods by the followinq test. We accept if x < 
X,p, otherwise we reject H^. Threshold x>j, determined by the siqnificance 
level of the test, and the typical values for x^* between 1 and 7. If we 
accept the hypothesis that the neiqhborhood is not homoqeneous, we qo on the 
next neiqhborhood. 

(4) Havinq established that the neiqhborhood is homoqeneous, and if the 
neiqhborhood is the first homoqeneous neiqhborhood in the imaqe, an initial 
estimate is created for the true mean of the first tarqet. iTiis estimate is 
obtained by takinq the average of all pixels in the neiqhborhood. Also, the 
number of pixels that were averaqed to obtain this initial estimate is 
recorded, alonq with processinq window coordinates. 

(5) If the homoqeneous neiqhborhood is not the first one in the imaqe, 

then there already exists some table of initial estimates. There are two 
possible outcoi.-.*s: (a) Tbe neiqhborhood belonqs to a tarqet, for which 

initial mean estimate has already been created; or (b) The neiqhborhood 
belonqs to a tarqet for which no initial estimate has been created. To check 
which of the two possible outcomes is true, the followinq statistic is 
calculated: 


Y. -ii 

i X 


+ in X, 


Li 


Li 


m"in (— + in X. 


• • • f ^ +£nX_.) 


where X is the local neiqhborhood mean, Xj^^ , X^2‘ 


X^j are initial mean 


estimates for the first j tarqets, and 1 < i < j. The statistic is nothinq 
more than the maximum likelihood criterion, developed in section 3. 

Now the two possible outcomes can be reformulated as (a) Neiqhborhood 
belonqs to the tarqet with mean Xj^^; and (b) Neiqhborhood belonqs to a new 


tarqet with mean X. 


and initial estimate is : X^ 


X. 


aiiia xiixuxaj. cnuxiuauc; xa> : 

Under the hypothesis that the neiqhborhood belonqs to tarqet with mean 
X, , Y; is approximately normallv distributed with mean 1 + InX, and standard 

bi X 

deviation of 1/oN. If the observation Y^ is within three standard deviations 
of its predicted mean 1 + fn Xj^_ , outcome (a) is likely. Otherwise, outcome 
(b) is likely. The measure of three standard deviations has been chosen 


heuristically. 
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(6) If in not within 3 standnrd doviationn (3//^) of its moan, a new 

initial tarq^'t pstinate - 5c is created. Number of pixels in the 

neiqhhorhood is also recorded, alonq witu processinc window coordinates, that 
define the current neiqhhorhood. Then the processinq window is moved to the 
next neiqhhorhood. 

(7) If Y^ is within 3 standard deviations of its mean, then the neiqhbor- 
hood is most likely to belonq to the tarnet with initial mean estimate of 

1 < i < j. If that's the case, the tarqet mean estimate can be updated by 
takinq a weiqhted averaqe of already existinq estimate X, and current 
neiqhhorhood averaqe S. 

It is important to assure that this weighted averaqinq required to obtain 
an estimate for a tarqet mean is done over non-overlapping neiqhborhoods, e.q. 
the neighborhoods that don't have common pixels. Since each neighborhood is 
defined by processing window of fixed size at a given position, and window 
coordinates have been recorded, it's easy to check whether the current 
neighborhood overlaps with all other neiqhborhoods which were used to obtain 
an estimate for a given tarqet mean X. . If the current neighborhood does 
overlap, its oray level averaqe cannot be used to update given tarqet mean 
estimate X, . 

(8) If the neighborhood does not overlap with all others that belong to 
the same target, and pixels from which were used to obtain an estimate for 
tarqet mean, then the initial estimate Xj^^ is updated by weighted averaging: 


X. ■ 

Li new 


X. • , j • H , . + XN 

Li old old 

H , . + N 
old 


(• 1 . 6 ) 


”new “ ”old 


(4.7) 


where la number of pixels that were used to obtain initial estimate, N is 

the number of pixels in the neighborhood. X^ and are initial 

estimates after updating for a tarqet mean, and number of pixels used to 
obtain that estimate, respectively. The current window coordinates are also 
recorded. 

This procedure is a recursive procedure, the result of which is a vector 
of target mean estimates. Its performance will be evaluated in section 5. 
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A TVst to Conponn.it.''* for Imnorfoct. Tnrnnt fiortn Fxtr.ictlon 


The procoduro doncribod in the previous section leads to tarqet mean 
estimates to be used as a inj>ut to maximum likelihood classifier. Fut these 
estimates may be imperfect. For instance, some targets may be missed, and 
others mav be estimated incorrectly. Hiis will increase the probability of 
classification error. However, there is a test to partially compensate for 
imperfect target mean estimates at the second (classification) stage. 

As before, let be our decision criterion; 



+ in X ■ 
Xj . 
1 


Suppose that the neighborhood gets assigned to arbitrary target i with mean 

X. . This test is similar to merging criterion described in the previous 

section. Under the hypothesis being true, Y^ is approximately normal with 

mean of 1 + in X. and variance of 1/aN. If the value of statistic 'i- is far 
‘'"i 

enough from its expected value 1 + in X. ("far enough” is a least three 
standard deviations: 3//aN ), then the center pixel of the neighborhood 
defined by processing window is assigned zero intensity. This test implies 
that even though the neighborhood belongs to a given target on the basis of 
maximum likelihood, its decision criterion's value is far apart, from the 
expected value of that decision criterion. This may occur because the 
neighborhood belongs to a target, the true mean of which is missing, or has 
been estimated incorrectly. Therefore we don't really know which target the 
neighborhood belongs to and thus assxgn zero intensity to its center pixel to 
show that target couldn't be determined. Similar testing is done in the area 
of digital communications. When at the receiver it is not possible to 
determine whether the received bit of information was zero or one because of 
high ambiguity, the output bit is assigned "don't care" condition. Practice 
of assigning "don't care" condition to cases when it is not possible to 
determine correct status is often called "erasures", and is described in the 
communication literature (26). The test developed here will be referred to as 
"erasure test" from here on. 
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5.0 Perfornanco 1 nation 


A . Verification of P«»rform.incn of thn Maxiwiin I.iknlihood Cl.aasifior 

A computer program based on the maximum likelihood classification alqor- 
ithm described in Section 3 was written and tested on SAR imaqe simulations 
and real SEASAT-A SAR imaqe 1271. Simulated images were of size 250 x 250 
pixels and contained 14 tarqets. This section discusses simulation results. 
The advantaqe of usinq simulated imaqes to test the performance of the maximum 
likelihood classifier is that the location of each tarqet and tarqet means are 
known a priori . Tbe description of simulation procedure is available in 
[21]. The basic idea of SAR imaqe simulation is that a noiseless image (where 
all pixel values are equal to values of their respective target means) is 
multiplied by a random variable, which incorporates desired speckle, e.g, 
number of looks per pixel, a. As a result, t’le noisy image arises, with the 
noise statistics fittinq the multiplicative noise model described 
previously. Such simulation is a good representation of a real SAR imaqe. 

The noiseless imaqe that serves as an input to simulation is called power 
map . The power map in Figure 2 was used for this experiment (the numbers in 
squares represent true tarqet mean); 
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Figure 2. Power map for SAR image simulation 


This image contains 14 targets 50x50 pixels each. The target means were 
selected such that a wide range of relative tarqet contrasts would be present 
(from 0.38 dB to 3 dB). This is done because the probability of classifica- 
tion error depends on relative tarqet contrasts (see Section 3.0. Also some 
of the same tarqets are separated and located in the different parts of the 
imaqe to illustrate the fact that tarqets with the same mean can be 
disjoint. Five simulations we^e made, based on this power map. Imaqes with a 
■ 2, 4, 6, 8, and 10 were created. (See Figure 3). Then they were processed 
with 9x9 sliding window bv the maximum likelihood classification algorithm. 


19 




























■I 4/«i^ page 

^0 WH,n;^roGwp„ 



Figure 3. 


Results of Maxlmuffl Likelihood Cl&ssl f Icatlon Processing (a) a 
(b) a - 4 (c) a - 6 (d) a - ft (e) a - 10 
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The results of processinq are available in Fiqure 4. JYom inspectinq the 
results, one can see that classification alqorithm's performance improves as 
number of looks per pixel qoes up, for a fixed window size. The black (zero 
intensity) points represent the edqe reqions selected by the test, described 
in Section 3.B. This edqe region exclusion test defined t ^st (but not all) 
local areas containing edqes. it also defined as edge areas some areas that 
are homogeneous; but that was not a major problem, because such areas only 
constitute a small proportion of the image. The classification is more 
adversely affected by missing edge areas than defining false edges. For any 
given target, the misclassif ied pixels were assigned to a target with mean 
either above or below the mean of the given target, as expected. Also, the 
highest number of misclassif ications occurred for a target with the lowest 
relative target contrasts T^ or Tj^ i + i» defined in Section 3.C. All of 
thesv* observations agree closely with the theory. ^ 

Figure 5 shows the plot of average probability of error -jUj ^ P^^ and 

maximum probability of error versus total number of samples, oN, for the 
simulations described above. These curves enable the operator to choose the 
right processinq window size /n x /TT to achieve the desired error performance 
for a given number of looks per pixel, a. 


THEORETICAL PROBABIUTY 
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Figure 5 . Predicted Probability of Classification 
Error for Simulations From Figure 2. 
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FvrtlviationofTiiroot M<»an Fxtraction T'nchniauo 


Ttie tarq«t m«an extraction alqorithm was applied to the simulations 
described previously. Since the true means of all taroets in these 
simulations are known a priori , it is possible to compare these means with the 
means, extracted by the alqorithm for different numbers of looks per pixel. 

The threshold on a chi-square test, x<i.» chosen to be 1.0, and the ranqe of 
intensities was divided into 10 intervals, each havinq equal expected 
counts. Ibe processinq window was chosen to be 13x13. Table 1 comnares 
results of the extraction procedure described above with known true means t 


True Mean 

Extracted Means 

a - 2 

a ■ 4 

a - 6 

a ■ 8 

a - 10 

5.0 

8.11 

4.93 

• “4.95 

^ • -5.20 

5.04 

10.0 

9.90 

18.27 

10.11 

- 

10.1 7 

20.0 

22.27 

20.85 

26.82 

27.71 

- 

30.0 

30.20 

29.88 

30.87 

30.69 

28.75 

40.0 

38.20 

40.01 

39.59 

40.73 

41 .22 

60.0 

63.21. 

60.22 

- 

57.71 

57.69 

80.0 

- 

80.97 

84.12 

80.10 

81 .45 

100.0 

100.34 

98.75 

101 .04 

- 

102.71 

120.0 

1 27.36 

115.50 

118.39 

123.54 

118.42 

140.0 

- 

138.87 

138.71 .. 

135.76 

140.37 

160.0 

. 

— 156.09 _ . 

. 154.08 

- 

159. <.7 

180.0 

170.43 

- 

185.84 

177.28 

- 

220.0 

212.82 

- 

- 

270.57 

- 

240.0 

- 

- 

- 

247.17 

246.89 


Table 1 


One can observe that extraction technique has been able to estimate most 
tarqet means correctly. However, there are some imperfections, like a few 
missinq tarqet mean estimates, and a few incorrect estimates. Effects of 
these imperfections can be partially compensated for by the test described in 
Section 4.H. Ttie results of applyinq maximum likelihood classifier with 
erasure test to the simulations utilizinq target means extracted by previously 
described alqorithm are qiven in Figure 6. The dark areas represent the 


24 



ORIGINAL PAGE 

black and white photograph 


I 

I 

I 

I 

I 

I 

r 

r 

r. 

r 

£ 

£ 

£ 

r 

1 . 

t 



FlRure A. Results of Processing Simulations with Target Means Estimated by 
Extraction Algorithm ("Erasure test" was used to partially 
compensate for imperfect extraction) (a) a - 2 (b) a - 4 

(c) a ■ 6 (d) 0*8 (e) a ■ 10 
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regions, where classifier failed to make a decision. One can illustrate the 
effectiveness of erasure test by the following examplei 

In the case of a ■ 2 simulation, the extraction routine missed a true 
target mean of S and extracted a false target mean of 8. Itie true target mean 
of 10 was found correctly and its estimate was 9.9. In the original image, 
targets with means of 5 and 10 are 50x50 squares next to one another in the 
bottom right-hand corner. Since true target mean of 5 was missing, the whole 
area which was supposed to belong to that target was assigned zero inten- 
sity. If no erasure testing was done, this whole area would have been 
assigned the value of false target mean of 8. The next square has a true 
target mean of 10, and most of its pixels were classified likewise. Only very 
few pixels were assigned the value of false label 8, most of them at the edge 
of areas with true target means of 5 and 10. thus, a number of classifica- 
tions to the false label was significantly reduced. 

Of course. If mean of 8 was not extracted, there wouldn't have been any 
misclassif ications at all. Therefore, it always degrades the performance more 
to extract a false target mean then to miss a true one. Yet the test, 
described in 4.B provides some degree of protection in both cases. 

C. Application of Complete Algorithm to Real SEASAT Image 

The SEASAT image (27] used in this experiment wac 512x512 image of 
agricultural scene. This image is shown in Figure 7. The number of looks per 
pixel, a, was approximately 3. The image was processed by maximum likelihood 
classifier with 9x9 processing window size with edge test and compensation 
test included. Target means were extracted manually, and automatically. 
Resulting images are given in Figure 8 and 9, respectively. Table 2 lists 
target means, extracted in both cases. Note that neither of the extraction 
techniques is perfect. 
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Figure 9. SEASAT-A-SAR Image 
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Figure 10. Result of Processing SEASAT-A-SAR Image (Target means extracted 
manually) 
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•11. Result of Processing SEASAT-A-SAR Imsge (Tsrget mesns estirasted 
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Real SFASAT Image Target Means 


Fatima ted Manually 
15.3 
25.0 

41.9 

48.9 
59.8 

107.8 


Fstimated Automatically 
18.05 
27.50 
40.96 
51 .28 
62.69 
72.57 


Table 2 

The performance in both cases is quite reasonable, especially compared with 
results of processing the same image using other techniques such as different 
types of qradient edge detectors (161. A major problem is high misclassif ica- 
tion rate at boundaries of targets (edges). This occurs because the statisti- 
cal model doesn't apply to edge regions, and the mean^/variance test described 
in Section 3.8 fails to select all the edge regions. This problem would be 
reduced if more sophisticated edge test was used, such as li)celihood ratio 
test, available in [16]. 
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6.0 Conclusions 

To aumnarize a n.Dcimun likelihoo*! claasif ication algorithm was dev«lop»»fl 
for SAR inaaos, and its theoretical performance was evaluated. Also, 
automatic extraction algorithm was developed to estimate target mean levels, 
and a test to partially compensate for imperfect extraction was introduced. 
The algorithms were tested, using simulations and SEASAT-A SAR imagery. 

The classification scheme developed here is the best classification 
possible based on the given statistical model. This statement can be made 
because the maximum likelihood approach minimized the probability of 
classification error (9]. Another advantage of maximum likelihood 
classification algorithm is a high computational efficiency. It takes only 
slightly longer to complete classification then to complete eguai weighted 
filtering of an image. 

There are two major problems. First, is high misclassif ication r.'*"* at 
edges. This problem is created by the fact that the edge test, described in 
Section 3.B does not detect all edge regions. The problem of high 
misclassif ication rate at edges can be overcome by applying a better edge 
detector, for instance, maximum likelihood edge detector, developed in t16). 
Another problem is the imperfect extraction of target means. Partially, this 
problem is overcome by applying compensation test, described in Section 4.B. 
The extraction algorithm is by no means optimal, although it does estimate 
most target means well (see Section 5.B). Further research is needed to 
improve extraction technigues. Finally, classification error can be reduced 
by applying a simple post-processing algorithm. This algo*'ithm would select 
pixels or small groups of pixels that have been assigned to a target other 
than the majority of surrounding neighbors. These "isolated" pixels would 
then be re-assigned the same target level of as its neighbors. 
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A DATA COMPRESSION TECHNU)UE EOR 
SYNTHPrriC APERTURE RADAR IMAGES 

ABSTRACT 

A data compression technique is developed for Synthetic Ai>erture Radar 
(SAR) imagery. The technique is based on a SAR image model and is designed to 
preserve the local statistics in the image by an adaptive variable rate 
modification of block truncation coding (BTC). A data rate of approximately 
1.6 bits/pixei is achieved with the technique while maintaining the image 
quality and cultural (point like) targets. The algorithm requires no large 
data storage and is computationally simple. 


Intro* h>ct liin 


1 . 

The porjK)B»5 of thin jviper Is to *i«scribe rin efficient imaqe data compres- 
sion techni*iue which h.»u hn»»n sjwci finally ilcslqneil for Synthetic Afvfrture 
Radar (SAR) imaqes. SAR has become an important class of imaqin*j sensor for 
both civilian and military applications. As with other imaqinq systems there 
Is a need to transmit and store SAR imaqes, and thus, there has been consider- 
able interest in efficient codinq alqorithms for SAR (1,2,3). 

The aim of data compression is to minimize the data rate while 
maintaining the information contained in the signal using as simple an 
algorithm as possible. Thus a desirable property of an image encoding 
algorithm is fidelity, i.e. the reconstructed (received) image should preserve 
of all of the "important" features of the sensed image. For example, in some 
radar applications cultural features which appear in SAR images as small 
bright features are important and should be faithfully reproduced. In the 
geologic analysis of SAR images texture (4) is important and thus should be 
preserved. Image fidelity (quality) is a difficult quantity to measure (5) 
because it is application dependent. Image data compression algorithms are 
also compared on the basis of their compressing capability, i.e. the number of 
bits per image sample in the coded image. Implementation complexity is also 
an important consideration in evaluating data compression algorithms (6). The 
technique described below preserves important image features (e.g. cultural 
features) at about 1.6 bits/sample and is simple to implement. 

Standard compression techniques fall into two broad categories — 
predictive and transform coding. Predictive coding is performed in the 
spatial domain and attempts to remove the local redundancies in the image. 
Transform coding is performed by an energy preserving transformation of the 
image into another image so that the maximum information is placed into a 
minimum number of transform components (6). Many different transforms, e.g. 
Fourier, Cosine, Karhune-Loeve, have been used. Transform coding tends to be 
more complex than predictive coding. It will be shown in Section 2 that 
predictive coding is not possible on SAR images, and further, it will be 
argued that transform coding is not a viable alternative because of the low 
correlation observed in SAR images even though it has been tried (1). 

The technique developed here is a modification of the Block Truncation 
coding (BTC) algorithm developed in (7). BTC is suitable for SAR images 
because it preserves the local statistics of the image. In SAR these statis- 


tics «ir^ Important. Tn BT<' th*» ltna<jft i« dlvlii*?i! into Hmall hlockB («.«]. 4x4) 
of picture elements (pixels) end for each bloc)c a one bit quartizer is applied 
such tliat the block can be reconstructed with the moments (e.<j. mean and 
variance) preserved. Clearly, in addition to the bit mask (quantized block), 
supplementary information is needed. In the BTC the supplementary information 
is the sample moan, x, and standard deviation, s^« For example, a BTC system 
using 4x4 blocks of picture elements (pixels) (with 8 bits/pixel) and P bits 
to code x and 8,^ results in a 4 to 1 compression or 2 bits/pixel. Reducing 
the number of mean and standard deviation code bits as well as further coding 
of the bit mask can result in further compression [7,8]. 

Direct application of the BTC algorithm [7] (2 bits/pixel) to SAR images 
produced reconstructed images that were of suitable quality visually and 
preserved cultural features. The'^e results will be presented in Section 6. A 
further reduction in bit rate was achieved by observing that the local mean 
and variance are proportional in SAR images of homogenous areas, and thus, it 
is required to transmit only the mean. This modification results in a 
5.3 to 1 compression or 1.5 bits/pixe) (using the above example). 

This modification did produce reconstructed images of acceptable quality, 
however, significant contrast was lost for cultural features. This weakness 
was overcome by developing an adaptive BTC algorithm. The adaptive BTC 
algorithm sends only the mean if the local area (block) fits the standard 
radar model. For those blocks where the model does not fit, both the mean and 
standard deviation are transmitted. 

Using this adaptive approach, a 5 to 1 compression or about 1.6 
bits/pixel was achieved with the quality of the original 2 bits/pixel BTC 
algorithm. Now a variable number of bits per block is required for the 
adaptive BTC technique. However, this modification to the BTC algorithm does 
not significantly increase its complexity. 

A statistical model for SAR Images will be reviewed in Section 2. The 
original BTC technique will be discussed in Section 3. Tbe modifications to 
the BTC algorithm for SAR images will be described in the following section. 
The BTC, modified BTC (mean only) and adaptive BTC were implemented and tested 
using SEASAT-A SAR imagery. These results are presented in Section 6. The 
SAR image data compression described here is simple, produces reconstructed 
images of adequate quality for many applications and tends to preserve 
cultural futures. 
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, A fo r Raila r ImrttjwH 

2 . 1 Poi nt Sta t isti cs 

An Imaqinq radtir illuminates areas of the terrain within its field of 
view and records the value of the (X)wer returned from nonoverlapping resolu- 
tion cells on the ground. A resolution cell is typically made up of a large 
number of scatterers, and under some mild assumptions we can model the signal 
received by the radar (before detection) as a narrowband Gaussian random 
process. Then, with a square-law detector, the value of the received jjower P 
from a resolution cell has the probability density function (9, 10, 11] 


fp(p) 


w"’ exp(-p/y) 


for p > 0 


( 1 ) 


where E{ P) ■ u* 

In most imaging radars several independent measurements of the reflected 
power for each resolution cell are obtained and are averaged to form the image 
intensity value Y(t^,t 2 ) for the resolution cell with a spatial location 
(t^,t 2 )« Th® probability density function (pdf) of Y is the gamma distribu- 
tion (9, 10) of the form 


fy(y) - ^ 


N-1 


>-N 


(U/N) exp(-yN/u) 

r(N) 


y > 0 


( 2 ) 


where N is the number of independent measurements (or "looks"), 


N 


Y - E P and P(N) - (N - 1)1. 

” i-1 

The mean value u of the power reflected from a resolution cell is propor- 
tional to the radar reflectivity X of the resolution cell, and we can assume 
that u “ X without any loss of generality. Since the radar reflectivity 
changes from resolution cell to resolution cell, we can model the reflectivity 
as a random variable X (or a random process X(t^,t 2 ^) write (2) as a 
conditional pdf of the form 




y^“ ^ x~^ exp(-yN/x) 


r(N; N 


-N 


(3) 


Wo asrtim<! a Sworllnn type II (10] tar<|et mo«1ol to th»* Htatlstioal 

characterl«ticH of the echo on a por-pixel baain and then let the moan 
reflectivity X vary to ino<.1oI the SAR image of a large hetero<icneouB scene. 
With an api>roprlate change of variable we obtain the relationship between X 
and Y: 


Y(t 


1 '^ 2 > 


xct^.tj) zct^.tj) 
_ 


(4) 


where Z has a standard chi-square distribution with 2N degrees of freedom 
[10], and X and Z are statistically independent. 

Note that represents the speckle noise in SAR images. Here we 

have explicitly shown that X and Z are functions of position, however, for 
notatlonal convenience the spatial dependence will be dropped. It can be 
easily shown that for a given X 


E(Y/Xl -p 


(5) 


Y/X N 


( 6 ) 


For our purposes, we can see that if we are considering an "homogeneous" 
target (i.e., E(X) ■ u), then we can predict the variance given the mean. 
This observation is the basis for the modified BTC algorithm. 

2 . 2 Autocorrelation Properties of SAR Images 

The feasibility of using either predictive cr transform coding for SAR 
image can be discussed in terms of the image correlations properties. On a 
local level, i.e. inside a homogenous area, the model (equation 4) indicates 
that adjacent pixels will be uncorrelated. This is quite different comparsd 
to images collected with noncoherent sensors. For noncouerert sensors, e.g. 
LANDSAT or Aerial photographs, pixels in homogenous areas are highly 
correlated. The low correlation of adjacent pixels has been observed 
previously [12] and eliminates predictive coding from consideration. 

However, transform coding can can operate over larger regions and and 
thus the regional correlation properties of SAR images needs to be 
considered. Figure 1 presents a typical autocorrelation* in the row and 


column directions of a f?AR Imaiio of terrain (s«o Flqur« 4 for tin* 

Note that thin aatocorrnlation functions decays very rapidly. The prnsttnee of 
some correlation indlcaten that transform coding is (Kissible, howovur the 
rapid decay implies that the size of the transform window muct he large, thus 
greatly increasing the memory and computational requirements of the 
compression algorithm. A sophisticated transform coding algorithm for SAR 
images has been reported [1). The technique described in (1] uses row/column 
deletation (resampling) and transform coding; however, the algorithm is 
computationally intensive and requires substantial amounts of memory. 

Reexamining the SAR image model (equation (4)) we notice that the mean 
and variance of each local area (within a larger homoneneous region) are 
redundant, and thus, a coding algorithm which preserves these image attributes 
would be suitable for SAR. In the next section, the BTC algorithm which does 
preserve these features is discussed, and it is modified to fit the above SAR 
image models. 


3. Review of Bloc)c Truncation Coding 

Let y and S^ denote the sample mean and standard deviation of a block 
(e.g. 4x4) of pixels in a SAR image. That is 


and 


Sy -{y^ - (y? (8) 

with 

" i-i ^ 

where 

m ■ number of pixels in the block 

and 

■ pixel intensity 


*7tie image contained 512 x 512 pixels and the autocorrelation function was 
obtained using FFT techniques. 


I ^ 
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In BTC a on*» bit qiiantl/,«r 1 h foi M.ich pixnl in tb»* block with tr>»' 

quantizin'! threahold H«t at y. That 1«, If > V *^*’**" *^*'«*f plx^l location in 
coded with a 1, otherwise with a 0. A bit mask in thus forned. This bit mank 
alonq with y and Sy tranBmitted/stored. 

At the receiver, a level A is asBigned to a point if that pixel location 
within the block contained a 0 and a level B if a 1 was contained. The lovelB 
A and B are selected to preserve the moments of the block. These levels can 
be simply found as [7, 13] 


A - y - S, 


and 


\/!ZI 

Y m-q 


where 


B ■ y • s, 


( 10 ) 


( 11 ) 


q m number of one's in the received block. 

For a 4x4 block and using 8 bits to code y and Sy results in 4 to 1 
compression ratio or 2 bits/pixel. 

4 . Modified Block Truncation Coding 

Based on the SAR image model, we can predict the standard deviation for a 
block based on the sample mean. Let the predicted standard deviation be 


where 


a - 


( 12 ) 


N - number of looks for the SAR. 

The BTC technique described in Section 3 is then applied at the source. Now 
only the sample mean and the bit mask are transmitted. At the receiver, the 
levels A & B are reconstructed using' 


(1 Vi 


N(m-q) 


and 


B - y (1 


-/ 


m-q 

Nq 


(13) 

(14) 


For a 4x4 block and using 8 bits to code the sample mean a compression ratio 
of 5.3 to 1 or a data rate of 1.5 bits/pixel is obtained. 


It wrtH todiul thfit thlH t**rhnl<|iu* prfxliicwl rncotta true ted of 

homogenaouM <irorts almont identical to t»u» 2 bit/pixel technique. However, 
contrest was lost on cultural features relative to the Initial HTC methfxl. 

Hie reason for this la obvious. Regions containing cultural features do not 
fit the SAR model of equation (4). To overcome this weakness an adaptive 
variable bit rate BTC was developed. 

5. Adaptive Variable Bit Rate Block Truncation Oadlng 

The goal of the adaptive BTC technique ia to use the modified BTC for 
those image blocks where it is appropriate, i.e. where the model fits, and to 
use the original BTC algorithm otherwise. A simple test based on the pre- 
dicted and sample variances was developed to Indicate if the data from an 

image block fits the SAR model. Specifically, if k*S ^ > a ^ (where k is a 

y P 

constant) then the model does not fit the data. That is, if the observed 
(sample) variance is "too much" greater than the predicted variance then it 
would be expected that the model does not provide an adequate description for 
the data. The proportionality constant, k, must be selected such that the 
probability of rejecting the model when the model is valid, Pp, is small, 
i.e.. 


P(k S ^ > o ^ [model is valid) 

\j n ' 


(15) 


From the theory of confidence intervals [14] we know that we can find a B such 
that 


m S 


m S 


P(- 


> B) » P(o. 


-) - P, 


(16) 


The above equation indicates that the probability that m S ^/B (here k ■ -5^) is 

9 y p 

greater than the unknown parameter o is p when the SAR mocel fits the 

, P ^ 

da'.a. Ttie predicted variance o ^ is calculated using equation (12). Ttiere- 

P 

fore, the SAR model is not appropriate if 


or 




mN 

B 


(17) 


( 18 ) 


whpr** Is Hi*lect.»*.l to forro h*. „n»all, n.q. 10"^. Note th.it y^/S ^ is 

an estimate for the number of looks usmi In the SAR processing based on the in 
pixels in tiie block. This ratio will be referred to as the local number of 
looks . 

Assuming that the pixel intern ties are Gaussian distributed (which is 

true for large N) the random variable n s ^/o ^ a distribution with m-1 

y P 

degrees of freedom. Thus, can be estimated. For example, let m ■ 16 (i.e. 
4x4 blocks) and N*»4 then for a Pp ■ .005, 0 ■ 32.8, k a .5 and mN/B * 2. 
this case we would expect that the sample standard deviations would be 
needlessly transmitted only for 1 in every 200 blocks. 

The adaptive BTC algorithm for each block is implemented as shown in 
Figure 2. At the transmitter, the sample mean and standard deviation are 
calculated and the bit mask is formed as specified in the original BTC 
algorithm. An estimate for the number of looks, y /S^ , is calculated next 
and compared to a threshold mN/B. If the local number of looks is less than 
the threshold, then a flag is set on and a data block Is sent which contains 
the bit mask, sample mean, sample standard deviation, and the flag bit 
indicating the presence of the standard deviation. If the local number of 
looks exceeds the threshold then it is highly probable that the SAR model is 
valid and thus the standard deviation present flag is set off and a block is 
sent which contains only the bit mask, sample mean, and the flag. At the 
receiver the flag is tested if it is on the original BTC reconstruction 
algorithm is applied, i.e. equations (10 and 11), otherwise the modified BTC 
algorithm (equations 13 and 14) are used to reconstruct the image. 

The adaptive BTC algorithm described above automatically adds 1/16 
bit/pixel overhead, the standard deviation present flag. It was found that 
for SEASAT-A SAR images that the quality of the original BTC technique was 
maintained using the adaptive approach at approximately 1.6 bits/pixel. 

The BTC algorithm produces reconstructed images which have a blocky 
appearance (13) when displayed with magnification. This characteristic also 
evident to when the BTC algorithm is applied to SAR images. However, the 
image reconstruction algorithm in the adaptive BTC technique can be modified 
to remove this blocky appearance. 

When the SAR image model, e.g. (2), is valid and only the mean is 
transmitted we know that a pixels marked with a 1(0) can be modelled by a 
conditional probability dersicy function, i.e., the p.d.f. given in equation 2 


condi tlorie.l on th»* ov».*nt Huit thn s.^mple is above (below) the mean. Synthetic 
samplinq (usiny psiMido random numbers) can be used in the reconstruction to 
map pixels marked with a 1(0) into a qray level bascsd on the appropriate 
conditional p.d.f. As will bo shown in the followiny section using pseudo 
random numbers in the reconstruction algorithm does effectively remove the 
blockness. Unfortunately, generation of pseudo random numbers is 
computationally intensive, thus this modification increases the computational 
complexity of the algorithm. In some applications this refinement will not be 
required . 

In the following sections the four compression algorithms, original, 
modified and adaptive BTC and adaptive BTC with pseudo random reconstruction 
will be compared. 

6. Results 

The compression algorithms, BTC, modified BTC, and adaptive BTC and 
adaptive BTC with pseudo random reconstruction have been applied to SEASAT-A 
SAR imagery. The purpose of this section is to discuss their performance. As 
mentioned previously, image quality is difficult to quantify. Here the 
performance evaluation is based on two criteria: 1) the faithful reproduction 

of cultural features, and 2) the general appearance of the reconstructed image 
relative to the original. 

The SEASAT-A SAR imagery used here had a resolution of 25 x 25 m with 
N “ 4. Each pixel intensity was represented by 8 bits (0-255 grey levels). 
(For more details about the sensor see [15].) For all the compression 
algorithms described here, the sample mean was coded using 8 bits/block. 

Also, a 4 X 4 pixel block was used in all cases. In the original BTC 
algorithm the standard deviation was also coded using 8 bits/block resulting 
in a 2 bits/pixel data rate. 

The modified BTC algorithm resulted in a 1.5 bit/pixel data rate. In the 
adaptive BTC algorithm, the standard deviation was coded using 7 bits/block, 
then allowing for the one overhead bit results in a maximum of 32 bits/block. 
Thus, in the adaptive BTC technique each block is coded into 25 or 32 bits 

depending on the statistics of the pixel intensities of the block. A 

threshold of 2.0 was used in all cases. The data rate of the adaptive BTC 

algorithm is variable. However, in most cases it was about 1.6 bits/pixel. 
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Thf r<?s|}ona«} of the firat three alqorlthms to a "j)oint” like tcirqet ia 
shown one-illmeiisionrtl ly in Fiqure 3. Fiyure 3a representa an inte.isity 
profile of 100 pixels from a SEASAT-SAR image. There is one bright feature in 
the center of this profile. The target-to-background contrast in Figure 3 is 
about 9dB. Figure 3b is the reconstructed profile using the original BTC 
technique. The target is still quite evident, the average background level 
has remained the same as expected and the target level has been reduced. The 
target-to-background contrast is about 7.5 dB, a loss of 1 .5 dB. The result 
of tne modified BTC algorithm is shown in Figure 3c. In this case, the 
target-to-background loss is about 3.4 dB. This loss is not considered 
acceptable. The adaptive BTC algorithm. Figure 3d, restored the profile to 
that given by the original BTC at a small cost in data rate from 1.5 
bits/pixel to 1.58 bits/pixel. There is little difference between the 
original and adaptive BTC results (Figures 3b and 3d), this observation is 
true for all the results presented here. 

A scene containing a variety of terrain features is shown in Figure 4. 
This scene is composed of 512 x 512 pixels. The three compressed images 
favorably compare to the original in terms of reproducing the terrain 
features. However, there is some difference in scale of the texture in the 
homogeneous regions caused by the block nature of the block coding technique. 
These differences are more easily seen in the agricultural scene shown in 
Figure 5. The texture patterns in the reconstructed images appear as speckle 
patterns. This might be attributed to the compression algorithm’s two level 
quantization of the pixels in each block. So the algorithm is mapping all the 
up fades to one value and the down fades to another. 

The upper left corner of the agricultural scene is shown magnified in 
Figure 6. At this scale the "blockyness" property (13) of the algorithm is 
clearly illustrated. However, the shape of most cultural features is pre- 
served. Specifically, the features identified as 1 , 2, and 3 on the original 
SEASAT-A image. Figure 6a, are preserved in shape in all three (6b, c and d) 
reconstructed images. Note the loss of contrast for the "point" like target 
(feature #2) between the original BTC (Figure 6b) and the modified BTC (Figure 
5c) reconstructed images. The result of the adaptive BTC algorithm using 
pseudo-random reconstruction is shown in Figure 6e. As expected this 
refinement reduced theblocky appearance of the reconstructed image. Itie 
properties of the BTC algorithm are also evident in Figure 7. Hiis scene 


) 


contdins ii watrtr Jxjdy, a dam, and a fjowor tratismlHsion line (the row of brlqht 
points near the bottom of the scone). Again, the shaiie of these features Is 
preserved and the modified BTC algorithm shows a loss of contrast for the 
point targets. 

7. Conclusions 

A data compression technique has been developed for SAR Images. The 
method was tested on SEASAT-A SAR data and found to produce Images with a 
suitable quality for a variety of applications. Itie algorithm developed here 
Is an extension of the BTC technique developed in [7]. The specific statis- 
tical properties of SAR data were used to Improve the data compression 
ratio. A compression ratio of 5 to 1 (data rate of 1.6 bits/pixel) was 
obtained. Further minor reductions in data rate might be possible by reducing 
the number of bits used to represent the sample mean and standard deviation as 
suggested in t18]. The benefit of such a reduction would have to be 
considered based on the application of the sensor. Also, further study is 
needed to evaluate the effect of changing the number of looks of the SAR on 
the quality of the reconstructed images. 

The technique presented here requires no large dati storage as opposed to 
transform coding methods [1], and is computationally simple so that a single 
chip implementation is possible [7]. As the SAR image formation moves closer 
to a real time operation and is thus performed on the sensor, data compression 
techniques as the one presented here will provide the system designer with 
additional trade-offs for transmitting and storing the data. 
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Figure 2a. Adaptive BTC Algorithm 
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Figure 6. 

Results of Coding for an Agricultural Scene: Magnified 
a) original, b) BTC, c) modified BTC, d) adaptive BTC 
and e) adaptive BTC with pseudo-random reconstruction 
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Figure 7. 

Results of Coding for General Scene 
a) original, b) BTC, c) modified BTC, d) adaptive BTC 
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■\nSTWAi’T 


An optimal frecpitMicy domain te'Xtural rdijf dft»‘i-tion f ilt»*r i‘i di vi'l<-)pt*d 
and its [ler f ormanco evaliiat«*d. For th»» qivcn ono-d im**n« i ona 1 t^xtnn- mod*-l, 
and filter bandwidth, the filter maximizes the amount of output ima*j<* enerqy 
placed within a qlven resolution Interval centered on the textural I'dqe. 

Filter derivation is based on relatinq textural edije detection to tonal edqi* 
detection via the complex lowpass equivalent r<presentation of narrowband 
bandpass signals and systems. The filter is specified in terms of translated- 
in-frequency prolate spheriodal wave functions. performance is evaluated 
using the asymptotic approximation version of the filter. This evaluation 
demonstrates satisfactory filter performance for ideal and non-ideal tex- 
tures. In addition, the filter can be adjusted to detect textural edges in 
noisy images at the expense of edge resolution. 


* • 


iNTi<oni'<"ri'>N A'li) o'.n'KV I r.v 


lift (»i*t 1 un in .til inport.int- first stf|> ui •*xt r»i«'f. inij uif orn.i * i nn t r )m 
an imaifr*. Many i1ft»*ction srli»>n'‘H h.ivi* boon nn[)l oy<‘H to <‘nh.inri' th<* 

boiinda r i <*s Ix-'twpon roifirinn of iiiff»*rr*nt avf*rarf»* 'fray tono . Th»*Ho tonal •*d<fo 
detectors are i naileipi.ite wher regions in an image are cha racter i zed by similar 
average gray tone, but different textural features. 

A textural ed<je detection filter is presented in this pafier which is 
optimal in the sense that, for the given model, a maximum amount of output 
image energy is placed within a given resolution Interval width and a given 
filter bandwidth. The resolution interval is centered on the textural edge in 
the input image. The filter is derived in the frequency domain, and is easily 
implemented on a digital computer using Fast Fourier Transform (FFT) tech- 
niques. 

The optimum textural edge detection filter is developed by treating the 
textural edge as a bandpass extension of a tonal edge. Hence, the optimum 
tonal edge detector derived by Shanmugan, Dickey and Green [1] (correspondence 
by Lunscher (2)), is related to the textural edge detection case via the 
complex lowpass equivalent representation of signals and systems. It should 
be pointed out that the development is carried out in one-dJ mension. However, 
symmetries required for extension to two-dimensions are retained. 

Section II presents a brief review of the optimum tonal edge detector. 

The textural model used in the development of the optimum textural edge detec- 
tor is then introduced in Section III. The mathematical form of the optimum 
textural edge detection filter and some one-dimensional examples are presented 
in Section IV. Concluding remarks are given in Section V. 


II. R!;vii'v OK Tin; optimum 'idnai. i'ihii: oktkotoi? 


Ths piir|H)s»* of this s»fction is to bri**fly r^-vii'w tlv optiniiim ton.il 
dpt**f'tor <1«Tiv*»fi liy Sh.inmiKj.in , ft al., Ml. For .i nivfn hrimlwilth , th** 

optimiun ton.il fdijf detector places a maxitnum amount of output imaqe enerqy 
within a qiven resolution Interval lenqth in the vicinity of tonal edqes. The 
tonal cdqe detector la inaensitive to textural edqes where the averaqe qray 
levels of the different textural reqions are equal. 

The derivation of the optimum tonal edqe detector is based on represent- 
inq the filter output (for a step edqe input) in terms of prolate spherlodal 
wave functions (for the derivation, see (1), (21). The exact one-dimensional 
form of the filter transfer function is qiven in Shanmugan, et al., (1) as 


H ( u)) I 

STEP , E ' 


(j)i/2n), |u)| < n 

0 elsewhere 


(1 ) 


Tl I 

where c ■ and 4 ;^ is the first order prolate spheriodal wave function. (The 
subscript STEP,E in Equation (1) denotes the Exact form of t*ie STEP edqe 
detector). For any given values of spatial bandwidth, £1, and resolution 
interval length, I, the transfer function in Equation (1) places the maximum 
amount of energy in I. The filter is difficult to implement in this form, 
because the values of 4 ;^ are tabulated. Application of approximations by 
Slepian and Streifer (1), yield the asymptotic approximation of the filter, 
which is in closed form, hence easy to implement. The resulting expression is 


^ 2 C(i) 

H (o)) - H_„^^(w) = K, cj exp( -) 

STEP,E STEP 1 


( 2 ) 


Combining the constants that appear in the argument of the exponent, and 
dropping the gain factor, , yields 


? 2 
, . , -Ku) . 2 -Ku) 

(w) = (ij(we ) = w e 

STEP 


(3) 


It should be noted that the parameters I and fl can no longer be indepen- 
dently specified. 
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;'hoK'<' tif K rtt'trt flu* li.in>lwi>)Lli of ttu* ,»nil 0I50 r<?sol'itlon 

intotvril lnn<Mh. As K i ncr«*.ist*3 , rusoliiflon Interval .siz«' incr»»<»Hos, r*n <1 

-/i <1 til ilncrtfoses . Not«> that ovon thomjh the asymptotic approxima- 
tion to the optimum transfer function is not strictly band 1 imi ted , Hj,^p.p((u) is 
effectively zero for spatial frequencies above a certain value, dependinq on 
the choice of K. The asymi>totic approximation will be used in the remainder 
of the deve lojxnent , 


III. n'XTIM^M, MMDn, 


Orif inInT«Mit ill f f ICO I t.y with textorol procoswintj is tho fact thst '’o 
ainql»' 'r st" mo.lcl »»xistH for chsr/ict«»rizin<j texture in Inaqes. The model 
used here in the (1»-ve lojiment of the optimum textural edqe detector capitalizes 
on the relationshi(> between texture and spatial frequency by representinq each 
texture as a sinusoid of different spatial frequency (i.e., fine textures 
contain qreater concentra t ions of enerqy at hiqher spatial frequencies than 
coarser textures lio) [3], (4), 151, Ih), (71, (8), [9]. 

In qeneral, a class of one-dimensional imaqes with n textures can be 
defined as 

q(x) - A(x) cos(u)^x + 0(x)) (4) 


where 


A(x) = a( 1 4- a(x) ) |a(x) I < 


(5a) 


and 


x 

0(x) - b / Bv X) dX (5b) 

GO 

The functions a(x) and 0(x) are t indom processes, represents the ith tex- 
ture, a and b are constants, and x is the spatial variable. Note that q(x) is 
allowed to be neqative. This can be viewed is subtracting off the mean level 
from in image, thus allowing neqative brightness or gray level. In this 
model, a(x) represents average gray level, and 0(x) represents the variation 
of spatial frequency within a texture. In other words, the envelope of q(x) 
can be thought of as the average gray level variation, while the underlying 
texture is represented by each different where the random change of tex- 
ture for a given is controlled by 0(x). Note that if time were the inde- 
pendent varial le, q(x) would he a double sidenand plus large carrier modulated 
waveform, with simultaneous frequency modulation. 

An ideal textuie is represented in this model by a sinusoid with constant 
spatial frequency and constant amplitude. Hence, a transition between two 


- - 


t,»'Xt<ir*‘S Clin r<^pr«;M«*nt<'d by .1 |hii<‘ siiiiiHoid .it otn* ••,|m'i.i 1 t r<M;in-ticy 

followed by ti piire sinusoid mother npiitidl frepuenry. For t be ide.il two 

texture case let 
A(x) - t 
0(x) - 0 

- OD < X < °° (infinite size) 

Thus, an iiTifiqe with two ideal textures and a textural edqe at x ■ 0 is repre- 
sented mathematically as 

f(x) ■ cos(ii)j^x), -® < X < « (6) 

where 

i » 1 tor X < 0 and 
i - 2 for X > 0. 

The optimum textural edge detector is derived using the ideal, two texture 
image, f(x). 
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!v. Mi'Tr-inM Trv.TiiMAi. nir.i-: iirn-rToi' in'sm/rs mjd m'KotJfiMjiT 

This •.iTtioM pr'.'seii's th<- ic.il form of ih** optinom •• f>i*‘ 

i on filter .uiil ifisrusnes th€> i>er f orm.uice of the filter for scvit.i 1 
(fiffercnt rl.issi*R of input Imiqes. The derivstion is only briefly skctchfrf 
here, the ifet.iils <ire qiven in Townsend [10). 

For a tw<i texture input imaqe with one texture represented by a sinusoid 
with frequency , and the other texture represented by a sinusoid with fre- 
quency ojji the transfer function of the optimum tonal edge detector is given 
by 

H (oil “ H ((d) + (7) 

OPT 1 2 


where 


H ((D) » H 

( (D - 

(D ) + H 

3 

+ 

3 

1 

STEP 

1 STEP 

1 

H^((d) = H 

( ID - 

(1) ) + H 

( (D + (D„ ) 


STEP 

2 STEP 

2 


and 


STEP 


-K(d 


(3) 


It is clear from Equations (7), (8), and (3), that the optimum textural edge 
detector is the sum of the responses of two bandpass "sub" filters, ((d) and 
H 2 ((d). Each "sub" filter is a translated-in-f requency version of the optimum 
tonal edge detector, h<-,^gp((D), discussed in Section II. Note that Hg,pgp((D) is 
translated to each of the two textural frequencies. 


The opti^’im textural edge detector is derived by recognizing that the 
two-ideal-texture input imaqe, f(x), given in Section III can be expressed as 
the sum of two truncated sinusoids, one at frequency (d.j , defined for -» < x < 

0 and the other at frequency id^, defined for 0 < x < +". But each of these 
two truncated sinusoids are bandpass at frequencies (d^ and (D 2 respectively. 
Each truncated sinusoid has a step function for its complex lowpass equivalent 
Ml). Because Hs.j.pp(t*J) is optimized for detecting step type edges, a bandpass 
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viTiion of ) I'lMi 1 1 ‘ ri'c) nn f r»*ijiii‘nry ^ is n[)tini)n tor i tui t'i>- 

i1i scont i mii ty ( modiil o t»-,f step function), in tli«* tr nnc-i ted siniisoir) rtt f rrspien- 
cy MO). iiml.irly, <i bandpass version of ( •.>) translate,) in frequency 

to ijy optinun for detectinq tbe <i isconti mii ♦ tlie truncated sinusoid at 

frequency The sun of the out(>uts of ttiese two indpass filters produces 

the optimized output. A bloclc diaqram of the filter structure for tlie two 
texture case is shown in Fiqure 1. 

A qualitative discussion is presented here to qain inslqht into how the 
filter worlcs. Fiqure 2 presents an example of the optimum textural edqe 
detector in the frequency domain. Note from the fiqure that the response at 
and Wj (the spatial freciuencies representinq the two ideal textures) is 
zero. Hence, does not respond to any input which has spectral enerqy 

only at these two frequencies. Therefore, the response to an input represent- 
inq eicher pure texture (in steady state) is zero. The textural edqe is 
characterized by a transition from one texture to the otlier. The Fourier 
transform of this boundary contains spectral enerqy at frequencies other than 
and <i> 2 * In particular, there is enerqy in the passband portions of 
Hqpt( » therefore filter response near the textural edqe is non-zero result- 
inq in a larqe amount of output imaqe enerqy in the vicinity of the textural 
edqe . 


The Fourier transform of the entire input imaqe is qiven by 
F(oj) = F^ (oj) + ^2^“’ 

where F,j(aj) amd F 2 (w) are the Fourier transforms of the truncated textures 
represent* J by sinusoids at w,, and W 2 respectively. Multiplication of F(uj) 
with riQp,p(a)) yields the transform of the output, G(o)), i.e., 

G(w) = F(oo) (w) (10) 

OPT 

but this is ec)uivalent to 

G(o)) = iF^(w) + F2 (co)] (H^(o)) + H2(w)) 

= F (w) H (o)) + F (w) H ((!))+ F (u) H (oj) + F (w) H,(u) (11) 

11 12 2 1 2 2 
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h\it 


(.*i) H (oj) = 0 (12) 

1 2 

cirid 

F 2 ( u ) H ^ ( u) ) = 0 (13) 

Substitvition of Equations (12) and (13) into EJquation (11) yields 
G(w) F^(u) H^(uj) + F2 (w) ^^2^***^ 

= (w) + G^lw) 114) 


Hence, 

q(x)=g^(x)+g^(x) (15) 

Equations (12) and (13) are true because of the spectral separation be- 
tween t)>e two sets of bandpass inputs and systems. In non-ideal texture 
cases, there can be considerable spectral overlap h>etween the Fourier trans- 
forms of the textures. The spectral overlap can cause non-:tero response of a 
system, H^(u)), for example, to a textxure not centered at , F 2 (to) for exam- 
ple. This could also occur if the bandpass bandwidth of H^(uj) is wide enough 
to pass a significant amount of energy due to F 2 (uj). 

Choosing the exponential parameter, K, such that the bandpass bandwidths 
of M^(u)) and H 2 (o)) are wider than the spatial frequency separation between 
and ti >2 results in non-zero response to the two textures. There is improved 
resolution at the expense of an increase in the "bacitground" level in the 
output image, thus decreasing edge visibility. The "background" refers to the 
out-of-resolution-interval gray level. Edge visibility describes the dif- 
ference in gray level between the in-resolution-interval and out-of-resolu- 
tion-inter val (background) portions of the output image. The spatial fre- 
quency separation of the textures affects the fjerformance of the filter, i.e., 
the greater the separation, the better the i>>rformance . 
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f 


!• Wii •■.i'.DWti lii 'ih.ri- ■: i.ri , , (1) t)wil tin* < )pi. i mii”i toti.il 

tor f'dil 1 !v usi*i) to i>nh.itii*t* t()n,il in imaip’s t:or rnpt.t'ii by adilitivo whitn 

oiUissi.tit Th<* *..in** •■'n*ory •tpplit’s to th<’ optimum tnxt'irfil o.lqft iU»teo- 

tor . Tin* I'xi xTn**M t 1 . 1 1 [ .ir.«m<*t**r, K, can hr; cliost»n to decrease the bandwidth of 
the "sub" filters to deort*ase the effects of the noise. The price jiaid for 
this is an increase in the resolution interval lenqth (10). The benefits of 
increased edqe visibility may more than offset the decrease in resolution. 

Fiqure 3 shows the result of implementinq the filter on a diqital compu- 
ter. Displayed are the input and output imaqes (one-dimensional) of the 
optimum textural edqe detection filter for an input with tv«ro ideal textures 
(one textural edqe). The textural edqe is clearly marked in the output imaqe. 


The transfer function, Hf^p,p((i)), can be generalized to n textures by 
simply adding more translated-in-f requency versions of ^s-rEp^w). Denote the 
generalized, n texture transfer function as hQp,j, ^^( u>) , defined as 


H (ui) 

0PT,n 


n 

V 

L 

i-=1 


H, (w) 


(16) 


where 


H. (uj) 
1 


(in- 

step 


U). ) 
1 


H (u) 
STEP 


0 ). ) 
1 


(17) 


and represents the frequency of the ith texture. Each of the n filters 
respond to transient energy where textural transitions occur but null out 
response to the ith texture in steady state. An example of a one-dimensional 
output imaqe for an input imaqe containing four ideal textures with three 
textural edges is shown in Fiqure 4. The normalized frequencies of the four 
different textures in the fiqure are .04n, .06 ti, .OSn, and .Iti, with each 
texture occurring once in the input imaqe. 

It should be pointed out that although each of the "sub" filters (i.e., 
H^(uj), H 2 (uj), •••) are narrowband bandpass about the respective textural 
frequencies, the overall system bandwidth and imaqe bandwidth are about equal, 
as shown in Fiqure 5. The total textural edqe detector bandwidth, BW, is 
written in terms of the tonal edge detector bandwidth as follows: 
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igure 4 Magnitude of the response of the textural edge detection 
filter due to an input image with four ideal textures and 
three textural edges. The normalized spatial frequencies 
of the four textures are .04 tt, .06-, .08:^, and .1 -:t. 
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The most qeneral case of the model used in this development is one in 
which each of the s))atial frequencies representi nij ttie different text\ires iri 
the imaqe are allowed to rardomly deviate about some aver.jqe frequency. This 
complication is introduced to allow for some of the irreqularity of a real 
texture. A one-dimensional exampla in which both the amplitude and spatial 
frequency vary in proportion to independent random processes is shown in 
Figure 6. In this example, the average normalized spatial frequencies repre- 
senting the two textures are .04n and .In respectively. In terms of the 
general model presented in Section III, a(x) and fl(x) are independent Gaussian 
noise processes, with unit variance. The bandwidths of the amplitude noise 
and frequency noise processes are .008n and .006n respectively. Note that the 
filter adequately marks the two textural edges in the image, but also responds 
to regions within each texture where the spatial frequency changes. Decreas 
inq the bandwidth of the noise modulating the frequency causes the spectral 
separation of the textures in the input image to increase. This results in 
improved performance of the filter at distinguishing textural edges from 
frequency deviations within a texture. 
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of the optLiTun textural edge detector response 
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A I r 'i'Jin.Jln . 1 <’’ i on li.ts li> imi I.‘7>- 1 i i| ,«m) 

which, (or tiu* ‘livini mn'lcl l).inhwi htii, .1 niximin .inonrit of 

<*nf*r<jy within a Bf>cci f if*<i resolution interval near the tex'utal er)i|e. 

Til'* textural <let»*ctor was deriveil by rnlatinh textural i><i«te .letertion to 

tonal edqe det«Ttion vi,i comjilex lowf>ass ectui valent representation. Hence, 
the optimum textural e(l<|o detector was found to lie a sum of translated-in-f re- 
quency versions of the optimum tonal edqe detector. Tliis form allows the 
filter to be adapted to multi textural imaqes. In addition, examples were 
presented which show the filter's insensitivity to tonal features in an 
image. The filter is adjustable; resolution can be traded for edge visibility 
in the case wher«. the input image has been corrupted by noise. 

The qualitative and complex nature of texture suggests that a totally 
general approach to modeling and classifying texture may never lie found. It 
has been an objective in this investigation to develop a filter which opti- 
mizes a certain criteria relating to textural edge detection. But, as always, 
simplifications and assumptions wore made indicating the need for further 
research. The model used in this development represented texture in terms of 
spatial frequency, and gray tone in terms of amplitude. One example of fur- 
ther research might be to base the development on a muie complex model which 
incorporates a statistical description of texture. in addition, further work 
is needed in extension of the one-dimensional filter to two-dimensions. 

This work has provided an approach to textural edge detection which can 
be implemented on digital hardware using the FFT. With the increased size and 
availability of digital computing facilities at a decreased cost, digital 
image processing methods will become more popular in the future. 
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